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1  .  INTRODUCTION 


1 . 1  Summary 


It  is  known  that  large  space  structures  will  be  subjected  to 
thermomechanical  loadings  and  environmental  conditions  which  are 
likely  to  degrade  the  constitutive  properties  of  the  structural 
materials,  thus  leading  to  possible  failure  of  these  vehicles. 
Therefore,  it  is  desirable  to  develop  new  analytical  models  which 
are  capable  of  accounting  for  these  degraded  properties  so  that 
design  procedures  can  be  improved.  There  are  three  important 
aspects  of  such  an  effort:  selection  and  development  of 
constitutive  models  for  degrading  materials  which  are  applicable 
to  large  space  structures,  construction  of  analytic  models  for 
predicting  the  dynamic  response  of  these  structures,  and 
experimentation  to  determine  the  precise  nature  of  the  material 
parameters  to  be  utilized  in  the  analytical  model.  These  three 
components  of  the  research  must  be  tied  together  into  a  single 
concise  effort  in  order  to  obtain  a  useful  model. 

This  research  proj  ect  i s  £  three  year  effort  to  develop  an 
analytic  model  capable  of  pr edl ct ing  the  response  of  space 
structures  with  degrading  mat  ar lal  properties  under  quasi-static 
as  well  as  dynamic  eye  lie  thermomechanical  loading  conditions. 
This  report  details  the  research  completed  during  the  second  year 
of  AFOSR  contract  no.  F4962(H83^0067 . 
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Statement  of  Work 


model  is  being  developed  for  predicting  the 
thermomechanical  response  of  large  space  structures  to  cyclic 
transient  temperature  loading  conditions.  The  research  is  being 
conducted  in  the  following  stages: 


1)  selection  and  specialization  of  thermomechanical 
constitutive  equations  to  be  utilized  in  the  analysis  of  large 
space  structures; 

2)  construction  (where  necessary)  of  coupled  energy  balance 
equations  (modified  Fourier  heat  conduction  equations)  applicable 
to  the  constitutive  models  selected  in  item  1); 

3)  casting  (where  necessary)  the  resulting  field  laws  into 
coupled  and  uncoupled  variational  principles  suitable  for  use 
with  the  finite  element  method; 

4)  finite  element  discretization  of  the  variational 
principles  for  several  element  types; 

5)  experimentation  to  determine  material  properties  to  be 
utilized  in  the  constitutive  models;  and 

6)  parametric  studies  of  the  quasi-static  and  dynamic 
response  of  large  space  structures  undergoing  thermomechanical  ly 
and  environmentally  degraded  material  properties. 

The  experimental  effort  (discussed  in  5)  is  ^el^hg  supported 
by  DOD  equipment  grant  no.  841542.  The  total  research  effort 
outlined  above  spans  a  three  year  period.  The  following  section 
details  results  obtained  during  the  second  year. 
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RESEARCH  DURING  SECOND  YEAR 


2 . 1  Summary  of  Completed  Research 

The  following  tasks  have  been  completed  during  the  second 
year  of  research: 

1)  development  of  generalized  constitutive  equations  for 
metal  matrix  composites  with  distributed  damage; 

2)  experimentation  to  determine  material  parameters  for  the 
model  developed  in  item  1 )  ; 

3)  completion  of  algorithmic  development  for  space 
structures  with  damage  induced  and  spacially  variable  stiffness 
loss ; 

4)  completion  of  parametric  studies  for  graphite/epoxy 
composite  space  structures  using  item  3); 

5)  completion  of  algorithmic  development  for  viscoplastic 
space  structures  with  thermomechanically  induced  heating; 

6)  completion  of  parametric  studies  for  aluminum  space 
structures  using  item  5);  and 

7)  completion  of  development  of  bounding  techniques  for 
hysteret ical ly  induced  temperature  ri3e  in  t her mo v i scop  1  as t  i  c 
space  structures. 

In  addition,  the  following  tasks  are  well  underway  at  the 
completion  of  the  second  year  of  research: 

1)  development  of  an  analytic  method  for  modeling  beam-like 
structural  components  with  damage  induced  stiffness  loss; 

2)  development  of  a  finite  element  model  for  beam-like 
space  structures  with  spatially  degrading  material  properties  and 
subjected  to  solar  flux  heating  and  radiation  boundary 
conditions;  and 

3)  development  of  a  model  for  predicting  the  structural 
response  of  space  structures  with  load  induced  damage  which 
causes  changes  in  structural  damping. 

The  items  briefly  outlined  above  will  be  detailed  further  in 
the  following  sections. 


2 . 2  Development  of  Constitutive  Equations 


As  stated  in  the  original  three  year  proposal  [1],  one 
objective  of  the  research  was  to  identify  constitutive  equations 
which  account  for  material  property  degradation  in  polymeric 
composites,  metal  matrix  composites,  and  high  strength  metal 
alloys.  Constitutive  equations  for  polymeric  composites  and  high 
strength  metal  alloys  were  found  to  exist  in  the  research 
literature,  and  candidate  models  were  chosen  and  reported  in  the 
first  annual  report  [2].  However,  no  appropriate  models  for 
metal  matrix  composites  were  found  in  the  literature.  Therefore, 
it  was  decided  to  develop  a  constitutive  model  for  metal  matrix 
composites  which  would  be  applicable  to  space  structures!!  3],  The 
current  state  of  this  model  development  will  be  detailed  in  this 
section. 


2.2.1.  Generalized  Model  Development  for  Metal  Matrix  Composites 

Because  metal  matrix  composites  are  expected  to  be  utilized 
commonly  in  space  structural  applications,  it  was  felt  that  some 
constitutive  model  development  was  warranted  for  this  class  of 
materials.  The  distinguishing  feature  of  metal  matrix  composites 
is  the  substantial  inelastic  (either  elastic^plastic  or 
viscoplastic)  nonlinearity  which  occurs  in  the  matrix.  On  the 
other  hand,  chopped  fiber  metal  matrix  composites  do  not  exhibit 
the  degree  of  layered  anisotropy  observed  in  laminated  continuous 
fiber  polymeric  composites.  Due  to  these  differences,  the 
internal  state  in  metal  matrix  composites  can  be  significantly 
different  from  polymeric  composites.  Accordingly,  a  generalized 
model  was  developed  for  this  material.  Although  the  model  is  an 
extension  of  previous  research  on  polymeric  composites  [4],  the 
mechanics  of  damage  development  are  totally  different.  The 
details  of  this  model  are  given  in  Appendix  6.1.  A  synopsis  of 
the  important  results  is  given  in  this  section. 

The  model  utilizes  the  thermodynamics  with  internal  state 
variables  (ISV)  [5]  to  develop  the  following  stress-strain 
relations : 

^ij  “  cijkl(£kl  “  *1kl  "  *4kl  "  ^kl5  »  (1> 
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stress  tensor. 
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thermal  strain  tensor. 

For  the  case  without  second  phase  material  and  no  damage,  the 
model  reduces  to  classical  plasticity  equations  [6]. 


Furthermore,  for  the  uniaxial  case,  the  model  reduces  to 
Kachanov's  damage  model  [7]  when  plastic  strain  is  negligible. 

Many  of  the  space  structural  applications  of  interest  in 
this  research  involve  uniaxial  stress  states,  such  as  in  the  case 
of  truss  structures,  and  Eu 1 er-Bernou 1  i  i  beam  structures  [8]. 
Therefore,  the  uniaxial  equations  have  been  studied  in  detail. 
These  are,  for  the  isothermal  case: 

°xx  *  ^£xx  “  xx  “  °^4xx  “  £xxJ  * 


where  E  is  Young's  modulus.  A  graphical  r epresentat ion  of  this 
relation  is  shown  in  Fig.  1,  where  it  is  shown  that  the  history 
dependent  damage  parameter  0(4 xx  can  be  determined  by  observing 
the  tangent  modulus  on  unloading. 

The  following  section  details  the  experimental  effort 
developed  under  the  current  contract  to  verify  and  characterize 
the  model,  as  well  to  develop  the  internal  state  variable  growth 
laws,  for  a  typical  chopped  fiber  metal  matrix  composite. 

2.2.2.  Experimental  Research  on  Metal  Matrix  Composites 

The  primary  objective  of  the  experimental  effort  is  to 
develop  a  technique  for  determining  and  evaluating  damage  in 
metal  metrix  composites.  This  technique  must  be  capable  of 
detecting  cracks  and  voids  (free  surfaces)  in  the  composite. 
These  cracks  are  generally  on  the  order  of  microns  in 
characteristic  dimension,  so  that  scanning  electron  microscopy  is 
required  to  measure  the  damage.  Specimens  have  been  loaded  to 
different  levels  and  the  damage  studied  at  each  increment.  Once 
the  amount  of  damage  is  determined  it  can  be  input  into  the 
general  constitutive  model  for  metal  matrix  compo 3 i t e s ( Se e 
Section  2.2.1.). 

The  material  used  in  this  study  was  obtained  from  ARCO 
Metals  Sllag  Operation  in  Greer,  S.C.  The  composition  of  the 
material  is  6061  Aluminum  with  a  twenty  percent  volume  fraction 
of  F-9  silicon  carbide  whiskers.  Plate  is  made  from  the 
materials  by  a  powder  metallurgy  process  and  cast  into  billets. 
The  billets  are  then  rolled,  extruded  or  machined  to  the  desired 
shapes.  The  SiC  whiskers  average  two  microns  in  diameter  and 
twenty  microns  long.  The  composite  has  a  T-6  temper.  Tensile 
test  coupons  have  been  machined  in  accordance  with  ASTM  E^S 
(Tension  testing  of  metallic  materials)  to  the  dimensions  shown 
in  Fig.  2.  For  the  initial  portion  of  the  study  all  specimens 
have  been  machined  with  the  same  orientation  with  respect  to  the 
plate  for  the  purpose  of  uniformity  (with  respect  to  the  SiC 
whisker  or ientat ion).  A  second  phase  of  the  testing  involves  the 
use  of  tensile  test  specimens  oriented  perpendicular  to  the 
Initial  specimens. 

A  mechanical  test  system  was  necessary  to  load  specimens  to 
different  levels  &  measure  responses  so  that  the  damage  at 
various  load  levels  could  be  determined.  An  Instron  model  1125 
screw-driven  test  system  with  2”  wedge  action  grips  was  used. 
Longitudinal  displacement  data  were  obtained  by  the  use  of  an  MTS 
model  632.1  1  B— 20  1"  gage  length  extensometer.  The  displacement 
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information  was  amplified  by  an  in-house  signal  conditioner  and 
load  information  was  amplified  by  the  Instron  controller.  Both 
load  and  displacement  information  were  plotted  by  an  analog  X-y 
recorder  and  input  to  a  computer  data  acquisition  system.  This 
data  acquisition  system  consisted  of  a  DEC  model  PDP  11/23  Plus 
computer  with  an  A  to  D  converter  manufactured  by  ADAC 
Corporation.  With  this  system,  load,  time,  and  longitudinal 
displacement  data  could  be  stored  and  a  real  time  plot  of  load 
versus  displacement  could  be  monitored  during  the  actual  tests. 

The  mechanical  tests  were  performed  by  first  loading  the 
specimen  in  the  Instron,  calibrating  and  zeroing  the  extensometer 
using  an  extensometer  displacement  calibrator,  calibrating  the 
load  cell  and  plotter,  then  securing  the  extensometer  to  the  edge 
of  the  coupon  with  rubber  bands. 

All  tests  were  monotonic  with  a  cross-head  speed  of  .05 
in/min.  Several  specimens  were  tested  to  failure  to  obtain  a 
data  base  about  the  ultimate  strength  of  the  material.  After  the 
mechanical  response  of  the  material  was  determined,  specimens 
were  loaded  past  yield  at  500  pound  increments  so  that  the  damage 
at  each  increment  could  be  determined. 

Once  the  metal  matrix  composite  coupons  were  loaded,  the 
next  step  in  the  process  was  to  prepare  them  for  examination 
with  the  scanning  electron  microscope.  Tensile  coupons  must  be 
sectioned  into  pieces  that  are  about  1/2in  X  1/2in.  This  was 
accomplished  by  using  a  Micro-Matic  precision  slicing  and  dicing 
machine  using  abrasive,  cutting  wheels.  Since  there  was  a 
question  about  the  degree  of  anisotropy  caused  by  the  fiber 
orientation  in  the  plane  of  the  plate,  two  sections  were  cut  from 
each  coupon  so  that  the  two  orthogonal  views  of  the  coupon  could 
be  viewed  under  the  microscope.  The  sectioned  pieces  were  then 
mounted  in  a  conductive  mounting  material  (Konductomet  I) 
manufactured  by  Buehler.  Mounting  was  accomplished  by  the  use  of 
a  Leco  PR-22  Pneumatic  Mounting  Press  with  a  core  cycle  of  nine 
minutes  at  4200  psi  &  300F.  The  Konductomet  I  material  is 
essentially  a  carbon  filled  phenolic  that  is  designed  for 
electron  microscopy  specimen  use.  After  mounting,  the  specimens 
must  be  polished  so  that  the  structure  can  be  seen  in  a  plane. 
Polishing  i3  very  critical  because  if  it  is  not  done  properly, 
detail  can  be  lost  or  polishing  induced  artifacts  will  be 
created.  The  SiC/Al  composite  presents  other  problems.  There  is 
a  vast  difference  between  the  hardness  of  SiC  and  aluminum. 
Unless  great  care  i3  taken  and  an  appropriate  grade  of.  abrasive 
compound  is  used,  an  uneven  terrain  is  created  by  the  removal  of 
the  softer  aluminum  matrix  leaving  exposed  SiC  whiskers.  The 
mounted  specimens  were  wet  sanded  in  one  direction  on  each  of  the 
following  grits  of  sandpaper:  240,  320,  400,  600.  After  each 
grit,  the  specimen  was  washed  in  water  to  remove  any  residue, 
then  sanded  on  the  next  finer  grit  in  a  different  direction  until 
all  remaining  evidence  of  the  previous  sanding  direction  was 
removed.  The  specimen  was  then  cleaned  ul trasonical ly  in  MEK. 
Diamond  abrasive  polishing  compound  was  chosen  for  the  fine 
polishing.  Both  the  polishing  compounds  and  cloths  were  obtained 
from  Buehler.  This  was  a  necessity  because  of  the  hardness  of 
the  SiC  particles.  Rough  polishing  was  accomplished  by  using  a 
polishing  wheel  with  Metadi  seven  micron  polish  on  a  nylon  cloth. 
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Each  specimen  was  polished  until  all  visible  scratches  from  the 
600  grit  sandpaper  were  removed.  The  specimens  were  then 
u 1 tr asoni ca 1 1 y  cleaned  in  MEK  and  hand  polished  with  Metadi  II 
three  micron,  then  Metadi  one  micron  heavy  polishing  compounds  on 
TEXMET  polishing  cloths.  After  each  change  in  polishing 
compound,  it  was  imperative  to  thoroughly  clean  all  residue  of 
the  previous  compound.  If  this  is  not  done  properly,  a  smooth 
polished  surface  is  impossible  to  obtain  because  the  residual 
large  pieces  of  polishing  compound  remove  large  pieces  of 
material  along  with  smaller  pieces  removed  with  the  finer  grain 
compound. 

After  the  polishing  technique  was  performed,  the  polished 
surface  of  the  composite  was  examined  under  an  SEM  to  verify  a 
uniform,  smooth  surface  was  present. 

Once  the  surface  was  polished,  chemical  etching  was 
performed  so  that  surface  details  could  be  easily  seen.  A 
solution  of  10  gr  NaOH  and  90cc  distilled  H20  at  1 50F  was  used 
for  approximately  15  seconds  to  etch  the  specimen.  The  specimen 
was  then  rinsed  and  ul trasonical  1  y  cleaned  in  distilled  water 
followed  by  ultrasonic  cleaning  in  MEK.  The  final  step  in 
preparation  of  the  sample  was  to  apply  a  thin  (“100A)  coating  of 
gold  by  vacuum  deposition.  This  is  necessary  to  achieve  a  good 
image  on  the  SEM.  The  SEM  operates  by  having  a  voltage  applied 
to  the  specimen.  When  the  voltage  is  applied  to  the  specimen, 
free  electrons  are  released  from  the  specimen.  It  is  these 
electrons  that  are  removed  from  the  orbitals  that  are  received  by 
a  sensor  that  forms  the  image.  Light  compounds  have  less  free 
electrons,  thus  the  image  formed  is  not  as  sharp  as  an  image  of  a 
heavy  compound.  Since  gold  is  a  heavy  compound,  a  thin  coating 
of  it  on  the  specimens  provides  an  electron  source,  without 
losing  surface  detail. 

After  the  specimens  were  sectioned,  mounted  in  the 
conductive  mount,  polished,  etched,  and  coated,  the  next  step  in 
the  process  was  to  examine  them  under  the  scanning  electron 
microscope.  A  Joel  JSM-25  II  scanning  electron  microscope  was 
used  in  this  study. 

Choosing  an  acceleration  voltage  for  the  specimens  was  a  key 
consideration.  If  the  voltage  was  too  low  not  enough  electrons 
were  free  to  form  an  image.  If  too  high,  charging  at  interfaces 
and  edges  occurred  leaving  bright  spots  on  the  image  and 
destroying  detail.  After  experimentation  a  value  of  25kv  was 
used.  Most  pictures  were  taken  at  2,000X.  This  magnification 
allowed  details  of  the  matrix,  whiskers,  and  interfaces  to  be 
seen  along  with  more  macroscopic  details  such  as  large  voids. 
The  final  product  was  a  series  of  SEM  photographs  of  the  surfaces 
of  orthogonal  edges  for  specimens  under  both  virgin,  loaded,  and 
failed  condition  (see  Fig.  3).  Along  with  SEM  photographs  of 
polished  surfaces  to  reveal  internal  structure,  photographs  of 
fracture  surfaces  were  taken  to  observe  the  effect  of  the  whisker 
reinforcement . 

Although  this  research  is  not  complete  at  this  time,  the 
photographs  are  currently  under  study  to  determine  the  damage 
state  at  each  load  level.  These  data  will  then  be  fed  into  the 
damage  model  developed  in  Section  2.2.1  in  order  to  evaluate  the 
material  constants  and  complete  the  model  description  for  the 
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Fig.  3.  SEM  Photographs  of  Chopped-Fiber 
Metal  Matrix  Composite  . 


chopped  fiber  metal  matrix  composite. 

The  above  section  details  the  M.S.  thesis  research  of  Mr.  E. 
W.  Nottorf. 

2 . 3  Space  Structural  Response  Due  to  St  if f ness  Loss 

Fibrous  composites  are  known  to  undergo  a  small  but 
significant  amount  of  stiffness  loss  due  to  load  induced 
microcracking  [3].  This  stiffness  loss  usually  occurs  over 
several  hundred  thousand  load  cycles.  Due  to  the  stress 
dependent  nature  of  the  damage,  the  stiffness  loss  is  spatially 
variable  and  concentrated  in  the  areas  of  high  stresses.  This 
spatial  change  in  the  material  properties  of  the  structure 
results  in  appreciable  changes  in  the  dynamic  response  of  the 
structure . 

A  part  of  the  current  research  has  been  to  develop 
approx  imate  -m-ethods  for  determining  this  long  term  change  in 
structural  response  (See  Appendix  6.2.).  The  procedure  developed 
here  is  to  subject  the  structure  to  a  dynamic  load  input  which  is 
in  phase  with  one  of  the  first  few  fundamental  modes  of  the 
structure  in  the  undamaged  state.  Utilizing  the  stress  field 
calculated  from  this  analysis,  it  is  possible  to  estimate  the 
spatial  dependence  in  the  stiffness  loss  that  a  space  structure 
would  undergo  by  approximating  the  constitutive  equations  which 
are  designed  to  model  long  term  material  degradation.  The 
approximation  enters  through  the  negligence  of  the  actual  long 
term  history  dependent  nature  of  the  damage  process.  This 
approximation  is  made  for  reasons  of  numerical  economy.  It  is 
theoretically  possible  to  obtain  the  actual  history  dependent 
response,  but  this  would  require  many  hours  of  computational 
time  on  a  large  mainframe  computer  for  realistic  space 
structures . 

Because  this  procedure  is  approximate,  the  results  should  be 
considered  as  qualitative  in  nature.  Consider,  however,  the  case 
of  a  cantilevered  space  boom,  as  shown  in  Fig.  4.  In  Fig.  5  the 
displacement  of  the  axial  dimension  of  the  boom  is  shown  for  the 
third  mode  as  a  function  of  damage  induced  stiffness  loss  in  the 
second  mode. 

It  is  clear  from  Fig.  5  that  damage  induced  by  cyclic 
loading  in  fibrous  composites  will  alter  considerably  the  dynamic 
response  of  space  structures  even  for  small  amounts  of  damage 
induced  stiffness  loss.  Mode  shapes  are  found  to  change 
dramatically,  with  node  points  shifting  substantial  distances  in 
the  structure.  These  results  indicate  that  active  control 
devices  which  are  placed  based  on  the  initial  undamaged  mode 
shapes  may  require  substantial  relocation  after  the  structure  has 
been  in  operation  for  a  few  months  or  years. 


2.3.1.  Analysis  of  History  Dependence  of  Structural  Response  of 
Simple  Space  Structures  with  Load  Induced  Stiffness  Loss 

Although  it  is  not  possible  to  construct  a  more  precise 
history  dependent  structural  algorithm  for  a  representative  space 
structure,  a  more  accuralte  method  is  under  development  at  this 
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time  for  a  single  beam  member  with  various  boundary  conditions. 
This  model  will  carry  out  the  actual  time  integration  for  the 
slowly  degrading  structure  as  damage  accumulates  for  each  cycle. 
Using  this  new  model  for  simplified  structures,  it  will  be 
possible  to  determine  the  actual  structural  response  for  a  load 
Input  of  several  hundred  thousand  cycles.  The  following  is  a 
brief  description  of  this  procedure. 

The  well  known  partial  differential  equation  for  the  free 
vibration  of  a  beam  is 


where  E  is  Young's  modulus,  I  is  the  moment  of  inertia  of  the 
cross-section,  A  is  the  cross-sectional  area,  p  is  the  mass 
density,  y  is  the  transverse  displacement,  x  'is  the  axial 
coordinate,  and  t  is  time. 

A  number  of  solutions  to  the  above  differential  equation  are 
available  in  the  literature  for  both  uniform  (constant  cross- 
section)  and  nonuniform  (variable  cross-section)  with  different 
boundary  conditions.  Most  of  the  solutions  are  for  beams  with 
homogeneous  material  properties.  These  solutions  have  been 
obtained  by  assuming  that  the  stiffness  of  a  structural  element 
is  constant  in  time  and  therefore  independent  of  loading  history. 
Neither  material  damage  nor  en vironmental  ly  caused  degradation 
are  considered  in  these  analyses. 

Due  to  the  occurrence  of  load  induced  and  history  dependent 
damage  in  composite  materials,  these  previously  obtained  results 
represent  unrealistic  approximations  of  the  actual  structural 
behavior  (See,  for  example.  Appendix  6.2).  In  particular,  the 
resonant  frequencies  and  mode  shapes  of  the  structure  can  be 
severely  altered  by  the  introduction  of  spacially  variable 
damage.  These  parameters  in  turn  can  have  a  substantial  impact 
on  the  active  control  algorithm  to  be  employed  for  control  of 
flexible  body  modes.  By  introducing  material  damage  and 
environmentally  caused  degradation,  the  stiffness  of  a  structure 
is  no  longer  a  constant,  since  it  will  change  substantially 
according  to  the  stress  distribution  and  the  history  of  external 
loading.  The  stiffness  loss  may  change  the  natural  frequencies 
and  mode  shapes  substantially.  With  the  material  damage  and 
environmentally  caused  degradation  involved,  the  differential 
equation  becomes  difficult  if  not  impossible  to  solve  in  closed 
form. 

The  concept  of  internal  state  variables  (ISV)  is  introduced 
to  represent  the  history  dependent  change  of  stiffness.  An 
internal  state  variable  D  is  utilized  as  a  local  ISV  representing 
the  damage  state.  Together  with  the  ISV  growth  law,  the  finite 
element  solution  technique  can  be  modified  to  account  for  the 
history  dependent  stiffness  of  the  beam  element,  with  resulting 
field  equations 

LmII m  \  4- L*0lw]  =  Lo\  ,  (u) 


[/] 


where 


n/vrj  =  J  MiMjfAC*-)**- 

0 

kjj  =•  j  E.(£/Tj-o)T  C*)aY^  ,Vi"  d-Y. 


,  (5) 


The  above  set  of  second  order  ordinary  differential  equations  for 
each  element  is  combined  to  represent  the  eigenvalue  problem  for 
the  beam  structure. 

The  occurrence  of  damage  will  cause  the  loss  of  stiffness, 
that  is,  the  stiffness  is  history  dependent.  Experimental 
results  indicate  that  the  time  scale  for  damage  and  degradation 
is  very  long  compared  to  the  first  fundamental  frequency  of  the 
structure.  Therefore,  the  mathematical  algorithm  is  treated  as 
linear  with  slowly  varying  coefficients.  In  this  research, 
particular  interest  is  being  placed  on  the  natural  vibration 
solution  of  a  beam  structure  with  history  dependent  stiffness  and 
the  investigation  of  the  possible  effect  of  material  damage  and 
stiffness  reduction  on  the  natural  frequencies  and  mode  shapes  of 
planar  beam  structures  with  various  boundary  conditions  (free-? 
free,  clampedffree,  clamped^clamped  and  simply  supported). 

The  research  also  focuses  on  the  investigation  of  the 
internal  state  variable  representation  of  the  damage  phenomenon. 
The  damage  in  a  composite  material  includes  a  sequence  of 
microstructural  and  macrostructural  events  such  as  microvoid 
growth,  matrix  cracking,  fiber  matrix  debonding,  interior  laminar 
cracking,  edge  delamination  and  fiber  fracture.  The  most 
significant  effect  of  damage  on  the  material  properties  is  that 
the  stiffness  will  be  substantially  changed  during  the  life  of 
the  component.  The  constitutive  equation  for  a  composite 
material  could  be  represented  as 

<r-  E ( £ ^  £T)  ,  (6) 

where  E  is  Young's  modulus,  which  will  change  according  to  the 
damage  D  as 

E  -  E0(1  “  D)  .  (7) 

The  subscript  o  represents  the  initial  condition.  The  damage 
D  is  an  internal  state  variable  describing  the  damage  phenomenon 
during  the  life  of  the  composite  structure,  which  is  governed  by 
the  internal  state  variable  growth  law 

D  -  f ( £ , T , D )  .  (8) 

Since  the  damage  phenomenon  is  not  fully  understood  at  this 
time  some  approximations  to  equation  (8)  will  be  made  in  the 
research  in  order  to  reasonably  predict  the  damage  behavior. 

The  above  section  details  the  Ph.  D.  thesis  work  of  Mr.  Y. 
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T.  Chang. 

2.3.2.  Development  of  a  Structural  Model  with  Temperature  Field 
Induced  Damage 

For  many  space  structures  the  primary  source  of  cyclic 
loading  will  be  due  to  thermal  strains  induced  by  solar  and  earth 
radiation.  The  determination  of  the  temperature  field  in  this 
analysis  comprises  a  difficult  matter  in  itself.  The  following 
is  a  description  of  a  model  currently  under  development  which 
will  account  for  thermal  effects  on  the  damage  process. 

The  anticipated  construction  of  large  space  structures  using 
composite  materials  has  stimulated  interest  in  the  relationship 
between  material  damage  and  structural  response.  Early  research 
concludes  that  during  the  normal  life  cycle  of  an  LSS,  damage,  in 
the  form  of  a  reduction  of  stiffness  up  to  25  %,  should  be 
expected.  In  addition,  as  little  as  5%  reduction  in  stiffness 
has  been  shown  to  significantly  alter  structural  response.  The 
current  research  attempts  to  study  in  more  detail  the 
relationship  between  damage  and  stress-strain  distributions  in  a 
simplified  space  structure.  Emphasis  is  placed  on  analyzing  the 
structure  under  a  set  of  accurate  t her mome chan i ca 1 1 y  induced 
loads.  The  full  impact  of  thermal  loads  has  previously  been 
largely  ignored.  Many  studies  have  analyzed  structures  composed 
of  truss  type  elements.  This  results  in  neglecting  bending  by 
the  members  themselves.  For  structures  made  with  materials  of 
high  conductivity,  this  is  not  inaccurate.  However,  composite 
materials  have  low  conductivities  and  therefore  undergo  large 
thermal  gradients  through  their  cross  section.  To  incorporate 
these  expected  thermal  moments,  the  models  to  be  studied  will  be 
composed  of  beam  type  elements. 

The  outlined  boundary  value  problem  is  complicated  by 
several  factors.  First,  a  one  way  coupling  between  temperature 
and  displacements  exists.  It  is  one  way  in  that  displacements 
depend  on  temperatures.  Secondly,  the  problem  is  nonlinear  due 
to  the  introduction  of  radiation  boundary  conditions.  Thirdly, 
there  are  constantly  changing  thermal  loading  conditions  due  to 
varying  earth-structure-sun  orientation.  Finally,  geometrical 
factors  such  as  shadowing  and  interelement  radiation  and 
conduction  exist.  These  factors  combine  to  create  a  highly 
complex  problem  (See  Fig.  6.). 

Several  simplifying  assumptions  to  the  geometry  and 
environment  of  the  structure  can  be  made  without  compromising  the 
usefulness  of  the  result.  For  a  structure  composed  of  long  thin 
members  which  are  sparsely  placed,  interelement  contributions  may 
be  neglected.  Each  member  can  be  treated  as  an  isolated, 
independent  body  absorbing  thermal  radiation  and  in  turn  emitting 
its  own  radiation.  For  a  structure  in  a  geosychr onous  or  other 
high  earth  orbit  the  effects  of  earth  emitted  radiation  and  earth 
reflected  solar  radiation  are  minimal.  Furthermore,  LSS  such  as 
solar  collectors,  antennas,  and  telescopes  will  have  space-fixed 
orientations.  This  will  result  in  steady-state  thermal 
conditions  due  to  solar  radiation  of  nearly  constant  direction 
and  intensity.  This  steady-state  thermal  condition  also  results 
in  negligible  axial  temperature  gradients  (See  Fig.  7.). 
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Fig.  7.  Model  of  Thermal  Problem  for 
Space  Structural  Component. 


The  solution  method  to  be  outlined  here  is  a  completely 
numerical  one.  The  algorithm  consists  of  five  parts,  as  shown 
in  Fig.  8.  Finite  elements  are  used  to  construct  the  temperature 
field  through  the  c r o s s - s e c t i o n  of  a  beam  member.  This 
temperature  distribution  is  then  converted  to  thermal  forces  and 
moments  using  the  following  equations: 


Pr=  J  feVATdA 
My  =  f  E  =<r  AT  Z  <M 
Mf=(E°<TATyJI\  '  <9> 

Here  E  represents  Young's  modulu3,c<T  is  the  coefficient  of 
thermal  expansion,  A  is  cross-sectional  area,  £T  is  the  change 
in  temperature,  and  y  and  z  are  the  distances  to  the  geometric 
centroid.  These  loads  are  first  calculated  in  the  local 
coordinates  of  the  cross-section  and  then  transformed  into  the 
global  coordinates  of  the  structure.  These  global  thermal  forces 
and  moments  then  serve  as  input  to  another  finite  element 
routine.  The  output  from  this  is  the  stress  and  strain 
distribution  in  a  beam  member.  The  final  step  is  to  predict  the 
distribution  of  damage  in  each  member  based  on  the  stress 
history.  Once  the  damage  is  known,  a  new  stress-strain  field  may 
be  generated  for  the  current  damage  state. 

The  beam  code  is  a  standard  transient  one  dimensional  code 
and  is  the  result  of  application  of  the  standard  finite  element 
formulation  to  the  governing  differential  equations: 
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where  uQ,  vQ,  wQ  -  centroidal  displacements  in  the  x,  y,  and  z 

coordinate  directions, 

P,  M„,  M_  -  axial  force  and  moments, 

/  * 

PT,  MyT,  M2t  -  thermally  induced  axial  force  and  moments 
E1  •  Young's  modulus, 

ft  ft 

Iyy  ,  IZ2  -  modulus  weighted  moments  of  inertia,  and 
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GENERATE  TEMPERATURES 
THROUGH  A  CROSS-SECTION 


modulus  weighted  cross-sectional  area. 


For  more  details  about  the  above  formulation  see  reference  8. 

In  order  to  place  equations  (10)  in  a  variational  form,  each 
must  be  multiplied  by  a  separate  test  function  and  integrated 
over  the  volume  of  a  typical  element.  The  resulting  equations 
may  then  be  discretized  in  the  standard  way. 

The  heat  transfer  code  is  less  general.  The  code  is 
developed  from  the  general  formulation  for  heat  transfer  in  a 
plane  geometry: 

ocvT  =  bx(k*$)  +  '&(ky)i£)  +Q 

y  '  ,(11) 

with  boundary  conditions 

ky\y^y  ~  ^ +  h(Ta.'T)  +i4(Trt-T^)%  (12) 

with  variables  defined  as  follows: 


thermal  conductivities, 
internal  heat  generation, 
direction  cosines, 
flux  normal  to  surface, 
absorptivity, 
film  coefficient, 
ambient  temperature, 
emisslvity, 

Boltzman's  constant, 


Tr  -  reference  temperature, 

p  -  density,  and 

Cy  -  specific  heat. 

Because  of  the  steady^state  conditions  the  governing  equation 

reduces  to  the  following: 


O  =  ix/'kx&c)  f  %/  (ky  }yj  +-  Q 


.  (13) 


All  nonlinearity  resides  in  the  boundary  conditions  (12).  To  set 
up  the  finite  element  formulation,  equation  (11)  is  cast  into  its 
variational  or  weak  form  by  multiplying  by  the  test  function 


and  integrating  over  an  element  volume/le: 

^Ti 


O  -  j^lh  (I*  ( kv  s y)J  ST 41 2  i-J^QSrcm 


(14) 


Integration  by  parts  and  the  substitution  of  boundary  conditions 
yields  the  following  equilibrium  equation  in  which  B  is  the 
boundary  of  an  element: 


f  r  /  XT  ^£T 

JnJ-txZx  - 


S  V  JdSl  = J^GiSTdn 


+f[ftot  +  hfa-T)  ^iTrWtroS 


(15) 


Discretization  of  equation  (15)  results  in  the  following  set  of 
equations  for  a  single  element: 


CKLe]t[KMeL])(T^ 


- 


(16) 


KNLI1  contains  the  unknown  temperatures  raised  to  the  3rd  power. 
Because  of  this  the  solution  of  the  equations  requires  that  K  ^ , 
be  approximated  by  a  set  of  known  temperatures  and  then  updated 
through  iteration  until  the  correct  temperature  field  is 
obtained . 

Equation  (16)  may  be  recast  in  the  following  form: 


ifcb]-  (XIM  +  [KmIIM  -  joi  -fo? 

Newton  Iteration  can  then  be  used  to  solve  (17): 

L  =  -jAr. Ji 

I t.-  ^  it  Tk )  +  a  *  kuj 


( 17) 


(18) 


where 


and  n-number  of  degress  of  freedom  per  element  and  AT  is  the 
temperature  correction  on  the  Kth  Iteration.  To  increase  the 
accuracy  of  the  solution,  load  steps  are  used.  Convergence  of 
the  solution  is  assumed  if  ["bf/ST]  remains  non-singular  and  the 
load  steps  are  small.  Equation  (18)  gives  the  final  form  of  the 
equations  to  be  solved: 


[g7khr?kt,=  !Q}k-([KL]  +  CK»J)kM 


(19) 


The  process  of  ultimate  failure  in  a  composite  is  preceded 
by  a  series  of  events  such  as  transverse  cracking,  de 1  ami nat i on , 
fiber  breakage,  and  fiber-matrix  debonding.  This  sequence  of 


microstructural  and  macrostructural  phenomena  is  termed  damage. 
Global  material  properties  such  as  stiffness  and  ultimate 
strength  can  be  substantially  altered  by  an  accumulation  of 
damage . 

This  research  models  damage  as  a  load  history  dependent 
reduction  in  stiffness  of  the  structural  members.  The 
distribution  of  the  damage  will  be  dependent  on  the  resulting 
stress  field.  The  members  most  heavily  stressed  will  incur  the 
greatest  degree  of  damage.  The  intensity  of  the  damage  will  be 
determined  through  the  use  of  a  power  law  degradation  of 
stiffness.  Therefore  the  degradation  of  the  modulus  is  given  by 

e'  =  E0  Ei  -  ,  (20) 

where  E'  is  the  reduced  modulus,  E0  is  the  initial  modulus,  (Tis 
the  stress  in  the  member,  cTmax  is  the  ultimate  strength,  and  n 
is  some  power. 

The  power  law  model  is  simplistic.  In  reality  the  models 
are  expected  to  be  quite  complex.  Work  is  currently  being  done 
to  construct  these  models  for  future  work  in  this  area. 

This  research  has  two  goals.  The  first  is  to  investigate 
the  sensitivity  of  stress-strain  distributions  in  beams  to  damage 
induced  by  t hermome chan i ca 1 1 y  induced  loads  typical  of  their 
proposed  environment.  The  second  is  to  provide  a  numerical 
algorithm  to  be  used  in  future  research  to  provide  thermal 
loadings  for  structural  members.  Expansion  of  the  algorithm  to 
include  interelement  contributions  and  transient  thermal  loads 
is  suggested  to  enlarge  the  geometries  and  environments  that  can 
be  studied. 

The  above  section  details  the  M.S.  thesis  research  Oi  Mr.  J. 
D.  Lutz. 

2 . 4  Space  Structural  Response  Due  to  Damping  Change 

In  a  vibrating  structure  mechanical  energy  is  continuously 
converted  to  other  forms  of  energy  via  irreversible  thermodynamic 
processes,  and  this  energy  loss  is  partially  exhibited  in 
structural  damping.  The  energy  dissipation  may  be  caused  by 
thermal  flux,  material  inelasticity,  friction,  creation  of  new 
boundaries  via  fracture,  chemical  processes,  or  other  sources. 
The  emphasis  of  this  research  is  to  develop  a  model  capable  of 
predicting  this  damping  for  composite  materials  undergoing  load 
induced  damage. 

A  damping  measure  which  is  commonly  used  is  the  quantity  of 
energy  D  dissipated  during  one  cycle  of  harmonic  motion.  The 
maximum  potential  energy  V  stored  in  the  structure  is  related  to 
the  energy  dissipated  by  the  loss  factor  given  by 

~D 

^  “  -2W  y  •  (21) 

When  the  structure  is  built  up  from  elements  with  known 
damping  characteristics  a  proper  damping  matrix  can  be  determined 


yielding  the  fundamental  structural  dynamics  problem 


LmIIx!  +  LclUl  +  Ct]U]- _  (22) 

in  which  x  is  the  vector  of  discrete  coordinates,  f(t)  is  the 
vector  of  forcing  functions,  and  CM],  [C],  and  [K]  are  the  system 
mass,  damping,  and  stiffness  matrices,  assembled  from  the 
individual  matrices  of  structural  finite  elements.  The  damping 
and  stiffness  depend  on  the  damage  in  the  structure. 

For  the  LSS  these  matrices  are  very  large,  and  the 
coordinates  become  highly  coupled.  Thus,  the  direct  solution  of 
the  equations  of  motion  becomes  time-consuming  and  expensive.  A 
very  common  way  to  analyze  damped  structures  is  to  use  a  modal 
approach,  which  is  permitted  when  only  the  response  of  the 
structure  within  a  certain  frequency  range  is  of  interest.  Then 
the  response  can  be  expressed  by  means  of  a  limited  number  of 
vibration  modes.  The  reduction  of  the  extent  of  the  calculations 
is  the  advantage  of  this  approach. 

The  modal  matrix  of  eigenvectors  can  be  calculated  from  the 
undamped  free  vibration  problem 

L/vOJxlf  CfcHxl-  to]  .  (23) 


Modal  equations  can  be  obtained  from  the  transformation 

■  <24> 

where  [$]  is  the  matrix  of  the  eigenvectors  and  q  is  the  vector 
of  time  dependent  generalized  coordinates.  In  general,  the 
damping  matrix  in  transformed  coordinates  is  nondiagonal,  so  that 
the  modal  equations  are  coupled. 

As  an  alternative  to  the  assumption  of  modal  damping,  the 
equations  of  motion  can  be  uncoupled  without  restrictions  on  the 
damping  matrix  [C],  other  than  that  it  be  symmetric.  By 
selecting  a  different  set  of  generalized  coordinates,  equation 
(22)  can  be  reduced  to 

+  YU}  ,  (25) 

where 

q(t)  -  [ xT( t ) : xT ( t ) ]T  is  a  2N  dimensional  state  vector,  and 

Y(t)  -  [fT(t):0T]T  is  a  2N  dimensional  generalized  force 
vector . 


The  associated  eigenvalue  problem  becomes 


-  O 


(26) 


which  can  be  reduced  to  the  form 


A 


1 


(27) 


* 

where,  assuming  that  M  is  nonsingular 


A 


-i  i  -i 

-M  C|-M  K 

- T  -  - 

.  X  °  _ 


The  solution  consists  #of  2N  gigenvalues  X,  and  2N  eigenvectors. 
Because  the  matrices  M  and  K  are  real,  if  is  an  eigenvalue, 
then  Tj_  is  also  |n  eigenvalue.  The  eigenvectors  are  orthogonal 
with  respect  to  M  and  K  . 

The  damping  factor  for  composites  has  been  observed  to 
increase  with  damage.  Plunkett  [11]  found  that  the  damping 
factor  could  be  calculated  from  the  crack  density  as  a  function 
of  strain  level  and  the  strain  volume  fraction. 

If  it  is  assumed  that  the  energy  dissipation  in  a  material 
depends  on  the  local  strain  state  and  its  time  history  the 
damping  factor  can  be  of  the  form 


-  f  (D)  ,  (28) 

where  D  is  the  history  dependent  damage  parameter.  An  internal 
state  variable  theory,  [4]  can  then  be  utilized  to  relate  the 
damping  factor  with  a  given  damage  state.  Thus,  response 
calculations  can  be  made  for  large  space  structures. 

The  research  detailed  in  the  above  section  is  the  Ph.  D. 
thesis  work  of  Mr.  S.  Kal yanasundar um . 


2 . 5  Temperature  Change  Due  to  Hysteretlc  Heating 

It  is  envisioned  that  many  space  structural  components  will 
be  fabricated  from  aluminum  alloys.  Furthermore,  it  has  been 
suggested  that  utilizing  these  components  in  their  post-yielded 
state  can  lead  to  enhanced  passive  damping  of  the  structure  due 
to  hysteretlc  strain  energy  conversion  to  heat.  For  viscoplastic 
structures  this  is  a  complex  nonlinear  two-way  coupled  problem  in 
that  the  conservation  of  energy  is  coupled  via  the  temperature  to 
the  other  field  equations  [9]. 


The  purpose  of  this  portion  of  the  research  has  been  to 
construct  a  model  capable  of  predicting  the  structural  response 
of  a  viscoplastic  space  structure  to  thermomechanical  loading. 
The  primary  emphasis  is  placed  on  the  prediction  of  the 
temperature  rise  accrued  in  typical  viscoplastic  members  due  to 
hysteretic  loss  during  cyclic  loading. 

The  general  formulation  of  this  problem  was  completed  during 
the  first  year  of  research  [10].  During  the  second  year  the 
model  has  been  applied  to  a  typical  space  truss  structure  made  of 
aluminum  members  (See  Appendix  6.4).  The  structure  is  assumed  to 
be  coated  with  special  high  emissivity  materials,  and  is 
subjected  to  various  levels  of  solar  and  earth  radiation.  The 
problem  is  complicated  not  only  by  the  thermomechanical  coupling, 
but  also  by  the  nonlinear  constitutive  equations  and  radiative 
boundary  conditions. 

As  shown  in  Fig.  9,  the  predicted  temperature  rise  for  an 
input  cyclic  loading  at  5  Hz  can  be  very  substantial  in  spite  of 
the  fact  that  the  radiation  boundary  conditions  provide  a 
structural  heat  loss  mechanism. 

The  algorithm  developed  in  reference  10,  although 
compr ehens i ve  in  nature,  is  very  computationally  expensive  to 
utilize.  In  fact,  it  is  so  costly  as  to  preclude  the  analysis  of 
the  response  of  a  structural  member  undergoing  several  thousand 
load  cycles.  Therefore,  alternative  procedures  were  considered 
for  analyzing  the  temperature  profile  as  a  function  of  time. 
This  research  has  resulted  in  the  completion  of  the  three  part 
paper  described  in  Appendices  6.5  through  6.7. 

The  principal  result  obtained  from  this  research  is  the 
determination  of  upper  and  lower  bounds  for  the  hysteretical ly 
induced  temperature  rise  as  a  function  of  load  history,  geometry, 
and  material  properties.  As  shown  in  Fig.  10,  these  bounds 
verify  the  results  obtained  in  Appendix  6.4. 

On  the  basis  of  these  results  it  is  concluded  that  material 
inelasticity  in  viscoplastic  structural  members  3hould  be  used  as 
a  passive  damping  mechanism  only  with  extreme  caution.  The 
temperature  rises  associated  with  this  damping  process  can  be  so 
severe  as  to  compromise  the  integrity  of  the  structure. 


TEMPERATURE  CHANGE 
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ABSTRACT 


A  continuum  mechanics  framework  is  utilized  herein 
to  construct  constitutive  equations  for  metal  matrix  composites 
with  damage.  Matrix  plasticity  and  microcracking  are  modelled  via 
the  concept  of  internal  state  variables.  Imposition  of 
thermodynamic  constraints  results  in  a  set  of  stress-strain 
relations  which  are  dependent  on  damage  and  plastic  strain.  These 
equations  are  specialized  to  a  one-dimensional  case  and  it  is 
then  demonstrated  how  the  material  parameters  may  be  determined 
experimentally. 


INTRODUCTION 


The  characterization  and  modelling  of  matrix 
plasticity,  cracks  and  other  forms  of  microscopic  damage  in  metal 
matrix  composites  is  developed  herein  using  the  principles  of 
continuum  mechanics.  The  theory  developed  incorporates  globally 
averaged,  history  dependent,  thermodynamically  constrained 
constitutive  relations  and  utilizes  tensor  valued  internal  state 
variables  to  model  history  dependent  energy  dissipative 
(irreversible)  phenomena. 

The  two  primary  energy  dissipative  phenomena 
considered  are  inelastic  strain  and  damage.  Depending  on  the 
complexity  of  the  material  under  consideration,  these  internal 
state  variables  can  be  scalar  or  tensor  valued  functions.  Since 
these  phenomena  are  inherently  history  dependent,  the 
specification  of  the  internal  state  variables  representing 
inelastic  strain  and  damage  and  their  associated  growth  laws 
introduces  history  dependence  into  the  boundary  value  problem. 

There  exist  many  phenomena  in  the  microstructure  of 
a  material  which  can  be  classified  into  energy 
dissipative/irreversible  processes.  These  include  crack 
formation,  dislocation  movement  and  arrangement,  grain  boundary 


sliding,  chemical  changes  and  frictional  losses  due  to  rubbing  of 
fractured  surfaces.  Therefore,  some  degree  of  clarification  is 
required  regarding  how  inelastic  strain  and  damage  are  to  be 
defined.  In  this  paper,  all  damage  will  refer  to  cracks  in  the 
material  microstructure. 

Results  obtained  from  previous  research  into  the 
constitution  of  crystalline  solids  [1,2]  indicate  that  for  the 
class  of  materials  into  which  metal  matrix  composites  is 
classified,  the  inelastic  strain  tensor  is  treated  as  a  second 
order  tensor  valued  internal  state  variable.  The  primary 
mechanisms  of  inelastic  strain  are  postulated  as  those  of 
dislocation  density  (drag  stress)  and  dislocation  arrangement 
(back  stress)  which  are  themselves  zero  order  (scalar)  and  second 
order  tensor  valued  internal  state  variables  respectively. 
Research  into  the  kinematics  of  crack  initiation  and  growth 
presented  in  this  paper  as  well  as  reference  [3]  postulates 
damage  to  be  a  second  order  tensor  valued  internal  state 
variable. 

The  damage  is  assumed  to  be  statistically 
homogeneous  within  a  representative  volume  element,  which  is 
assumed  to  be  small  in  comparison  to  the  body  of  interest.  Under 
the  condition  of  small  scale  statistical  homogeneity,  all 
continuum  based  conservation  laws  are  assumed  to  be  valid  on  a 
global  scale  in  the  sense  that  all  changes  in  the  continuum 


problem  resulting  from  internal  damage  are  reflected  only  through 
alterations  in  the  constitutive  behavior.  Therefore, 
microstructural  phenomena  such  as  cracks  qualify  as  damage  and 
their  effects  on  the  performance  of  a  material  can  be  reflected 
through  the  constitutive  relations.  However,  macroscopic  and/or 
nonhomogeneous  damage  states  such  as  large  scale  surface  cracks 
are  treated  as  boundary  effects  which  must  be  reflected  in 
conservation  laws  via  changes  in  the  external  boundary  conditions 
rather  than  in  the  constitutive  relations. 


DEVELOPMENT  OF  THE  CONSTITUTIVE  RELATIONS 


We  now  proceed  to  develop  specific  constraints  on 
constitutive  behavior.  In  this  analysis  the  following  assumptions 
are  made: 


a.  thermomechanical  coupling  is  non-zero, 


electrical  and  magnetic  effects  can  be 
neglected  [4]  as  these  can  be  either 
controlled  or  or  their  effect  on  the  body  of 
interest  calculated, 


c.  infinitesimal  deformations,  and 


d.  in  the  absence  of  damage  or  at  constant  damage 
state  the  material  behavior  is  initially 
linear  thermoelastic  and  isotropic. 


Field  Parameters 


Consider  a  body  occupying  a  closed  region  V  with 
material  points  Xj^  -  In  an  analagous  method  to  that 
proposed  by  Coleman  and  Gurtin  [5],  it  is  postulated  that  the 
following  state  variables  are  required  to  fully  characterize  the 
state  of  all  material  points  Xj  within  a  body  at  all  times  t: 


a. 


b. 


c. 


d. 


e. 


displacement  field 


stress  tensor 


body  force  per  unit  mass 


heat  flux  vector 


internal  energy  per  unit 


ui  “  UjCXj ,t)  (1 ) 

°ij  "  °ij(xk’t)  (2) 

ft  -  f^Xj.t)  (3) 

qi  “  qi(xj,t)  (4) 

mass  u  -  u(Xj,t)  (5) 


f. 


heat  supply  per  unit  mass 


r  -  r(Xj,t)  (6) 


g 


entropy  per  unit  mass 


(7) 


s  -  S(Xj.t) 

h.  absolute  temperature  T  <■  T(Xj,t)  (8) 

and 

i.  »ijk  ■  k  - ' . n  <9) 

where  a^j  are  a  set  of  n  second  order  tensor  valued  internal 
state  variables  necessary  to  model  inelastic  deformation  and 
damage . 


Field  Equations 


The  following  pointwise  field  equations  ,  written 
in  differential  form,  are  assumed  to  hold  for  all  media 
undergoing  infinitesimal  strains: 

a.  conservation  of  linear  momentum 

“u.j  *  >ri  •  »at  (,0) 

b.  conservation  of  angular  momentum  (assuming 
negligible  body  moments): 


(11) 


°ij  “  °ji 


c. 


kinematics: 


eij  "  1/2  (ui,j  +  uj,i) 


(12) 


d. 


conservation  of  mass: 


dp/dt  ■  0 


(13) 


e.  conservation  of  energy: 


eu  •  '  qj,j  *  <,r 


(14) 


f.  second  law  of  thermodynamics 


pY  =  ps  -  pr/T  +  (q./T)  . 

J  9  J 


(15) 


where  Y  is  the  specific  entropy  production  rate.  Equation  (15) 
together  with  the  Clausius-Duhem  inequality: 


Y  2  0 


(16) 


is  assumed  to  hold  for  all  processes.  It  should  be  noted  that 
although  equation  (15)  introduces  one  additional  equation,  it 


simultaneously  introduces  a  new  variable,  namely  Y.  Therefore, 
equation  (15)  alleviates  no  degrees  of  freedom.  However,  although 
inequality  (16)  cannot  itself  specify  any  degrees  of  freedom,  it 
will  impose  constraints  on  the  allowable  form  of  the  constitutive 
equations. 


In  the  foregoing  analysis  the  body  force  vector, 
inertial  effects  and  body  moments  will  be  assumed  to  be 
negligible  and  the  conservation  of  mass  is  trivially  satisfied. 


Helmholtz  Free  Energy 


mass: 


We  now  define  the  Helmholtz  free  energy  per  unit 


h  =  u  -  Ts  (17) 

where  h  is  the  arithmetic  sum  of  the  total  internal  energy  (u)  of 
a  body  and  the  internal  entropy  production  (irrecoverable  energy 
loss)  (Ts)  dissipated  during  a  process.  Thus,  the  Helmholtz  free 
energy  represents  the  net  available  internal  energy  within  a  body 
following  any  process. 


It  follows  that: 


u  ■  h  +  Ts  +  Ts 


(18) 


Equation  (18)  can  be  substituted  into  equation  (14)  and  this 

*  result  then  combined  with  equation  (15)  and  inequality  (16)  to 

give: 

pTy  -  -ph  -  pTs  *  °klekl  -  ^jgj/T  -  0  (19) 

where 

8j  ■  t,j  <20) 

The  field  equations  (10),  (11),  (14),  (15)  and  (16) 
are  combined  with  conditions  applied  on  the  boundary  B  of  a  body 
to  specify  the  complete  field  problem. 

I 

Constitutive  Relations 

* 

Suppose  that  the  displacements  u  and  temperature  T 
are  specified.  Then  the  method  of  Coleman  and  Noll  [6]  can  be 
used  to  obtain  the  spatial  and  time  distribution  of  the  body 
force  vector  (f)  and  the  heat  supply  (r)  from  the  conservation  of 

C  linear  momentum  and  energy  equations  respectively,  assuming  that 

the  displacement  and  temperature  fields  are  specified. 


Constitutive  relations  can  then  be  constructed  [7]  from  the 


remaining  state  variables  and  their  spatial  derivatives.  The  form 
of  the  constitutive  relations  used  including  internal  state 
variables  and  their  associated  growth  is  postulated  as  follows: 


°ij  ’  °ij(ekl’T’gk'aklP)  (21) 

u  «  uUkl.T,gk,aklP)  (22) 

s  -  s(ekl.T,gk,aklP)  (23) 

qi  ■  qi(ekl*T’gk’“klP)  (24) 

and 

V  ■  aki,<*krT-v«kip>  (25> 


The  form  of  equations  (21)  through  (24)  implies 
that  all  constitutive  equations  are  evaluated  in  the  specified 
3tate  (Xj,t).  For  this  reason  j ,  u,  s  and  q^^  are  termed 
observable  state  variables  since  they  can  be  determined  from 
equations  of  state  for  all  times  t  even  though  there  is  implicit 
history  dependence  via  the  internal  state  variables  aklP.  In 
addition,  from  the  definition  (17)  and  the  form  of  equations  (22) 
and  (23)  it  follows  that: 


Note  that: 


1)  the  constitutive  relations  (21)  through  (24)  are 
equations  of  state  and  therefore  cannot  contain  time  rates  of 
change  of  state  variables; 

2)  the  principles  of  local  action  and  equipresence  [8,9] 
hold  for  equations  (21)  through  (24),  but  need  not  be  satisfied 
in  the  growth  laws  [equation  (25)]; 

3)  the  superscripts  p  and  q  in  equation  (25)  range  from 
one  to  the  number  of  internal  state  variables  required  to  fully 
characterize  the  inelastic  response  and  the  damage  state  of  the 
body; 

4)  the  growth  laws  by  their  definition  embody  time 
dependence  and  can  therefore  contain  both  spatial  and  time 
derivatives  of  state  variable  arguments.  If  equations  (25)  are  at 
all  times  integrable  in  time  then  the  following  is  an  equivalent 
expression  to  (25): 

aij1<^xm,t^  */  dt'  (26a) 

**—00 

where  t  is  the  time  of  interest  and  t’  is  a  dummy  variable  of 

[/ 

integration.  Therefore,  it  is  apparent  that  are  not  directly 
observable  at  any  time  and  must  be  considered  as  internal  or 
hidden  state  variables. 


Internal  State  Variables  Considered 

When  a  body  denoted  by  B  is  subjected  to  some 
traction  or  deformation  history  as  shown  in  Figure  1 ,  it  will 
undergo  a  thermodynamic  process  which  will  in  general  be 
irreversible  to  some  degree. 

This  irreversibility  is  introduced  by  such 
phenomena  as  dislocation  movement  and  growth,  fracture  (both 
micros  and  macroscale),  friction  (due  to  rubbing  of  fractured 
surfaces),  grain  boundary  sliding,  deformation  twinning  and 
chemical  changes.  However,  due  to  the  relative  magnitudes  of 
these  energy  dissipative  processes  the  two  primary  internal  state 
variables  postulated  to  be  of  interest  in  this  analysis  are  those 
of  inelastic  strain  (a^j)  representing  dislocation  phenomena  and 
damage  (a^)  representing  microcracking.  Although  there  may 
exist  some  interrelationship  between  these  variables,  they  are 
considered  to  be  independent  as  they  represent  different  physical 
phenomena. 

The  two  microstructural  mechanisms  which  are 
considered  to  have  the  greatest  prominence  [10]  at  temperatures 
less  than  one  half  the  melting  temperature  of  the  material  under 
consideration  are: 

a.  Drag  Stress  (a2):  representing  1ocally  averaged 


dislocation  density  and  producing  isotropic 
hardening;  and 

b.  Back  Stress  (a^j):  representing  the  locally 
averaged  dislocation  arrangement  and  producing 
kinematic  hardening  (Bauschinger  effect). 


Factors  Influencing  the  Growth  of  Inelastic  Strain. 


The  primary  independent  internal  state  variables 
influencing  the  growth  of  damage  are  considered  to  be  the 
historical  inelastic  strain  and  damage  states,  back  stress  and 
drag  stress.  The  growth  law  for  the  inelastic  strain  is 
therefore  postulated  as  follows  [11]: 


X1  i j  "  ^ekl,T,a1kl’0t2,o3kl’0t4kl^°ij  ' 


(27) 


where : 


a2  “  °2(ekl,T,0t1kl’Ql2,a3kl,a4kl) 


(28) 


a3ij  “  a3ij<ekl,T,ailJ,a2,a3ij»04kl)  (29) 


X  -  scalar  valued  function  of  state 


and 


deviatoric  stress  tensor 


It  should  be  noted  that  the  back  stress  and  drag 
stress  represent  mechanisms  which  physically  exist.  However, 
dislocation  arrangement  and  density  are  microstructural 
phenomena.  Although  they  influence  the  growth  of  the  inelastic 
strain,  they  do  not  directly  enter  into  the  stress-strain 
relations  because  they  are  not  kinematic  quantities. 

Factors  Influencing  the  Growth  of  Damage. 


Damage  is  defined  in  this  analysis  as  intergranular 
mechanisms  such  as  grain  boundary  sliding,  microvoid  growth  and 
microstructural  cracks. 

As  the  arrangement,  density  and  growth  of 
dislocations  will  effect  the  growth  parameters  (rate,  direction) 
of  any  damage  existing  in  the  microstructure,  the  primary 
internal  state  variables  influencing  the  growth  of  damage  are 
considered  to  be  inelastic  strain,  the  previous  damage  state, 
back  stress  and  drag  stess.  Therefore,  the  growth  law  for  damage 
is  postulated  as: 


a4ij  “  0Uij(ekl,T,o1kl'a2’a3kl’°4kl) 


(30) 


r 


for  isotropic  conditions: 

a4ij  "  f*ekl,T,a1kl,a2,a3kl'a4kl)6iJ  (31  J 

It  should  be  noted  that  although  damage  is 
obviously  a  directionally  related  quantity  and  therefore 
tensorial  in  nature,  it  is  difficult  to  distinguish 
phenomenologically  between  damage  and  drag  3tress  since  both  can 
interpreted  as  stiffness  reducing  mechanisms. 

Functional  Form  of  the  Helmholtz  Free  Energy. 

As  the  inelastic  strain  and  damage  are  postulated 
to  be  the  primary  independent  internal  state  variables  affecting 
the  constitution  of  the  materials  considered  herein  and  both  are 
functions  of  the  back  stress  and  drag  stress,  the  Helmholtz  free 
energy  is  postulated  to  have  the  following  form: 

h  "  h  (ekl,T,8k,01kl,a4kl)  (32) 


Thermodynamic  Constraints. 


Thermodynamic  constraints  on  the  form  of  the 
constitutive  relations  [equations  (21)  through  (24)]  can  be 
accomplished  using  the  Coleman-Mizel  procedure  [12]  under  the 
assumption  of  the  Claus ius-Duhem  inequality  [equation  (16)].  In 
this  procedure,  the  Helmholtz  free  energy  [equations  (32)]  is 
implicitly  differentiated  with  respect  to  its  arguments  as 
follows: 


h  -  dh/de.  ♦  dh/dT  T  +  dh/dg  g 

kx  kl  '  m  m 

*  dh/Soit  ij“t  ij  * 


(33) 


Substitution  of  equation  (33)  into  the  inequality  (19)  results 
in: 


pTY  *  I-0kl(emn,T,gm,amn  ^  ~  p^h/^ekl(emn,T,gm,amn  )]ekl 

-  (34) 
‘  8ljk('WT«.-0 

’  VWr>«-<Sml)*J/T  *  ° 


where  i  assumes  the  values  of  1  and  4  representing  the  inelastic 
strain  and  damage  respectively.  On  the  basis  of  this  technique  it 
can  be  shown  that  the  satisfaction  of  the  First  and  Second  Laws 
of  Thermodynamics  will  lead  to  the  following  conclusions  [2]: 

h  -  h(ekl,T,a1kl.a4|<1)  (35) 

Ojj  -  pdh/de^  (36) 

s  *  -  dh/dT  (37) 

u  -  h  +  Ts  (38) 

and  aiJq  -  flijq(ekl»T*°1kl»02,a3kl'a4kl)  (39) 

Note  that: 


1  )  thermodynamically  constraining  the 
constitutive  relations  removes  the  dependence  of  the  Helmholtz 
free  energy  on  the  temperature  gradient  gk. 

2)  the  above  thermodynamically  constrained  forms 
of  the  constitutive  relations  [equations  (36)  through  (38)]  show 
an  intrinsic  dependence  of  stress,  entropy  and  internal  energy  on 
the  Helmholtz  free  energy. 


Following  the  Coleman-Mizel  procedure,  inequality 


(34)  reduces  to: 


pTY  -  -pdh/do^  niJk  -  Qjgj/T  *  0  (40) 

where  the  first  term  is  called  the  internal  dissipation  and  the 
last  term  is  dissipation  due  to  heat  conduction. 

Coleman  and  Gurtin  [5]  utilize  inequality  (40)  to 

prove  that: 


■  Vj  *  °<  lgj'  >  (,M) 

at  an  asymptotically  stable  equilibrium  state  where  is  the 
positive  semidefinite  thermal  conductivity  tensor.  It  will  be 
assumed  herein  that  these  conditions  hold  and  the  higher  order 
terms  in  equations  (41)  will  be  neglected. 

Therefore,  the  specification  of  the  Helmholtz  free 
energy  function  [equations  (35)]  will  completely  define  the 
constitutive  equations  (21)  through  (23)  and  equations  (41)  will 
define  the  constitutive  equations  (24).  Thus,  if  the  internal 
state  variable  growth  laws  (39)  can  be  determined  the  field 
problem  will  be  completely  specified. 


The  Local  Averaging  Process 

Constitutive  relations  (21)  through  (24)  are 
theoretically  pointwise  in  nature;  that  is,  they  are  applicable 
to  fixed  infinitesimal  material  points.  However,  there  is  no 
practical  method  to  construct  experiments  on  material  points 
since  at  the  microscopic  level  the  continuum  assumption  becomes 
invalid.  Therefore,  to  obtain  constitutive  relations  which  are 
applicable  to  a  continuum,  it  is  considered  acceptable  to 
construct  constitutive  equations  by  subjecting  local  specimens  to 
surface  deformations  or  tractions  which  lead  to  spatially 
homogeneous  stress  and  strain  fields  so  that  some  local  average 
of  the  pointwise  observable  state  variables  can  be  determined 
directly  from  the  effects  on  the  boundaries  of  the  specimen. 

As  shown  in  Figure  2,  the  scale  of  the  smallest 
dimension  of  a  local  specimen  is  generally  constructed  to  be  at 
least  an  order  of  magnitude  larger  than  the  scale  of  the  largest 
material  inhomogeneity.  This  sizing  assists  in  preserving  the 
continuum  assumption  while  still  averaging  out  the  effects  of 
point  defects  such  as  crystal  lattice  dislocations.  Conversely, 
the  scale  of  the  largest  dimension  of  a  typical  specimen  should 
be  small  compared  to  the  scale  of  the  boundary  value  problem  of 
interest.  This  constraint  is  necessary  in  order  to  preserve  the 
notion  that  the  constitutive  equations  are  indeed  pointwise  in 


nature 


Dimension  of  Local  Specimen 


The  local  rather  than  pointwi.se  constitutive 
equations  that  result  from  experimentation  are  assumed  to  be  of 
the  same  form  as  pointwise  equations  (21)  through  (25).  For 
example,  in  the  uniaxial  test  described  in  Figure  2  it  is 
customary  to  define  the  following: 

o^dXgdx^  (42) 


and 

T  ST(a1(a2,a3)  (44) 

where  L  is  the  local  specimen  gage  length,  A  is  the 
cross-sectional  area  in  the  Xg-x^  plane,  and  (a^a^a^)  is  some 
arbitrary  point  on  the  surface  of  the  specimen.  Utilizing  these 
quantities,  it  is  then  hypothesized  that: 


f  £11 o*1 


where 


°1 1  ( e1 1  ,T,°klP)  "  °11  (e11  ,T,otklP) 


S,.,P  -  1/V  /  a,  P 


®kl  dxidx2dx3 


(45) 

(46) 


and  all  quantities  with  bars  represent  the  locally  measured  state 
variables. 


Although  equations  (45)  represent  an  often  used  way 


V* 


of  relating  pointwise  to  experimental  results,  the  local 


averaging  process  is  nevertheless  prone  to  shortcomings  since  the 
definitions  (42)  through  (44)  all  represent  nonunique  relations 

k 

between  pointwise  state  variables  and  o  9  e..f  T,  a,,  and  their 

ij  ij 

_  _  _  u 

locally  defined  counterparts  5  ;,  t..,  T,  a, .  .  There  are  in  fact 

ij  ij  ij 

an  infinite  number  of  distributions  aklP(x1 .Xg.x^)  which  will 
result  in  identical  values  for  However,  assuming  that  the 

scale  of  any  inhomegeneities  is  small  and  that  the  distribution 
of  oij^  is  random,  the  specimen  will  be  statistically  homogeneous 
and  the  relation  between  oj<1P  and  aklP  will  be  approximately  one 
to  one.  Note  that  for  clarity  the  overbar  will  be  excluded  from 
expressions. 


Determination  of  the  Constitutive  Equation  Form. 


To  evaluate  the  constitutive  form  of  the  stress 
tensor  and  the  entropy  from  the  equalities  given  in  equations 
(36)  and  (37),  the  pointwise  Helmholtz  free  energy  per  unit  mass 
is  expressed  as  a  second  order  Taylor  Series  expansion  in  the 
independent  state  variable  arguments  as  follows: 
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+  1/2Jijkla4ija4kl  +  Kij£ijAT  +  Lija1ijAT 


+  Mija4ijAT  +  Nijkleija1kl  +  °ijkleij°4kl 


+  Pijkla1ij°4kl 


If  second  order  effects  are  neglected  then  the 
stress  tensor  can  be  evaluated  from  equation  (36)  by 
differentiating  equations  (47)  with  respect  to  the  strain  tensor 
as  follows: 


°u  -  <’ah/Jeij 
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Now  define: 


pBij  "  °ij 


”P^ij AT  "  Cijklekl 


rpNijkla1kl  "  Cijkla1  kl 


(48a) 


"p0ijkl°4kl  "  Cijkl°4kl 


Substitution  of  equations  (48a)  into  equations  (48)  yields: 


R  T 

°ij  *  °ij  +  Cijkl(ekl  "  “ikl  "  a4kl  +  ekl  } 


(49) 


where 


a^.  -  residual  stress  tensor,  and 


e.  .  -  thermal  strain  tensor, 

ij 


The  expression  for  entropy  can  be  readily  obtained 
by  differentiating  equation  (47)  with  respect  to  temperature. 


A  coupled  heat  equation  can  also  be  obtained  by 
substituting  equation  ( 1 8 )  into  (14),  imposing  the  conditions  of 
equations  (35)  through  (37)  and  then  subtituting  the  Helmholtz 
free  energy  expression  into  the  result  [3]. 


Description  of  the  Internal  State. 


Consider  an  arbitary  region  denoted  by  B  as  shown 
in  Figure  3-  (The  body  B  is  assumed  to  be  small  on  a  scale  to  a 
boundary  value  problem  of  Interest).  Now  consider  an  arbitary 
local  element  denoted  as  L  with  external  surfaces  SI  which  are 
chosen  normal  to  a  set  of  Cartesian  coordinate  axes  (x^.x^x^)  as 
shown  in  Figure  4.  Although  the  element  L  is  an  intrinsic  part  of 
the  body  B,  it  can  be  thought  of  as  being  removed  from  B  and  the 


newly  created  surfaces  subjected  to  appropriate  boundary 
conditions  so  that  the  response  of  the  element  L  to  an  imposed 
boundary  condition  is  identical  to  that  of  the  element  when  it  is 
located  in  B.  The  scale  of  L  is  chosen  so  that  its  dimensions  are 
small  compared  to  the  dimensions  of  B.  However,  the  dimensions  of 
L  are  also  large  enough  to  guarantee  statistical  homogeneity  of 
the  material  properties  and  existing  defects  even  though  the 
total  surface  area  of  the  defects  may  be  the  same  order  of 
magnitude  as  SI . 

Internal  surfaces  resulting  from  fracture  are 

labelled  S2.  The  volume  of  the  element  is  defined  to  be  V,  which 

includes  the  volume  of  any  damage  (cracks  and/or  voids)  denoted 

as  V  . 
c 

In  order  to  describe  the  internal  state,  we  first 
consider  the  kinematics  of  a  typical  point  0  with  neighboring 
points  A  and  B  as  shown  in  Figure  5.  Before  deformation  lines  OA 
and  OB  are  orthogonal  as  shown  in  Figure  5a.  After  deformation  we 
imagine  that  the  lines  Joining  O', A*  and  B'  are  as  shown  in 
Figure  5b  and  just  at  the  instant  that  deformation  is  completed, 
a  crack  forms  normal  to  the  plane  of  AOB  through  point  O'  as 
shown  in  Figure  5c.  Furthermore,  point  O’  becomes  two  material 
points  O’  and  0"  on  opposite  crack  faces  and  points  A'  and  B' 
deform  further  to  points  A'1  and  B''.  It  is  assumed  that  all 
displacements,  including  crack  opening  are  infinitesimal  so  that 


a)  Point  O  Prior  to  Deformation, 


b)  Point  O  After  Deformation  and  Prior 
to  Fracture  Process, 

c)  Point  O  After  Fracture. 


%  <*  '  •  •  «  *  » 

1  .  . 


the  observer  at  an  appropriate  observation  distance  from  point  0 
sees  only  the  deformation  A'  '  O'  B".  The  strain  associated  with 
this  deformation  is  appropriately  called  an  observable  state 
variable.  However,  the  strain  of  interest  is  associated  with 
A"0,,B".  Therefore,  it  is  essential  to  construct  an  internal 
state  variable  which  will  relate  these  two  strain  descriptions. 

We  therefore  construct  the  vectors  u  connecting  O' 
and  0*'  and  nc  describing  the  normal  to  the  crack  face  at  O'  as 
shown  in  Figure  6.  It  should  be  noted  that  uG  can  be  used  to 
construct  a  pseudo^strain  representing  the  difference  in  rotation 
and  extension  of  lines  A"0"B". 


The  rate  of  change  of  surface  energy  released  per 
unit  local  volume  due  to  cracking  in  L  [3]  is  given  by: 


°ijuinj 


(50) 


Assuming  that  any  cracks  in  the  body  recover  elastically,  that 
is,  close  upon  the  removal  of  the  load,  then  we  further  assume 
that  this  process  is  reversible  and  that  tractions  f  can  be 
applied  at  point  O'  which  will  close  the  crack: 


u>l\ 


(51) 


Using  Cauchy's  formula  Equation  (46)  can  be  placed  in  the  form: 


crack  face 
tangent 


crack  faces 


.-v  ."V 


UL  “ 
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dS 


(52) 


where  the  subscripts  denote  quantities  associated  with  the  actual 
crack  geometry. 


Guided  by  the  fact  that  uG  and  nc  describe  the 
kinematics  of  the  cracking  process  at  point  0,  we  now  define  the 
following  second  order  tensor  valued  internal  3tate  variable  to 
describe  damage: 


c  c 

“ij  “  ui  nj 


(53) 


We  now  define  the  locally  averaged  internal  state 
variable  describing  damage  to  be: 


“HIJ  ■  k/» 


LVL  /  Vj 
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dS 


(54) 


where  Sc  -  surface  area  of  cracks,  and 


k  -  a  constant  resulting  from  the  simplification  of 
equations  (48)  to  equations  (49). 


Mathematical  Model  for  Damage 


The  damage  model  considered  in  this  paper 
classifies  damage  as  a  displacement  related  mechanism;  current 
models  [13]  classify  damage  as  a  stiffness  reducing  mechanism. 
For  the  case  of  inf initesimal  strains,  the  following  section 
demonstrates  that  the  model  proposed  herein  is  consistent  with 
these  existing  models. 

Kachanov's  reduced  stiffness  model  Cl 3]  is  given 
by: 


°u  ‘  Bu  *  Cijki’(cki  -  “iki '  £ki  >  (55) 


Equations  (55)  can  be  shown  to  be  equivalent  to  equations  (48a) 
by  first  equating  as  follows: 


Now  let: 


Cijkl(ckl  "  °1kl  "  °4kl  '  ekl  } 


Cijkl'(ekl  '  °1 kl  ’  ekl  5 


Cijkl'  ’  Cijkl  "  ^ijklmn^mn 


(55a) 


(56) 


then  (55a)  becomes: 
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UNIAXIAL  FORM  OF  THE  MODEL 


For  the  purpose  of  this  analysis,  the  temperature 
field  is  assumed  to  be  significantly  less  than  1/2Tm.  Therefore, 
if  thermal  strain  is  neglected,  the  uniaxial  isothermal  form  of 
equation  (44)  is  given  by: 

o  -  oR  +  E(e  -  a.j  -  a^)  (58) 

Consider  the  input  strain  profile  shown  in  Figure 
7.  The  stress^strain  curve  resulting  from  a  uniaxial  tension  test 
of  a  metal  matrix  coupon  using  this  strain  profile  is  given  in 
Figure  8.  De.io^e  the  slope  of  the  elastic  portion  of  the  load-up 
curve  as  (E)  and  the  slope  of  the  unloading  curve  as  (Et). 

The  slope  (Efc )  is  defined  as: 


Fi  gure  7 


r 


Efc  -  do/de  (59) 

Thus,  differentiating  Equation  (58)  with  respect  to  the  strain 
tensor  gives: 

-  E(  1  -  da.j/de  -  da^/de)  (60) 

If  the  specimen  is  unloaded  following  the 
application  of  a  traction  or  deformation  field  the  inelastic 
strain  remains  constant  as  the  strain  energy  (oe)  has 
substantially  decreased  due  to  the  decrease  in  the  applied  stress 
(o).  This  is  shown  diagramatically  in  in  Figure  9. 

Therefore,  for  unloading: 

da.j/de  -  0  (61) 

and  Equation  (60)  becomes: 

Et  -  E(  1  -  da4/de)  (62) 

Now,  if  the  aluminium  matrix  is  linear  elastic  on 

unloading  then  Et  is  constant  on  unloading  and  the  implication  is 
that: 

da^/de  -  constant  (on  unloading) 


(63) 


The  curve  corresponding  to  this  result  is  given  at  Figure  10. 


EXPERIMENTAL  VERIFICATION  OF  MODEL 


In  order  to  qualitatively  verify  the  supposition 
that  the  inelastic  strain  tensor  (a^j)  can  be  regarded  to  be  an 
internal  state  variable,  consider  the  example  of  a  uniaxial  bar 
subjected  to  applied  displacements  such  that  the  end  tractions 
will  be  evenly  distributed.  It  is  customary  to  deduce  the 
inelastic  strain  in  an  experiment  of  this  t  ype  by  utilizing  the 
output  from  a  loadcell  to  determine  the  stress  and  then  make  use 
of  equations  (58)  to  determine  the  elastic  strain.  However, 
equations  (58)  entail  the  damage  tensor  (a^j).  Therefore  an 
experimental  procedure  is  required  to  distinguish  between 
deformation  due  to  inelastic  strain  (back  stress  and  drag  stress) 
and  that  due  to  damage.  The  procedure  proposed  is  demonstrated  by 
equations  (61)  and  (63)  and  illustrated  in  Figures  9  and  10. 

Therefore,  it  is  proposed  that  the  inelastic  strain 
and  damage  resulting  from  an  applied  stress  or  deformation  field 
can  be  distinguished  experimentally;  the  inelastic  strain  is 
obtained  from  the  unrecovered  strain  and  the  damage  is  related  to 
the  change  in  slope  between  the  loadup  and  unloading  portions  of 


the  stress-strain  curve. 


CONCLUSIONS 


The  foregoing  analysis  has  developed  a  constitutive 
model  for  the  inelastic  deformation  of  a  continuum  with  damage. 
Following  the  application  of  thermodynamic  constraints  on  the 
form  of  the  constitutive  relations  and  a  kinematically 
justifiable  definition  of  the  Helmholtz  free  energy,  a 
constitutive  equation  relating  the  stress  tensor  to  the  total 
strain  tensor,  the  inelastic  strain  tensor  and  the  damage  tensor 
was  developed.  Experimental  results  from  uniaxial  tests  carried 
out  on  aluminium  matrix  material  with  silicon  carbide  fibres  has 
determined  that  a  global  measure  of  microstructural  damage  could 
be  directly  related  to  the  slope  differential  between  the  load-up 
and  unloading  curves  on  the  uniaxial  stress-strain  plot  [14]. 
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EFFECT  OF  DEGRADATION  OF  MATERIAL  PROPERTIES 
ON  THE  DYNAMIC  RESPONSE  OF  LARGE  SPACE  STRUCTURES 


S.  Kalyanasundaram*,  J.  D.  Lutz*,  W.  E.  Haisler**,  and  D.  H.  Allen*** 
Texas  A AM  University 
College  Station,  Texas  77843 


Abstract 

In  this  paper  the  effect  of  degradation  of 
material  properties  on  structural  frequencies  and 
mode  shapes  of  Large  Space  Structures  (LSS)  Is 
Investigated.  The  difficulty  and  cost  of 
maintenance  of  LSS  make  It  a  necessity  to  design 
these  structures  to  operate  with  a  certain  amount 
of  load-induced  damage.  This  damage  Is  commonly 
observed  in  fibrous  composite  media. 

Sensitivity  studies  conducted  on 
representative  space  truss  structures  indicate  that 
degradation  of  material  properties  may  have  a 
significant  effect  on  the  structural  mode  shapes 
and  frequencies.  For  even  small  amounts  of 
reduction  In  stiffness  (10*),  frequencies  and  nodal 
locations  may  change  significantly.  It  Is  clear 
that  these  effects  must  be  taken  Into  consideration 
when  designing  control  systems  for  Large  Space 
Structures. 


Introduction 

Due  to  economic  constraints,  it  Is  projected 
that  advanced  high  strength-to-weight  ratio 
aerospace  materials  will  be  utilized  in  future 
generation  space  structures.  Such  materials 
Include  polymer  and  metal  matrix  fibrous 
composites,  which  are  known  tp  jjndergo  a  certain 
amount  of  load  Induced  damage,  ’  These  materials 
are  also  expected  to  undergo  a  certain  amount  of 
environmentally  Induced  damage  or  degradation,  thus 
resulting  in  significant  stiffness  losses. 

Experimental  research  on  advanced  composite 
materials  Indicates  that  the  material  may  undergo 
up  to  15  percent  loss  in  stiffness  due  to 
thermomechanical  fatigue,  which  causes  a  variety  of 
damage  modes  In  the  structure.  Additional  loss  of 
stiffness  may  be  attributed  to  elevated  temperature 
and  chemical  changes  due  to  solar  radiation  and 
other  environmental  effects.  This  reduction  In 
stiffness  affects  the  dynamic  response  which  In 
turn  is  critical  In  the  development  of  control 
systems  for  LSS.  In  this  paper,  sensitivity 
studies  will  be  presented  which  Investigate  the 
effect  of  stiffness  loss  on  structural  frequencies 
and  mode  shapes. 

The  advent  of  the  space  shuttle  has  made 
possible  the  development  of  LSS.  Control  systems 
for  stabilizing  and  maneuvering  these  very  large 
space  structures,  especially  those  for  precise 
pointing,  will  require  extension  of  current 
technology. 
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Although  large  size  by  itself  does  not  arouse 
concern,  structural  flexibility  resulting  from 
minimizing  the  structural  weight  In  non- 
gravltatlonal  fields  may  present  problems. 
Extremely  large  structural  flexibility  may  result 
In  large  amplitudes  and  low  frequencies  (.01  to  10 
Hz)  which  may  create  new  complications  for  control 
designers. 

As  an  example  of  the  precision  required1  ,  a 
typical  radiometry  application  may  utilize  a  200 
meter  antenna  with  an  effective  beam  width  of  0.01 
degrees  and  have  requirements  limiting  the 
vibratory  beam  shift  to  less  than  0.005  degrees  and 
dynamic  surface  distortions  to  less  than  1mm. 
Maneuvering  or  maintaining  the  altitude  of  such  a 
satellite  leads  to  flexible  body  motion  which  must 
be  well  predicted  and  controlled. 

The  Importance  of  Interaction  between  control 
systems  and  vibratory  response  has  caused 
considerable  research  In  LSS  control  systems. 

The  current  practice  of  guaranteeing  a  large 
separation  between  modal  frequencies  and  the 
bandwidth  of  control  will  not  be  adequate  in  future 
applications.  The  combination  of  large  size  and 
payload-weight  restrictions  will  drive  structural 
frequencies  down  and  the  need  for  more  accurate 
pointing  will  drive  the  control  system  bandwidth 
up.  When  sufficient  frequency  separation  becomes 
impossible,  there  exists  a  need  for  adaptive 
control  systems.  This  leads  to  further  research  In 
the  design  of  structural  control  systems  actuator/ 
sensor  placement,  and  distributed  sensing  and 
actuation  as  opposed  to  co-located  sensors  and 
actuators . 

Techniques  for  achieving  modal  control  of  LSS 
will  require  a  more  accurate  knowledge  of  modal 
characteristics.  Optimum  sensor  and  actuator 
placement  will  be  greatly  Influenced  by  modal 
effects  which  must  be  known  to  a  greater  degree  of 
precision. 


Problem  Summary 

In  order  to  investigate  the  possible  effects 
of  material  degradation  on  the  dynamic  response  of 
LSS,  a  representative  space  truss  structure  has 
been  selected  In  the  shape  of  a  long  boom  as  shown 
In  Fig.  1.  Using  several  loading  histories,  stress 
distributions  have  been  obtained  for  each  truss 
member.  The  resulting  stress  distributions  can  be 
used  In  a  material  damage  model  to  define  material 
degradation  and  resultant  stiffness  reductions. 
Using  the  reduced  stiffness  properties,  modal 
analyses  have  been  conducted  on  the  structure  to 
show  the  effect  of  material  degradation  on  natural 
frequencies,  mode  shapes  and  nodes.  Details  of  the 
finite  element  model,  material  degradation  model, 
and  numerical  results  are  presented  below. 


Fig.  1  Space  Truss  Structure 


Model  Description 
Material  Degradation  Model 

The  process  of  ultimate  failure  of  composite 
materials  Is  preceded  by  a  sequence  of 

mlcrostructural  and  macrostruci  jral  events  which 
are  termed  as  damage.  These  events  may  be  due  to 
transverse  cracking,  delamlr^atlpn,  fiber  breaking 
and  fiber-matrix  debonding.  The  mechanical 

response  of  the  structure  Is  affected  by  this 
damage.  Global  material  properties  like  stiffness 
and  residual  strength  may  be  substantially  alceced 
during  the  life  of  the  structural  components.  ” 
Some  of  the  analytical  studies  Car  modeling  damage 
include  ^  shear  lag  concept,  fracture  based 
concepts,'  _1)  and  internal  state  variable 

theories.1  Although  Important  progress  has 

been  made,  current  understanding  of  damage  Is  not 
complete. 


Damage  in  polymeric  composites  is  modelled  In 
this  paper  as  a  load  history-dependent  reduction  In 
stiffness  In  each  structural  element.  The  internal 
state  variable  theory  (ISV)  Is  used  for  modeling 
mechanical  behavior  and  _tbq  stress  strain 
relationship  Is  of  the  form,1  i 
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In  this  case,  the  ISV  are  assumed  to  be 
second  order  tensor  valued  and  to  enter  only 
through  the  modulus  tensor.  C’ ,  Is  the 
effective  modulus  tensor  given  by  J 
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variables  which  are  given  by  the  following  set  of 
ISV  growth  laws, 


mn 


Q  mn(ekl,T,akl  5 
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At  low  homologous  temperatures  these 
materials  are  assumed  to  be  rate  insensitive  so 
that  the  above  model  will  result  In  quasi-elastic 
(rate  independent)  equations  in  which  Inelasticity 
Is  reflected  only  through  the  slowly  degrading 
modulus  tensor.  Experimental  evidence0' 

Indicates  that  the  time  scale  for  degradation  of 
C'  ...  is  very  long  compared  to  the  frequencies  and 
mode  shapes  of  representative  structures.  It  Is 
therefore  sufficient  for  many  space  structural 
applications  to  treat  equations  (1)  In  the  degraded 
state  only. 


The  stress-strain  relationship  for  the  truss 
elements  Is  a  one-dimensional  approximation  of 
equations  (2)  given  by 
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where  oxx  and  exx  are  the  uniaxial  stress  and 

T 

strain,  e  Is  the  thermal  strain,  and  E'  Is  the 
axial  stiffness  of  the  truss  element  given  by 


E*  -  E  (1  -  o)  (5) 

where  E  is  the  undegraded  axial  stiffness  and  a  Is 
a  scalar  valued  parameter  representing  the 
Integrated  effect  of  all  damage  modes  such  as 
matrix  cracking,  interlaminar  fracture,  fiber 
breakage,  and  fiber-matrix  debonding. 


Experimental  research  on  composite  materials 
lnd. cates  a  power  law  degradation  of  axial 
stiffness  as  a  function  of  stress  history.1  ’  i 
Hence  the  damage  ISV  growth  law  Is  assumed  to  be  of 
the  form 


o  -  k.  (o/o  )n  (6) 

1  max 

where  k.  and  n  are  material  parameters,  o  Is  the 
maxlmvn1 stress  In  the  structure,  and  o  Is  the  axial 
stress  in  each  truss  element.  For  constant  stress 
amplitude,  equation  (6)  may  be  integrated  In  time 
to  give  the  following  approximation 

c(t,)  -  k',  Coa^/o^]"’  (7) 

where  k'  and  n'  are  material  parameters  which  may 
be  time  dependent. 

A  power  law  form  of  damage  Is  used  herein  for 
simplicity  and  for  an  Initial  attempt  at  modeling 
the  structural  response  with  damage.  In  reality 
the  damage  laws  will  be  more  complex  1  «jnd  are 
currently  being  developed  for  future  work.1 

Finite  Element  Model 

Figure  1  Illustrates  the  geometry  of  the 
representative  space  truss  used  to  simulate  an 
antenna  boom.  This  structure  Is  sixty  feet  long 
with  10  bays,  six  feet  long  by  three  feet  wide. 
The  finite  element  model  has  124  space  truss 
elements  and  44  nodes.  In  the  initial  undegraded 
configuration,  the  material  properties  are  the  same 
for  all  members  with  the  following  values: 
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Material  type:  Graphite  epoxy  (Hexel) 

Young's  modulus  E  -  21,5x10°  gsl 

Cross  sectional  area  <  1.0  in 

Density  -  0.065  lazin'3  _< 

Coefficient  of  thermal  expansion  «  2x10  in/ in/  F 

Reference  temperature  ■  89.6  F 

Each  truss  member  is  idealized  with  a 
standard  six  degree  of  freedom  truss  element  of 
constant  cross  section.  Because  the  structure  is 
idealized  as  linear  with  slowly  varying  material 
properties,  conventional  linear  finite  element 
methodology  may  be  used  tp-write  global  equations 
of  equilibrium  of  the  form1  * 

[M]{q]  ♦  [K]{ql  -  «3)  (8) 

where  [M]  is  the  mass  matrix,  [K]  is  the  stiffness 
matrix,  (qj  is  the  nodal  displacement  vector,  and 
{ Q }  is  the  nodal  force  vector.  The  stiffness 
matrix  [K]  is  dependent  on  the  spatially  variable 
damage  state  a  which  varies  from  element  to 
element.  Standard  eigenvalue  extraction  may  be 
performed:  ir.  this  case,  subspace  iteration  was 
used  to  obtain  the  first  five  frequencies  and  mode 
shapes. 

Spatial  Distribution  of  Degradation 

The  spatial  distribution  of  degradation  and 
stiffness  reduction  of  LSS  will  be  complex  and 
dependent  on  loading  and  environmental  history. 
For  the  present  investigation,  wherein  material 
degradation  is  assumed  to  be  a  function  of  stress 
history,  it  was  necessary  to  make  some  assumptions 
about  the  corresponding  stress  history  and  spatial 
distribution  of  stresses  within  the  LSS. 

Two  approaches  were  used  to  obtain  candidate 
stress  histories/distributions  for  predicting  the 
stiffness  degradation.  In  one  case,  the  stress 
distribution  was  obtained  for  an  assumed  thermal 
load  hlstory/dlstributlon.  Secondly,  a  modal 
approach  was  used  wherein  it  was  assumed  that 
primary  degradation  occurred  in  the  first  two 
bending  modes  of  the  structure.  After  computing 
the  mode  shapes  for  the  first  two  undegraded 
bending  modes,  the  nodal  displacements  were  used  to 
compute  a  corresponding  stress  distribution. 

In  each  case,  the  degradation  model  given  by 
equation  (7)  was  then  used  to  obtain  degraded 
properties  for  each  truss  member  assuming  that  the 
element  stressed  the  most  was  degraded  a  specified 
percentage.  The  resultant  structure  with  degraded 
properties  has  spatially  variable  stiffness  that 
varies  from  element  to  element.  Mode  shapes  and 
frequencies  were  then  computed  with  varying  maximum 
percentages  of  degraded  properties. 


Discussion  of  Results 

Natural  frequency  and  mode  shape  responses 
have  been  obtained  for  several  stress-induced 
degradation  test  cases  as  described  above  for  the 
representative  space  tr^s  structure  shown  in  Fig. 
1.  This  particular  truss  structure  geometry, 
representing  a  segment  of  a  boom,  is  similar  to 
ones  being  used  for  other  PAC0SS  related  work. 
Assuming  the  boom  is  fixed  on  one  end  (at  x»0),  the 
five  lowest  frequencies  (for  the  virgin  structure) 
are  equal  to  3.4  Hz,  4.5  Hz,  4.6  Hz,  19.2  Hz,  and 
20.3  Hz.  The  first  mode  is  a  combined  torslon- 
lnplane  shear  mode,  the  next  two  modes  are  bending 


modes  about  the  z  and  y  axes,  respectively,  and  the 
fourth  mode  is  a  pure  torsion  mode. 

The  first  case  considers  the  boom  structure 
shown  in  Fig.  1  (which  is  assumed  to  be  fixed  on 
one  end)  with  a  thermal  gradient  over  the  cross- 
section.  It  is  likely  that  one  surface  of  the 
space  structure  will  become  significantly  hotter 
than  the  other  surface  due  to  solar  heating, 
attitude  of  the  structural  elements  and  shadowing 
effects.  To  investigate  the  effect  of  this  thermal 
gradient  through  the  depth  of  the  truss,  the 
stresses  in  each  element  were  calculated  by 
specifying  a  temperature  of  122°F  for  the  members 
on  the  top  surface,  80.6°F  for  the  members  on  the 
bottom  surface  and  100°F  for  the  diagonal  members 
connecting  the  top  and  bottom  surface.  With  this 
thermally-induced  stress  distribution,  the  axial 
stiffness  of  each  element  was  degraded  by  using 
equation  (7).  The  maximum  level  of  degradation 
(loss  of  stiffness)  was  set  to  a  prescribed 
percentage  for  the  element  with  the  highest  stress 
and  remaining  elements  were  degraded  according  to 
their  stress  level  by  using  equation  (7).  The 
value  of  n’  in  equation  (7)  was  assumed  to  be  0.75. 

In  Fig.  2  the  first  three  natural  frequencies 
are  plotted  for  different  levels  of  damage.  The 
effect  of  damage  on  the  natural  frequencies  is 
clear.  Increasing  the  level  of  damage  reduces  the 
stiffness  of  the  space  truss  and  this  in  turn 
drives  the  natural  frequencies  down  significantly 
even  for  modest  damage  states.  For  a  maximum  loss 
of  25S  in  axial  stiffness  (for  the  highest  stressed 
members),  the  first  three  natural  frequencies  are 
reduced  by  about  8X.  Since  mode  shapes  are 


Fig.  2  Effect  of  Damage  on  Natural  Frequencies 


Fig.  3  Effect  of  Degradation  on  Second  Mode 


Important  for  designing  the  control  systems  of  the 
large  space  structures.  It  Is  desirable  that  they 
be  constant  with  time.  Although  it  was  found  that 
there  was  no  appreciable  change  in  the  first  mode 
shape  between  the  undegraded  and  degraded  cases, 
higher  modes  were  altered  due  to  material 
degradation.  Figure  3  Is  a  plot  of  the  z 
displacement  for  the  second  mode  shape  along  the 
length  of  the  space  truss  (z-0,  y-0).  Significant 
changes  In  the  mode  shape  and  node  locations  as  a 
function  of  percent  degradation  are  observed.  The 
sign  of  the  modal  displacement  is  reversed  near  the 
free  edge  for  the  degraded  and  undegraded  cases  and 
the  location  of  the  node  (zero  displacement) 
changes  appreciably.  Figure  4  is  a  similar  plot  of 
the  y  displacement  along  the  length  of  the  space 
truss  for  the  third  mode. 

The  value  of  n'  In  equation  (7)  was  varied 
from  0.25  to  1.0  to  study  its  effect  on  the  mode 
shapes.  It  was  found  that  the  trend  In  mode  shape 
changes  was  similar  for  different  n'.  Figure  5 
Illustrates  this  point.  Here  the  z  displacement 
for  the  second  mode  is  plotted  along  the  length  of 
the  space  truss  for  different  values  of  n'  (maximum 
reduction  in  axial  stiffness  was  201).  The  plot 
indicates  that  increasing  n*  (l.e.,  decreasing  the 
nonlinearity  of  the  degradation  model)  tends  to 
Increase  the  changes  in  the  modal  displacements. 
Such  nonlinearity  becomes  increasingly  important 
when  stresses  vary  spatially  over  the  structure, 
l.e.,  some  members  are  highly  stressed  compared  to 
others. 


Fig.  4  Effect  of  Degradation  on  Third  Mode 
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Fig.  5  Effect  of  Material  Degradation  Exponent 
on  Second  Mode  Shape 


The  next  two  sample  cases  consider  the 
situation  where  we  assume  that  primary  degradation 
occurs  In  the  first  two  bending  modes.  For 
simplicity,  It  Is  assumed  that  damage  occurring  In 
one  mode  does  not  affect  the  damage  In  any  others, 
l.e.,  no  damage  induced  coupling  of  modes.  In 
reality,  this  may  not  be  the  case  and  will  be 
considered  In  future  research. 

In  the  first  case,  we  consider  the  case  where 
degradation  has  occurred  In  the  first  bending  mode, 
l.e.,  degradation  is  based  on  stresses  calculated 
from  the  modal  displacements  corresponding  to  the 
second  mode  shape.  Figure  6  shows  the  resulting 
first  three  natural  frequencies  for  different 
levels  of  damage.  For  a  maximum  reduction  In 
stiffness  of  25%  the  first  three  natural 
frequencies  decrease  by  8.6%,  9.2%  and  7.6%, 
respectively.  There  Is  little  change  in  the  first 
mode  shape  for  the  degraded  and  undegraded  cases. 
Figure  7  Is  a  plot  of  the  z  displacement  for  the 
second  mode  shape  along  the  length  of  the  space 
truss  and  shows  that  the  modal  displacements  change 
quite  drastically  for  the  degraded  structure.  The 
displacement  at  the  free  edge  Is  nearly  30  times 
the  magnitude  of  the  undegraded  case  for  a  maximum 
damage  of  25%  (the  sign  of  the  displacement  Is  also 
reversed)  and  the  location  of  nodes  also  change 
considerably.  Figure  8  Indicates  similar  changes 
In  the  y  displacement  for  the  third  mode  shape. 
The  fourth  mode  (torsional)  Is  relatively 
unaffected  by  the  degradation  of  material  stiffness 
properties.  This  Is  as  expected  because  the 
present  analysis  assumed  that  primary  degradation 
occurred  In  the  first  two  bending  modes.  Different 
results  would  be  expected  If  significant  stiffness 
reduction  occurred  In  the  primary  torsion  mode. 


Fig.  7  Effect  of  Degradation  on  Second  Mode 
Assuming  Second  Mode  Damage  State 


Fig.  6  Effect  of  Damage  on  Natural  Frequencies 
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Fig.  8  Effect  of  Degradation  on  Third  Mode 
Assuming  Second  Mode  Damage  State 


V  V 


Results  have  also  been  obtained  for  the  case 
where  damage  Is  assumed  to  occur  In  the  third  mode 
(second  bending  mode).  As  in  the  previous  examples 
there  is  no  appreciable  change  in  the  first  mode 
shape  between  the  undegraded  and  degraded  cases. 
The  z  displacement  corresponding  to  the  second  mode 
shape  is  plotted  in  Fig.  9  for  different  levels  of 
damage.  The  displacement  at  the  free  edge  is  very 
large  in  the  damaged  states  as  compared  to  the 
undegraded  state.  Figure  10  illustrates  similar 
results  for  the  third  mode  shape.  These  results 
show  that  the  mode  shapes  and  node  points  may 
change  significantly  for  even  small  damage  amounts. 


Conclusions 

This  study  has  attempted  to  investigate  the 
possible  effects  of  material  damage  and  stiffness 
reduction  on  the  modal  response  of  LSS.  Large 
space  structures  constructed  of  fibrous  composites 
will  experience  some  stiffness  reductions  produced 
by  load-induced  and  environmentally-induced  damage 
of  the  material.  To  what  extent  this  will  occur  is 
uncertain  at  this  point  but  even  small  damage 
amounts  appear  to  be  significant. 

The  present  work  has  shown  that  load-induced 
degradation  of  material  properties  may  have  a 
significant  effect  on  the  structural  frequencies 
and  mode  shapes.  For  the  representative  boom 
structure  considered  here,  even  small  amounts  of 
material  stiffness  degradation  (10S)  produce 
frequency  and  node  shifts  which  appear  to  be 
significant.  It  is  not  Inconceivable  that  mode 
shapes,  node  locations,  and  frequency  distributions 
will  change  over  the  plant  design  life  in  such  a 
way  that  the  structure  response  is  very  much 
different  from  the  virgin  structure.  Such  changes 
in  plant  response  would  require  "robust"  control  of 
a  nature  which  may  not  be  possible  with  present 
technology.  Consequently,  it  is  important  that 
these  effects  be  taken  into  consideration  when 
designing  the  control  systems  for  large  space 
structures. 
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Fig.  9  Effect  of  Degradation  on  Second  Mode 
Assuming  Third  Mode  Damage  State 


Although  preliminary,  this  study  suggests  the 
need  for  a  more  accurate  knowledge  of  the  physical 
nature  of  material  degradation  in  fibrous 
composites,  its  Influence  on  structure  stiffness, 
and  how  material  degradation  will  affect  the  long¬ 
term  modal  characteristics  for  large  space 
structures. 


Acknowledgement 


This  work  was  sponsored  by  the  Air  Force 
Office  of  Scientific  Research  under  Contract  No. 
F49620-83-C-0067. 


144  III  It*  440  441  444 

HORIZONTAL  POSITION 


References 

1.  Herzberg  R.J.,  Johansen  K.F. ,  and  Stroud  R.C., 
"Dynamics  and  Control  of  Large  Satellites," 
Astronautics  and  Aeronautics,  Vol.  16,  Oct.,  1978 
pp.  35-39. 

2.  Skelton  R.e. ,  "Algorithm  Development  for  the 
Control  Design  of  Flexible  Structures,"  NASA-CP* 
2258,  1982. 

3.  Skelton  R.E.,  "Large  Space  System  Control 
Technology  Model  Order  Reduction  Study,"  NASA-CP- 
2118,  1979. 


o>ho  oamass 

0«  IK  0 AMASS 
A-  14M  0 AMASS 
♦•IS*  0 AMASS 
AMASS 
»>f**OAMAS* 


10  Effect  of  Degradation  on  Second  Mode 
Assuming  Third  Mode  Damage  State 


■  . * ,  »  „  • .  * 
.IV.v.v.V.y' 


•  *  •  «\ 


J|  UsAi,  V 


1 


4.  Kroner  C.C. , "Optimum  Damper  Location  for  a 
Free-Free  Beam,"  NASA-CP-21 68. 

5.  Relfsnider  K.L.,  Schultz  K.,  and  Duke  J.C.. 
"Long-Term  Fatigue  Behavior  of  Composite 
Materials,"  Long-Term  Behavior  of  Composites,  ASTM 
STP  813,  1983,  pp.  136-159. 

6.  Relfsnider  K.L. ,  and  Jamison  K.,  "Fracture  of 
Fatigue-Loaded  Composite  Laminates,"  Int.  J. 
Fatigue,  pp.  187-197,  October  1982. 

7.  Masters  J.E.,  and  Relfsnider  K.L. ,  "An 
Investigation  of  Cumulative  Damage  Development  In 
Quasi-Isotropic  Graphite/Epoxy  Laminates,"  Damage 
in  Composite  Materials,  ASTM  STP  775,  1982,  pp.  40- 
62.  , 

8.  0  Brlen  T.K.,  "Characterization  of  Delamlnation 
Onset  and  Growth  in  a  Composite  Laminate,"  Damage 
in  Composite  Materials,  ASTM-STP  775,  1982,  pp. 
140-167. 

9.  Whitney  J.M.,  "Fatigue  Characterization  of 
Composite  Materials,"  Fatigue  of  Fibrous  Composite 
Materials,  ASTM-STP  723,  1981,  pp.  133-151. 

10.  Hlghsmith  A.L. ,  Stinchcomb  W.W.,  and 
Relfsnider  K.L. ,  "Stiffness  Reduction  Resulting 
from  Transverse  Cracking  In  Fiber-Reinforced 
Composite  Laminates,"  VPI-E-81.33.  Virginia 
Polytechnic  Institute,  November  1982. 

11.  Chou  P.C.,  Wang  A.S.D.,  and  Miller  H. , 
"Cumulative  Damage  Model  for  Advanced  Composite 
Materials,"  Material  Laboratory  Air  Force  Wright 
Aeronautical  Laboratories,  Report  No.  AFWAL-TR-82- 
4089. 


12.  Allen  D.H. ,  Groves  S.E.,  and  Schapery  R.A. ,  "A 
Damage  Model  for  Continuous  Fiber  Composites  -  Part 
I:  Theoretical  Development,"  Texas  A&M  Mechanics 
and  Materials  Center,  Report  No.  MM  5023-84-17, 
Feb.  1985. 

13.  Groves,  S.E.,  Allen,  D.H.,  and  Schapery,  R.A., 
"A  Damage  Model  for  Continuous  Fiber  Composites  - 
Part  II:  Model  Applications,  "Texas  A  AM  Mechanics 
and  Materials  Center,  in  preparation. 

14.  Schapery,  R.A.,  "On  Constitutive  Equations  for 
Viscoelastic  Composite  Materials  with  Damage," 
National  Science  Foundation  Damage  Workshop,  May  4- 
7,  1980. 

15.  Mlkulas  M.M.,  Bush  H.G.,  and  Card  M.F. , 

"Structural  Stiffness,  Strength  and  Dynamic 

Characteristics  of  Large  Tetrahedral  Space  Truss 
Structures,"  NASA  TM  X-74001 ,  March  1977. 

16.  Bathe  K.J. ,  and  Wilson  E.L.,  "Numerical 
Methods  In  Finite  Element  Analysis,"  Prentice-Hall 
Inc.,  1976. 

17.  Zlenklewicz,  O.C.,  "The  Finite  Element 
Method,"  Third  Edition,  McGraw-Hill  Book  Company, 
1978. 


7 


Mechanics  and  Materials  Center 
TEXAS  A&M  UNIVERSITY 
College  Station,  Texas 


PREDICTED  AXIAL  TEMPERATURE  GRADIENT 
IN  A  VISCOPLASTIC  UNIAXIAL  BAR  DUE 
TO  THERMOMECHANICAL  COUPLING 


BY 

D,  H.  Allen 


84-15 


NOVEMBER  1984 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

unclassified _ 

2a.  security  classification  authority 

NA 


2b.  QEC LASSI F ICATI ON/OOWNG  R  A  01 NG  SCHEDULE 

NA 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

MM  4875-84-15 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Bb.  OFFICE  SYMBOL 
(if  applicable) 

Aerospace  Engineering  Dept. 

NA 

6c.  ADORESS  (City,  State  and  ZIP  Code) 

College  Station,  TX  77843 

8a.  NAME  OF  FUNOING/SPONSORING 

iSf®^orce°8ff ice  of  Scientifi 
Research 

8b.  OFFICE  SYMBOL 
(If  applicable ) 

Sc.  ADDRESS  tCity.  State  and  ZIP  Code) 

Washington,  D.C.  20332 

11.  title  /include  Security  ciauificattont  Predicted  Axial  Temp, 

lb.  RESTRICTIVE  MARKINGS 

_ NA _ 

3.  OISTRI8UTION/AVAI  LABILITY  OF  REPORT 

Unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

NA 


Air  Force  Office  of  Scientific  Research 


7b.  AODRESS  (City.  State  and  ZIP  Code > 


Washington,  D.C.  20332 


10.  SOURCE  OF  FUNDING  NOS. 


itU  PW  iiWWiw  Rjl 


12.  PERSONAL  AUTHOR(S) 

David  H.  Allen 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO. 

13a.  type  of  REPORT 

Interim 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM _  TO  . 


14.  DATE  OF  REPORT  (Yr..  Mo..  Day ) 

November  1984 


15.  PAGE  COUNT 


FIELD 

GROUP 

COSATI  CODES 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  numben 

The  thermomechanical  response  of  a  uniaxial  bar  with  thermoviscoplastic  constitution 
is  predicted  herein  using  the  finite  element  method.  After  a  brief  review  of  the 
governing  field  equations,  variational  principles  are  constructed  for  the  one  dimensional 
conservation  of  momentum  and  energy  equations.  These  equations  are  coupled  in  that  the 
temperature  field  affects  the  displacements  and  vice  versa. 

Due  to  the  differing  physical  nature  of  the  temperature  and  displacements,  first 
order  and  second  order  elements  are  utilized  for  these  variables,  respectively.  The 
resulting  semi-discretized  equations  are  then  discretized  in  time  using  finite  differenc¬ 
ing.  This  is  accomplished  by  Euler's  method,  which  is  utilized  due  to  the  stiff  nature 
of  the  constitutive  equations. 

The  model  is  utilized  in  conjunction  with  with  stress-strain  relations  developed 
by  Bodner  and  Partom  to  predict  the  axial  temperature  field  in  a  bar  subjected  to  cyclic 
mechanical  end  displacements  and  temperature  bourdary  conditions.  It  is  found  that 


20.  oistribution/availability  of  abstract 
UNCLASSIFIEO/UNLIMITEO  SI  SAME  AS  RPT.  □  OTIC  USERS  □ 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

David  H.  Allen 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22b.  telephone  number 

( Include  Area  Codei 


[_  DO  FORM  1473,  83  APR 


EOITION  OF  1  JAN  73  IS  OBSOLETE. 


22c.  OFFICE  SYMBOL 


security  classification  of  this  page 


S1CUWITV  CLASSIFICATION  OP  THIS  PAGE  _ _ __ — — — 

spacial  and  time  variation  of  the  temperature  field  is  significantly  affected  by  the 
boundary  conditions. 


2/3 


AD-8162  139  A  MODEL  FOR  PREDICTING  THERHOHECHANICAL  RESPONSE  OF 
LARGE  SPACE  STRUCTURE.  .  (U)  TEXAS  A  AND  H  UNIV  COLLEGE 
STATION  MECHANICS  AND  MATERIALS  RE.  .  D  H  ALLEN  ET  AL. 
UNCLASSIFIED  JUN  85  MM-4875-85-11  AFOSR-TR-85-1816  F/G  22/2 


NL 


V  / 


L*  |2j8  |2.5 

US  ■“ 

“  Li  |2_2 

Sr  i3.6  mbs 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANOABOS  - l»*3 -  A 


'AW/itAVv  i aw. 


PREDICTED  AXIAL  TEMPERATURE  GRADIENT 
IN  A  VISCOPLASTIC  UNIAXIAL  BAR  DUE  TO 
THERMOMECHANICAL  COUPLING 


by 


D.H.  Allen 
Assistant  Professor 
Aerospace  Engineering  Department 
Texas  A&M  University 
College  Station,  TX  77843 


4875-84-15 


November  1984 


ABSTRACT 


The  thermomechanical  response  of  a  uniaxial  bar  with  theraoviscoplastic 
constitution  is  predicted  herein  using  the  finite  element  method.  After  a 
brief  review  of  the  governing  field  equations,  variational  principles  are  con¬ 
structed  for  the  one  dimensional  conservation  of  momentum  and  energy  equations. 
These  equations  are  coupled  in  that  the  temperature  field  affects  the  displace¬ 
ments  and  vice  versa. 

Due  to  the  differing  physical  nature  of  the  temperature  and  displacements, 
first  order  and  second  order  elements  are  utilized  for  these  variables,  respec¬ 
tively.  The  resulting  semi-discretized  equations  are  then  discretized  in  time 
using  finite  differencing.  This  is  accomplished  by  Euler's  method,  which  is 
utilized  due  to  the  stiff  nature  of  the  constitutive  equations. 

The  model  is  utilized  in  conjunction  with  stress-strain  relations  devel¬ 
oped  by  Bodner  and  Partom  to  predict  the  axial  temperature  field  din  a  bar  sub¬ 
jected  to  cyclic  mechanical  end  displacements  and  temperature  boundary  condi¬ 
tions.  It  is  found  that  spacial  and  time  variation  of  the  temperature  field 
is  significantly  affected  by  the  boundary  conditions. 

TABLE  OF  SYMBOLS 

t  -  time 

P  -  axial  internal  resultant  force 
Px  -  axial  externally  applied  force  per  unit  length 
x  -  axial  coordinate  dimension 
a  -  axial  stress  component 
A  -  cross-sectional  area 

Tx  -  end  traction  in  units  of  force  per  unit  area 


s  -  surface  area 


Sc  -  area  of  the  longitudinal  surface  of  the  bar 
e  -  axial  strain  component 
u  -  axial  displacement  component 

-  internal  state  variable  representing  axial  inelastic  strain 
G  -  Young's  modulus  in  the  axial  coordinate  direction 

a  -  coefficient  of  thermal  expansion  in  the  axial  coordinate  direction 
T  -  temperature 

T  -  reference  temperature  at  which  no  deformation  is  observed  at  zero  load 
R 

-  internal  state  variable  representing  drag  stress 
q  -  axial  component  of  heat  flux 

k  -  coefficient  of  axial  thermal  conductivity 
Cv  -  specific  heat  at  constant  elastic  strain 
p  -  mass  density 

r  -  internal  heat  source  per  unit  mass 
L  -  length  of  the  bar 

INTRODUCTION 

It  is  well  known  that  mechanical  and  thermodynamic  coupling  are  signif¬ 
icant  in  metallic  solids  [1-11].  The  author  has  recently  developed  a  model 
capable  of  predicting  this  coupling  effect  in  thermoviscoplastic  metals  [12]. 
In  the  previous  paper  a  cyclic  strain  control  loading  on  a  sample  of  IN100  at 
1005°K  (1350°F)  was  used  to  predict  a  temperature  rise  of  approximately  3.7°K 
per  cycle  when  the  strain  amplitude  was  2 %  and  the  specimen  was  adiabatically 
insulated. 

The  focus  of  the  current  research  is  to  consider  the  effect  of  thermal 
boundary  conditions  on  this  same  process.  The  introduction  of  these 


conditions  causes  the  strain  and  temperature  fields  to  be  inhomogeneous  even 
though  the  stress  field  is  homogeneous  if  the  bar  is  prismatic.  This  spacial 
variation  in  the  field  variables  causes  the  process  to  be  difficult  to  model 
because  the  thermomechanical  constitutive  equation*?  are  highly  nonlinear  stiff 
differential  equations.  In  this  paper  the  finite  element  method  is  utilized 
to  spatially  discretize  the  dependent  variables  displacement  and  temperature, 
and  the  finite  difference  method  is  employed  for  timewise  discretization. 

This  process  results  in  a  set  of  highly  nonlinear  algebraic  equations. 

Since  the  thrust  of  this  research  is  to  obtain  accurate  results  without 
regard  to  numerical  efficiency,  the  results  are  obtained  via  the  relatively 
inefficient  but  accurate  method  of  simply  utilizing  successively  smaller  time 
steps  along  with  refined  spatial  mesh  to  obtain  a  convergent  and  therefore 
accurate  solution  for  the  temperature  and  displacement  fields  both  spatially 
and  as  a  function  of  time  for  a  cyclically  imposed  end  displacement. 

The  physical  interest  in  the  problem  is  to  determine  the  effect  of 
temperature  boundary  conditions  on  the  predicted  temperature  rise  in  a  bar  sub¬ 
jected  to  cyclic  mechanical  loading.  It  is  found  from  the  analysis  that  the 
introduction  of  these  nonadiabatic  boundary  conditions  causes  significant  axial 
temperature  gradients.  Since  nonadiabatic  conditions  cannot  be  avoided  in 
experimental  research,  it  is  concluded  that  experimental  tests  of  this  type 
should  be  viewed  with  caution  when  their  purpose  is  to  construct  constitutive 
relations. 


PROBLEM  SOLUTION 

Field  Problem  Description 

The  following  field  equations  are  given: 

a)  equilibrium  [13], 

3P  -  -p_  (x)  , 

3x  x 


(1) 


3 


where  the  axial  resultant  P  is  defined  by 


P 


adA 


and  (2) 


T  ds 
x 


(3) 


b)  strain-displacement  relation 


£ 


3u 

3x 


(4) 


c)  thermomechanical  constitution, 

O  -  E[e  -  ox  -  a (I  -  Tr)]  ,  (5) 

30^ 

— ■  %(e,  T,  a^)  ,  i  -  l,z  ,  and  (6) 

o--kI  >  (7) 

where  z  is  the  total  number  of  internal  state  variables;  and 

d)  conservation  of  energy 

[(E€  -  Ea1  +  EoIr)  ^  +  B»2  I  H]  -  EO,  T  If  -  PCV  |£  -  +  Pt  -  0  .  (8) 

The  conservation  of  mass  is  satisfied  trivially  and  the  second  law  of  thermo¬ 
dynamics  has  been  previously  shown  to  be  satisfied  by  the  above  equations  [14- 
lb].  It  should  be  noted  that  equilibrium  equation  (1)  satisfies  equilibrium 
in  the  axial  coordinate  direction  only  in  an  average  sense  over  the  cross-section. 
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The  above  6+Z  equations  (excluding  definition  (3))  define  a  nonlinear 
initial-boundary  value  problem  (together  with  appropriate  thermal  and  mech¬ 
anical  initial  and  boundary  conditions)  in  which  the  following  dependent  vari¬ 
ables  are  sought  as  functions  of  x  and  t:  o,  £,  u,  q,  T,  P,  and  a^. 

For  convenience  the  domain  is  defined  to  be  of  length  L,  so  that  boundary 
and  initial  conditions  are  of  the  form: 


u (x , 0)  =  Ug  *  known 
T (x, 0)  =  Tg  *  known 


initial  conditions 


(9) 


and 


essential 

boundary 

conditions 


/  u(0,t)  =  u° 

u(L,t)  =  u^ 
C 

|  T(0,t)  =  T° 
^  T(L,t)  = 


known  or  P(o,t)  =  P^ 
known  or  P(L,t)  =  P^ 
known  or  q(0,t)  =  q® 
known  or  q(L,t)  =  q^ 


known  \ 
known  I 


known  I 
known  / 


natural 
boundary 
conditions.  (10) 


It  is  now  assumed  that  a  =  a(x)  so  that  equation  (2)  reduces  to 


P  -  OA 


(11) 


Therefore,  substituting  (4)  into  (5)  and  this  result  into  (11)  gives 


P  -  EA 


[s  -  °i  - a(I  -  v] 


(12) 


The  above  result  is  now  substituted  into  (1)  to  obtain 


(13) 


5 


which  represents  the  differential  equation  relating  displacements  and  temper¬ 
ature  to  the  applied  load  px  (x) . 

Equations  (4)  and  (7)  are  next  substituted  into  energy  balance  law  (8) 
and  this  result  is  integrated  over  the  cross-sectional  area  A  to  obtain 


[(ES-E“i+  “i)-5 

+  Afe(kf)--AOr 


+  Ea2  T 


AEa  T 


v  U  A  V  x 

3t3x  v  Tt 


where  it  has  been  assumed  that  all  field  variables  depend  on  x  and  t  only. 
Now  define 


q  dA 


-/k£ 


dA  -  -k  A 


Careful  inspection  of  equations  (13)  and  (14)  will  indicate  that  these 
equations,  together  with  internal  state  variable  growth  laws  (6)  and  initial 
and  boundary  conditions  (9)  and  (10),  represent  a  well-posed  boundary  value  prob¬ 
lem  in  terms  of  the  2+z  dependent  variables  u,  T,  and  ot^. 

Solution  Procedure 

The  field  problem  is  to  be  solved  analytically  using  the  semi-discretized 
finite  element  technique  with  timewise  finite  differencing.  In  order  to  ac¬ 
complish  this,  differential  equations  (13)  and  (14)  must  first  be  vritten  in 
a  suitable  variational  form. 

Variational  Equations 

Consider  first  equation  (13).  This  governing  equation  is  integrated  against 

a  suitably  smooth  test  function  v  -  v(x)  over  the  domain  of  some  element  fl  : 

x  <  x  <  x  . .  : 
e  .  e+1 


l. 


6 


[s{“[&-“i  v] }  +pk] 


0 


Integrating  by  parts  results  in 


Xe+1 

-  /  v]' 


[3u 

3x  “  al  ' 


a(T  -  Tr) 


v  p  dx 
rx 


Substituting  equation  (12)  into  the  boundary  term  thus  results  in 


-  /"uglfi.,. 

V  " 


a(T  -  Tr)  dx  * 


“v(Xe+l)  P(xe+1)  +  V(xe)  P(xe) 


x 

-f 


v  p  dx 
rx 


Now  consider  equation  (14).  Once  again  the  governing  equation  is  inte¬ 
grated  against  a  suitably  smooth  test  function  w  *  w(x)  over  the  domain  of  the 


element  ft  : 

e 


I  »[A  (Ei-  +  “>)f 


fe  (*£) 


+  Ea  T 


32u  3T 

-***  3^  -  ApCv  3T  +  A 


+  A  pr  >  dx  ■  0 


Integrating  the  heat  flux  term  by  parts  results  in 


3w  3t  /  3u  \  3al  2  3T 

•-kAl$3l  +  wA  (E3f-ECtl+EaTR)3r  +  EaT3T 


•EaT  3B^  fdX  "  W(xe+1}  Q(xe+1>  -  W(xe)  Q(xe} 


fc  I  - 


pr  J  dx 


where  equation  (15)  has  been  substituted  into  the  boundary  terms. 


Finite  Element  Spacial  Discretization 

Quadratic  displacement  and  linear  temperature  fields  are  now  chosen  within 
each  element: 


u(x,t)  -  Z  u*  ,  xg  <  x  <  x  j  ,  and 
i*l 


T(x,t)  *  Z  T?  »  x  <  x  <  x 
i-1  1  6 


e+1 


6  6  6  6 

where  u^  *  u^(t)  and  =  T^(t)  are  the  nodal  displacements  and  temperatures, 
respectively,  and  ^  *  i|/  (x)  and  are  quadratic  and  linear  shape 

functions,  respectively  [17]. 

Appropriately,  v  and  w  are  endowed  with  the  properties  of  u  and  X,  respec¬ 
tively,  so  that 
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(23) 


w  =  <j>  i  -  1,2 
i 

Substitution  of  equations  (12)  and  (21)  through  (23)  into  variational  principle 
(18)  results  in 


The  above  may  be  written  in  the  following  compact  form 


3 

Z 

J-l 


2 

u?  +  Z 


i  -  1,2,3, 


where 


EA 


dx  dx 


dx 


i  *  1,2,3;  j  =  1,2,3 


(25) 


(26) 


d*i  e 

EAa  dx 


i  -  1,2,3;  j  =  1,2 


(27) 


^  d r  (_ai  +  aTR)dx 


-p  (xi) 


dif  p  dx 
Y1  x 


i  =  1,2,3. 


(28 


Similarly,  substitution  of  equations  (21)  through  (23)  into  equation  (20)  re¬ 
sults  in 


Equations  (29)  may  be  written  in  the  following  form: 


3  2  Xe+L  T  /  2 

i,  ru u! +  .i  ru  Tje + .  /  ♦;*  i, 


/  2  \  /  3  due  an  /  2  dT®  \ 


e+l  3a. 

-  /  ♦*  A(-E0l  +  Eaxs)  ^ 


Q(x±)  , 


where 


*e+l  dij/®  3  a. 

K® .  =  /  AE  <fr*  -r-J-  -r-i  dx 

ii  w  Ti  dx  3t 


i  =  1,2;  j  =*  1,2,3;  and 


*e+l  d<j>*  d<j>e 

Ki  s-  J  “irsT1*1  ;  1  ■  !•*  j  - 1.2- 


Finite  Difference  Timewise  Discretization 

Time  dependence  in  equations  (6)  and  (30)  is  handled  via  finite  differ¬ 
encing.  Although  higher  order  approximations  may  be  used,  Euler  forward  dif- 

0  0 

ference  approximations  are  now  entered  for  the  time  rate  of  change  of  a^,  T^, 
and  ue. 


^r~  (x,t)  3  [oj  (x,t  +  At)  -  (x , t) ] /At ,  k  =  1,  ...,z 


(t)  s  [T*  (t  +  At) 
dt  m 


( t) ] /At , 

in 


m  =  1 » 2  » 


and  (34) 


ar  (t)  3  [Um  (t  +  At) 


u'  (t)]/At, 

m 


n  -  1,2,3. 
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Substitution  of  (33)  through  (35)  into  finite  element  equations  (30)  gives 


r 


li 


I 


t 


|r 


\ 


Z  K?.  uf  +  Z  S?,  T® 
j=l  3  j=l  «  3 


y»  e+i 
+  /  A<j) 

X 


H2  1  f  2  /Te(t  +  At)  -  Te(t)  \  ] 

,;,'W]U.<'  ■■  '  >*-J 


„  T  l  e  el  T  1  +  ^ 

[A  j  jj  [jA 


r  2  /Ot+AC)  -i^  ] 
- 3« - )  Vj 


+  pr  >■  dx 


/e+1  ou. 

A4>®  t-Ea1  (t)  +  EaTR]  (t) 


3a, 


-Q<xi) 


i  *  1 , 2 


The  above  may  be  written  as  follows: 


z  r  ue  +  z  r 

j*l  x3  J  j  =  1  iJ  J 


2  2 


+  I  Z  C,  T®  if  +  Z  D,,  Te 

k«l  -,*1  &  j  k  j=1  ij  3 


2  3  2 

+  Z  Z  E,.  ,  if  u®  +  I  G,.  T® 
k-1  j-1  ikj  k  J  j-1  ij  j 


+ -u  tj  ■  Fi 


i  -  1,2, 
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(36) 


(37) 


where 


/e+1 

. i e  Ea2  ,e  ,e  . 

A*i  IT  *k  d*  • 


i  =  1,2;  j  =  1,2;  k  =  1,2  ,  (38) 


~e+l 

Dil  "  "  /  A4>i  ft"  fTt(c)(4,l)2  +  T2(t)  44]  dx’  i  =  1.2  ,  09) 


~e+l 

=  -  J  A d>*  ~  [T®(t)<D®4)*  +  Tj(t)  (4>|)2]  dx. 


i  -  1,2  ,  (40) 


Eikj  “ 


e+l 

I  A*e  Ett 

■  j  It 


\  ar  dx  • 


i  -  1,2,;  k  =  1,2;  j  *  1,2,3 .  (41) 


.  e  Ea  e-v  ,e 

"i  aF  [ulCt)  *1 


.  *♦!  .  e  s*2  . 

1  3*  +  u2(t)  *1  3*  +  u3(t)  *1 


e  »*3  1  . 

1  3^J  dx> 


1,2  ,  ( 


/Xe+1  r  3dje  3^0  3^1 

<  ft  4^4  IT +  4^4  *r  +  ^>*2  zr\ dx  * 

x  *- 

e 

«  e 

=  -  J  A4>®  J—  <j)j  dx  ,  i  *  1,2;  j  =  1,2; 

xe 

;  -  /  <  sr  ( j,  T3e<t)*l)+  Dr  d* 

X  ■»  J 


i  -  1,2  ,  (43) 


and  (44) 


x  . , 

/e+l 

A 4>®  [-EOj 


(t)  +  EaTR]  (t)  dx 


-Q(xJ 


1,2  . 


Equations  (37)  may  be  written  equivalently  as  follows: 


Z*ii 

j=l  J 


u, 


2 

+  Z 
j-1 


’ij 


n 


6  — 0 

where  and  are  as  defined  previously,  and 


2 

+  Z 
k-1 


Eikj  Tk 


(46) 


and  (47) 


2 

+  Z 

k=l 


C,  ..T.®  +  D.  . 
ijk  k  ij 


+  G, .  +  H 


ij 


(48) 


The  above  equations  may  be  adjoined  with  equations  (25)  to  obtain  the  following 
set  of  nonlinear  equations  for  each  element. 


i 

f 

i 


i 


5x1 


(49) 


where  all  nonlinearity  is  contained  in  [S],  {Fe},  and  {F®}. 


Global  Assembly  and  Boundary  Conditions 

Global  assembly  is  accomplished  in  the  standard  way  using  the  Boolean 
matrix  [17].  Interelement  continuity  is  guaranteed  by  setting 


'2 

r 

-  0 

I 

and  (50) 

Y2 

.♦r1 

-  0 

• 

(51) 

Boundary  conditions  are  implemented  in  the  standard  way:  1)  essential 
boundary  conditions  are  handled  by  placing  one  on  the  diagonal  of  the 


appropriate  row  and  zeros  off  diagonal  in  the  stiffness  matrix,  and  the  speci¬ 
fied  value  of  the  essential  variable  on  the  right  hand  side;  and  2)  natural 
boundary  conditions  are  implemented  directly  to  the  right  hand  side. 


Solution  of  the  Nonlinear  Algebraic  System 

Initial  conditions  are  used  for  the  first  time  step.  The  time  step  At 
is  supplied  for  each  load  increment  and  boundary  conditions  are  incremented 
directly  from  supplied  input  functions. 

The  internal  state  variable  is  handled  in  equations  (23)  and  (45)  by 
using  equations  (35).  is  initialized  according  to  reference  18.  The  non¬ 
linear  stiffness  matrix  [S]  is  initialized  using  nodal  temperatures  and  displace¬ 
ments  from  the  previous  time  step.  The  displacements  and  temperatures  at  time 
t+At  are  then  estimated  directly  and  without  iteration  by  utilizing  equations  (49) 
for  very  small  time  steps. 
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EXAMPLE  PROBLEMS 


In  order  to  completely  define  an  example  problem  it  is  necessary 
to  specify  internal  state  variable  growth  laws  (6).  Numerous  models 
have  been  proposed  for  crystalline  metals  [18,19].  Since  it  is  not  the 
purpose  of  this  research  to  compare  these  models,  a  relatively 
established  model  proposed  by  Bodner  and  Partom  [20]  has  been  chosen. 
This  model  contains  two  internal  state  variables:  the  inelastic  strain 
(o^)  and  the  drag  stress  (a2).  The  growth  laws  for  these  variables  are 
given  by 

•  2  D  cr 

a  i  "  VT  o  ]^]  exP 

and 

•  •  /  a2  -  Z  \r 

a  2  -  m<zi  -  <*2)  ,  (53) 

where  DQ,  n,  m,  Z  Z^,  and  r  are  experimentally  determined  material 
constants . 

For  the  purpose  of  modeling  the  temperature  gradient  in  a  specific 
component,  a  hypothetical  problem  has  been  chosen  using  material 
properties  representative  of  Inconel  100  at  1005°k  ( 1 350  °F) .  The 
material  and  geometric  properties  are  given  in  Table  I.  The  geometry  is 
representative  of  a  cylindrical  uniaxial  bar  which  is  2.50  inches  long 
and  0.25  inches  in  diameter.  It  has  previously  been  shown  that  Bodner 
and  Partom' s  model  accurately  predicts  the  stress-strain  behavior  of 
INI 00  under  uniaxial  loading  conditions  for  both  monotonic  and  cyclic 
strain  controlled  loadings  [12,18]. 
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5.032  MPa/°K 


PC 

v 

a  -  13.14  x  10-6  in/in/°K 
—6  2 

k  ■  21.0  x  10  MPa jn 

sec4K 

E  -  146.86  x  103  MPa 

A  -  7.12557  x  10~5  m2 

T  -  1005 °K 

K 

L  *  .06350  m 

Dq  ■  10  x  103  in/in 

n  -  0.70 
m  -  2.57 
Z1  -  1015.0 

Zj  -  600.0 

r  -  2.66 
*  0 . 0019 


Table  Material  Properties  for  IN  100  at  1005°K  (1350°F) 


Utilizing  the  material  properties  described  above,  the  following 
effects  have  been  studied  using  the  model  developed  herein: 

1)  the  effect  of  variation  of  strain  rate  on  the  time  dependence 
of  temperature  at  the  midpoint  of  a  monotonically  extended  uniaxial  bar 
which  is  insulated  on  the  longitudinal  boundaries  [Figs.  1-4]; 

2)  the  spacial  variation  of  temperature  for  the  case  above  [Fig. 
51; 

and 

3)  the  effect  of  end  temperature  boundary  conditions  on  the 
temperature  at  the  center  of  a  uniaxial  bar  which  is  held  at  fixed 
temperature  at  the  end  points  and  subjected  to  cyclically  imposed  end 
displacements  [Figs.  6  and  7]. 

Examples  1  and  2  are  constructed  primarily  to  determine  the  effects 
of  thermomechanical  heating  on  the  stress-strain  behavior  of  uniaxial 
constitutive  specimens.  It  is  found  in  examples  1  and  2  that  if  a 
specimen  is  mounted  in  an  experimental  apparatus  which  has  massive  grips 
simulating  a  fixed  temperature  boundary  condition  there  can  be 
substantial  axial  temperature  gradients  induced  in  a  time  dependent 
boundary  layer  near  the  ends  of  the  specimen.  On  the  other  hand,  these 
boundary  conditions  do  not  appear  to  substantially  affect  the  predicted 
stress-strain  behavior,  especially  when  the  strain  measurement  is  taken 
between  the  thermal  boundary  layers  near  the  grips.  Therefore,  it  would 
appear  that  the  standard  procedure  for  obtaining  stresses  and  strains  in 
uniaxial  bars  is  not  substantially  affected  by  mechanically  induced, 
axial  temperature  gradients  when  the  grips  are  at  fixed  temperature 
equivalent  to  the  Initial  specimen  temperature  and  the  bar  is  loaded 
monotonicalaly.  However,  it  should  be  noted  that  massive  grips  which 
are  mounted  outside  a  furnace  could,  by  their  much  lower  temperature 
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than  the  initial  speciman  temperature,  may  induce  significant  error  in  predicted 
strains  if  the  strain  is  measured  by  dividing  relative  displacement  by  some 
gage  length. 

The  final  example  demonstrates  that  under  cyclic  loading  conditions  the 
above  conclusions  may  not  necessarily  be  true,  especially  when  the  specimen  is 
subjected  to  high-cycle  fatigue  and  at  high  strain  rates.  There  is  definitely 
a  trend  towards  an  increasing  mean  temperature  in  the  bar,  and  this  mean  tempera¬ 
ture  is  strongly  affected  by  the  thermal  boundary  conditions  as  well  as  the 
loading  rate.  Although  it  would  be  interesting  to  determine  the  mean  temperature 
rise  in  a  cyclic  fatigue  test,  the  current  algorithm  precludes  this  analysis  due  to 
the  extremely  large  computer  times  necessary  to  predict  only  a  few  cycles  of 
response  (approximately  43.8  CPU  minutes  on  an  Amdahl  470/V6  for  the  example 
demonstrated  in  Figs.  6  and  7). 

Example  3  also  demonstrates  another  interesting  phenomenon  which  may  be 
significant  in  large  space  structures.  If  the  bar  is  perfectly  insulated  the 
mean  temperature  rise  per  cycle  for  the  relatively  slow  loading  rate  shown  in 
Fig.  6  is  3.7°K,  whereas  if  the  ends  of  the  bar  are  held  at  a  fixed  temperature  of 
1005°K,  the  mean  rise  is  1 .0°K  per  cycle.  Faster  loading  rates  show  less  difference 
between  the  adiabatic  result  and  the  fixed  end  temperature  result.  Since  many 
of  these  structures  are  expected  to  be  extremely  flexible  truss-like  configurations, 
a  typical  metallic  member  which  undergoes  some  yielding  (which  might  be  desirable 
in  order  to  induce  natural  damping)  might  in  fact  undergo  substantial  enough 
heating  during  vibrational  response  such  that  the  material  properties  could  be 
adversely  affected,  thus  resulting  in  a  material  related  failure  of  the  structure. 
However,  further  investigation  is  needed  on  this  last  Issue  since  it  is  expected 
that  the  primary  form  of  heat  flux  off  of  space  structures  will  be  via  radiation 
on  the  longitudinal  surfaces  of  the  truss  member.  Since  the  current  analysis 
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has  treated  these  surfaces  as  insulated,  no  general  statements  can  be  made  at 
this  time  regarding  thermomechanical  heating  in  space  structures. 
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CONCLUSION 


t 


The  current  research  has  attemped  to  demonstrate  the  effects  of 
mechanical  loading  on  one-dimensional  temperature  gradients  in  a  class 
of  viscoplastic  media.  Due  to  the  nonlinearity  and  stiffness  of  the 
field  equations,  it  was  necessary  to  utilize  a  numerical  algorithm. 
This  algorithm  has  been  shown  to  be  very  inefficient  for  solving  even 
one-dimensional  examples.  Therefore,  it  is  apparent  that  significant 
refinement  of  the  procedure  will  be  necessary  before  multi-dimensional 
analyses  can  be  performed  by  this  method.  Specifically,  it  would  be 
significant  to  determine  the  effect  of  transverse  temperature  gradients 
on  the  stress-strain  behavior  of  constitutive  specimens.  Furthermore, 
the  effects  of  thermal  boundary  conditions  on  the  longitudinal  surface 
needs  attention.  The  author  is  currently  studying  a  perturbation 
technique  for  more  efficient  solution  of  these  issues. 

The  above  points  notwithstanding,  the  current  research  demonstrates 
some  important  results.  These  are: 

1 )  The  axial  temperature  gradient  in  a  viscoplastic  uniaxial  bar 
is  strongly  affected  by  the  thermal  boundary  conditions  on  the  ends. 

2)  The  end  temperature  boundary  conditions  can  cause  temperature 
gradients  which  are  substantial  enough  to  induce  spacial  variations  in 
stress  and  strains  which  invalidate  the  standard  procedure  of  using 
average  quantities,  although  when  grips  are  mounted  within  a  furnace  at 
spacially  constant  temperature,  it  appears  that  the  standard  procedure 
is  accurate. 
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3)  There  is  a  trend  toward  increasing  average  temperature  in 
cyclically  loaded  bars;  whether  or  not  this  effect  is  significant  is 


I 

strongly  dependent  on  the  thermal  boundary  conditions  and  the  loading 

S  rate. 
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Abstract  s  -  surface  area 


This  paper  focuses  on  the  effect  of 
theraoaechanlcally  Induced  heating  on  the 
response  of  a  single  member  of  a  space  truss 
structure  which  behaves  viscoplastlcally.  The 
governing  equations  are  given  for  a  typical  truss 
member,  wherein  material  Inelasticity  is 
reflected  in  constitutive  equations  via  a  set  of 
internal  state  variables,  each  characterized  by  a 
history  dependent  growth  law.  The  governing 
equations  are  coupled  In  the  sense  that 
temperature  and  displacement  are  dependent  on 
each  other.  This  difficulty,  together  with  the 
fact  that  the  Inelastic  constitutive  equations 
are  nonlinear  and  numerically  stiff,  requires 
that  a  computationally  complex  semldiscretized 
finite  element  spatial  technique  be  utilized  to 
obtain  a  solution.  This  procedure,  detailed 
herein,  la  utilized  to  predict  the  response  of  a 
typical  metallic  space  truss  member  under 
vibrational  or  cyclic  loading.  Particular 
Interest  Is  placed  on  the  temperature  rise  In 
such  a  member  due  to  hysteretlc  loss  during 
structural  vibrations  and  In  the  presence  of 
complex  thermal  boundary  conditions 
representative  of  space  conditions.  Example 
cases  are  constructed  for  a  typical  cylindrical 
bar  of  6061 -T6  aluminum  both  with  and  without 
special  coatings.  Results  Indicate  that 
significant,  possibly  even  catastrophic,  heating 
can  occur  due  to  thermomechanlcal  coupling. 


Nomenclature 


t  -  time 

P  -  axial  Internal  resultant  force 
p  -  axial  externally  applied  force  per  unit 
length 

x  -  axial  coordinate  dimension 
o  -  axial  stress  component 
A  -  cross-sectional  area 

T^  -  end  traction  in  units  of  force  per  unit  area 
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S0  -  area  of  the  longitudinal  surface  of  the  bar 
e  -  axial  strain  component 
u  -  axial  displacement  component 

“  Internal  state  variable  representing  axial 
Inelastic  strain 

E  -  Young's  modulus  in  the  axial  coordinate 
direction 

a  -  coefficient  of  thermal  expansion  In  the 
axial  coordinate  direction 
T  -  temperature 

Tr  -  reference  temperature  at  which  no 
deformation  is  observed  at  zero 
load 

Oj  -  internal  state  variable  representing  drag 
stress 

q  -  heat  flux  vector 
q  -  axial  component  of  heat  flux 
k  -  coefficient  of  axial  thermal  conductivity 
Cy  -  specific  heat  at  constant  elastic  strain 
p  -  mass  density 

r  -  Internal  heat  source  per  unit  mass 
L  -  length  of  the  structural  element 
Dq,  n,  m,  Z1 ,  Zj,  ZQ,  r  -  material  constant  used 
in  Bodner  and  Partom's  model1 
qo  -  flux  on  longitudinal  boundary 
e  -  absorbing  portion  of  perimeter  of  an  element 
normal  to  longitudinal  axis 
qs  -  solar  radiation  flux 

qE  -  earth  radiation  flux 


l 


<*3  -  absorptivity 
Fg  -  earth  radiation  view  factor 
\9  -  Incident  angle  of  solar  radiation  on 
structural  component 

1E  -  Incident  angle  of  earth  radiation  on 
structural  component 

o3  -  Stefan-Boltzmann  constant  »  5.775  x  10 
MPa  m/sec/CK)4 
Tq  -  deep  space  temperature 

Introduction 


conditions  on  the  longitudinal  surface  of  the 
truss  component  are  significant  enough  to  carry 
off  all  heat  generated  due  to  hysteretlc  loss. 

The  paper  first  reviews  the  governing  field 
equations,  then  briefly  discusses  the  procedure 
whereby  a  numerical  algorithm  Is  constructed  for 
modeling  the  problem.  This  Is  followed  by  a 
detailed  discussion  of  the  implementation  of 
thermal  boundary  conditions.  Finally,  example 
results  are  obtained  for  representative  space 
structural  components. 


Governing  Field  Equations 


The  governing  field  equations  were  presented 
In  reference  5  for  quasl-statle  conditions.  For 
problems  Involving  Inertial  effects,  the 
governing  equations  are  as  follows: 

a)  equilibrium^' 


It  is  well  known  that  In  vlacoplastic  metals 
a  certain  amount  of  mechanically  Induced 
hysteretlc  mechanic  energy  loss  Is  converted  to 
heat,  thus  resulting  In  a  temperature  rise  in  the 
medium.  In  recent  research^’ J  a  model  has  been 
developed  for  predicting  this  effect  by  utilizing 
thermodymanic  constraints  together  with 
constitutive  equations  of  Internal  state  variable 
type4.  Furthermore,  It  has  been  shown  that  In  a 
perfectly  insulated  uniaxial  bar  ,  as  well  as  In 
a  uniaxial  bar  with  Insulated  longitudinal 
surface  and  fixed  end  temperature5,  significant 
temperature  rise  can  occur  in  the  component 
during  cyclic  loading. 

The  purpose  of  the  current  research  la  to 
simulate  the  response  of  a  typical  metallic  space 
truss  structural  element  (see  Fig.  1)  In  the 
postylelded  state  and  to  determine  if  significant 
heating  occurs  when  this  component  Is  subjected 
to  cyclic  mechanical  loading.  This  problem  Is  of 
Interest  because  a  certain  amount  of  material 
inelasticity  Is  desirable  in  order  to  produce 
passive  structural  damping.  The  factors  of 
Interest  in  this  simulation  are  the  effects  of 
thermal  boundary  conditions  and  loading  rate  on 
the  thermal  response.  In  particular,  It  is  of 
Interest  to  determine  If  radiative  boundary 


Fig.  1.  Typical  Space  Truss  Structural  Element. 


1Z.  “  p,(x)  (!) 

3x  * 

where  the  axial  resultant  P  is  defined  by 


c 

bl  strain-displacement  relation 

c  -  3u  (4) 

3x 


c)  thermomechanlcal  constitution, 


(5) 

a."  2  D  a  exp  [-/n4-l\/a2\  2n"| 

lvf°w  J 

(6) 

i2-m(Z1-a2)oi1-A1Z1 

(7) 

q»-k  .3T 

(8) 

3x 


where  o.  and  o_  are  the  internal  state  variables 
( ISV)  representing  inelastic  strain  and  drag 
stress,  respectively,  In  the  constitutive  model 
developed  by  Bodner  and  Partom  .  Several  other 
constitutive  models  have  been  developed  for 
vlscoplastlc  metals,  and  these  are  reviewed  In 
references  1  and  8.  Finally, 

d)  conservation  of  energy^’5, 

[  (Ee-Ea.+EaT_)3al  +  Ea2T3T]-EaT3e  -pc  3T-3q+pr-0 
1  R  3t  3C  3t  3t  Jx 

(9) 

Conservation  of  mass  is  satisfied  trivially 
(under  the  assumption  of  small  motions  In  a 
closed  system) ,  and  the  second  law  of 
thermodynamics  has  been  previous ^y  shown  to  be 
satisfied  by  the  above  equations^’  . 

The  governing  equations  are  adjoine^  with 
appropriate  Initial  and  boundary  conditions  such 
that  a  well-posed  boundary  value  problem  Is 
constructed  In  terms  of  the  following  dependent 
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varlabl«a  which  are  sought  as  functions  of  x  and 
t:  o,e,u,q,T,P,a. ,  and  a,.  Due  to  ISV  growth 
laws  (6)  and  (7;  (as  well  as  radiative  boundary 
conditions),  the  problem  is  nonlinear. 


1  e  e 

u(x,t)«  (OVjOO  *e<  x<3te^.j 

T(x,t)-|  T*  (t)$*(x)  xe<x<xe+1 


(15) 

(16) 


Solution  Procedure 


As  described  in  detail  In  reference  5  for 
the  quaslstatic  problem,  the  solution  is  obtained 
using  the  semi -discretized  finite  element 
technique,  wherein  finite  elements  are 
constructed  spatially,  and  finite  differencing  la 
used  in  time.  The  result  is  a  time  marching 
algorithm  which  is  reviewed  here  briefly. 

First,  equations  (4)  and  (5)  are  substituted 
Into  (2)  and  this  result  Is  substituted  Into  (1) 
to  give  the  following  equllbrlum  equation: 


Ht- 

Ol-ad-Vll 

l  1.3* 

Ji 

-  Px<*> 


(10) 


Next,  equation  (4)  is  substituted  Into  equation 
(9)  to  obtain  the  coupled  energy  balance  law: 


BE3u-Ea.+EaT  ^  +Ea2T3Tl  -EaT32u 

,  3x  ♦  /  3t  3t  J  3t3x 


(11) 


-PC  3T  -  V-q  +pr  «  0 
V  3t 

The  result  is  a  set  of  two  coupled  partial 
differential  equations  In  terms  of  axial 
displacement  u»u(x,t)  and  temperature  T«T(x,t). 

Variational  Principles  and  Finite  Element 

Discretization 

Selecting  a  suitably  smooth  test  function 
v-v(x)  over  the  domain  of  some  element  Q_ 

X  <  X<  X 

e  e  . 
variations 


.  ,  one  may  construct  the  following 
■Tal  principle  from  equation  (12): 


■/**““  ‘-"VI 


dx 


-v(x^i)P(x 


e+1 


)+v(xJP(xJ-4  efl  VP  d* 


(12) 


where  the  boundary  terms  result  from  the  standard 
integration  by  parts. 

The  variational  principle  for  heat  equation 
(11  )  is  constructed  by  first  Integrating  this 
equation  against  a  test  function  w-w(x)  on  Qa  to 
obtain 


w|[  (Ee-EaI-tEaTIl)3aH-Ea2T  3T  J- 


EaT  3  u 
3x3 1 


e 

-PC., 


3T  +  V-q  +  pr  l  dV 


3t 


}' 


(13) 


Integrating  the  flux  term  by  parts,  assuming  that 
nonaxlal  components  of  flux  are  negligible,  and 
substituting  equation  (8)  will  thus  result  In 

I  w  )a  fE3u-Ea,+EaTD\  3“l  +  Ea2T  3T  1 

I  U\3*  R;~ 

£T  (dx 

3M  x 

/•  e+ 1 

-  I  c  w  q  dx 

lxh  H '  ") d* 


e  2 

-AEaT  3  u  -  ApC  3T  -  IcA  3w 
3t3x  v3t  3x 


-w(Xe+i,Aq(Xe+i)  + 


w(x@)Aq(xe)- 


(14) 


+  /  wA  |pC 

Variational  equations  (12)  and  (14)  are  now 
discretized  by  assuming  the  following 
displacement  and  temperature  fields  in  a  typical 
element  (superscripted  e): 


where  u®  and  T®  are  nodal  displacements  and 
temperatures,  respectively,  and  *®  and  $®  are 

quadratic  .and  linear  shape  functions, 

respectively’.  Furthermore,  v  and  w  are  endowed 
with  the  properties  of  u  and  T.  Note  that  a 
higher  order  element  must  be  used  for 
displacement  than  temperature  due  to  the  fact 
that  temperature  produces  strain  rather  than 
displacement. 

Timewise  discretization  is  Implemented  via 
the  following  backward  finite  difference 
equations: 

dTe  (t)s(Te(t)-Te(t-4t)]/At  Of  1,2  (17) 

m  m  m 

dt 

due  (t)=[ue(t)-ue(t-At)]/At  m-1,2,3  (18) 

m  mm 

“dt 


The  above  equations  require  small  time  steps  in 
order  to  guarantee  numerical  accuracy.  However, 
they  are  uncondltonally  stable  which  is  necessary 
because  IS7  growth  laws  (6)  and  (7)  are 
numerically  stiff1  . 

Substitution  of  equations  (15)  through  (18) 
into  the  governing  field  equations  in  variational 
form  will  result  in  the  following  algebraic 
equations: 


2x3  2x2 

•  ^  w  ^  ' 

5x5 

where  [Ke], 
described  In 


(fl-  W 

5x1  5x1 


[Se], 
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[seL 


and 


(Fe) 


(19) 


are  as 


reference  5,  and 


fyi 

-  9C 

*  X 


dx 


(20) 


where  Is  as  defined  in  reference  5.  The  last 
term  in  the  above  equation  accounts  for  thermal 
flux  boundary  conditions  on  the  longitudinal 
surface  of  an  element. 

After  global  assembly  and  imposition  of 
boundary  conditions  equations  (19)  can  be  solved 
In  a  time  marching  scheme  In  order  to  obtain  the 
nodal  displacements  and  temperatures  as  functions 
of  time. 

Global  assembly  of  the  element  equations  is 
accomplished  in  the  standard  way  using  the 
Boolean  matrix’. 


Imposition  of  Boundary  Conditions 

For  a  typical  space  truss  structural 
element,  the  boundary  conditions  are  assumed  to 
be  of  the  following  type: 


u(0,t)  •  u°  -  known 
u(L,t)  "  -  known 
T(0,t)  «  T°  -  known 
T(L,t)  •  T^  •  known 
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and 


11 


qc~as[qs  COS  X8+FEa_aE)qsc08Xa+FEqEC0#Xe] 

+  a#c(T4-TD4)  (22) 

where  the  first  term  In  the  above  equation  Is  the 
solar  radiation  flux  absorbed  by  the  body,  the 
second  term  is  the  solar  radiation  flux  reflected 
by  the  earth  and  absorbed  by  the  body,  the  third 
term  is  the  earth  radiation  flux  absorbed  by  the 
body,  and  the  last  term  Is  the  flux  radiated  by 
the  member  to  apace. 

The  above  boundary  condltons  may  be 
implemented  to  the  ittscretlzed  global  equations 
in  the  standard  way  .  Although  equation  (22) 
technically  includes  the  unknown  temperature 
field,  the  component  temperature  is  treated  as  a 
known  quantity  In  this  term  for  each  time  step. 
This  approximation  is  acceptable  due  to  the  fact 
that  the  numerical  stiffness  of  constitutive 
equations  (6)  and  (7)  requires  extremely  small 
time  steps  In  order  to  obtain  an  accurate 
solution.  _  . 

EXAMPLE  PROBLEMS 

A  typical  structural  element  has  been 
modeled  with  properties  shown  In  TABLE  1 .  The 
material  properties  were  obtained  experimentally 
in  the  Mechanics  and  Materials  Center  at  Texas 
A&M  University12  for  A1  5086  at  room  temperature, 
which  is  similar  to  A1  6061-T6. 

Sample  cases  were  constructed  for  various 
cyclic  loading  rates  for  two  different  sets  of 
thermal  boundary  conditions,  as  described  In 
TABLE  2.  Both  cases  are  considered  to  be  "worst 
cases"  in  that  the  component  Is  In  a  maximum 
radiation  flux  condition  at  the  maximum 
equilibrium  temperature  during  one  orbital  cycle. 
The  two  cases  differ  in  the  emnissivlty  and 
absorptivity  values  for  the  component  due  to 
differences  In  surface  treatment  of  the 
component.  For  case  I,  the  component  is 
anodized,  and  for  case  II,  the  component  Is 
painted  with  high  emmlsslvlty  ITTRE-S1 3GL0  white 
paint13. 

He  now  consider  two  elements  In  a  large 
space  structure  (see  Fig.  1).  Both  elements  are 
constructed  of  the  same  material  and  are 
geometrically  Identical.  However,  element  one  is 
painted  with  the  high  emmlsslvlty  paint  described 
above  and  Is  in  full  view  of  both  earth  and  sun, 
whereas  element  two  Is  anodized  and  Is  in  view  of 
earth  only.  For  this  case,  as  described  in  Table 
2,  the  components  have  identical  equilibrium 
temperatures  TP  -  295*K(obtalned  by  setting  q  «0 
in  equation  (22*)  )• 

In  both  cases  the  structural  members  have 
been  subjected  to  50  cycles  of  loading  at  three 
different  frequencies:  1  Hz,  5  Hz,  and  25  Hz. 
These  frequencies  have  been  selected  as 
representative  of  resonant  frequencies  in  a 
representative  space  structure.  For  example,  a 
typical  structure  analyzed  In  reference  14  has 
resonant  frequencies  of  4.1  Hz  and  3.4  Hz  in  the 
first  two  modes.  Because  the  resonant  frequency 
of  the  first  mode  In  the  structural  element 
Itself  is  240  Hz,  Inertial  effects  may  be 
neglected  in  these  examples. 

Results  for  the  cases  described  above  are 
shown  lri  Figs.  2  through  8.  In  Figs.  2  through  4 
the  cyclic  stress-strain  curve  is  shown  at  the 


4 


STRESS  (KSI) 


Location  x-L/2  for  CASE’  I  and  at  all  three 
loading  rates.  It  Is  found  that  In  all  cases  the 
specimen  reaches  cyclic  saturation  after 
approximately  five  cycles.  Thereafter,  the 
hysteretlc  energy  loss  per  cycle  becomes  a 
constant  value. 

In  figs.  5  through  7  the  temperature  rise  Is 
plotted  for  both  cases  at  all  three  loading 
rates.  As  expected,  the  amount  of  temperature 
rise  Increases  dramatically  with  loading  rate. 
For  example,  after  50  cycles  the  total 
temperature  rise  at  x-L/2  Is  17.5*K(1  Hz), 

62.5°K(5  Hz),  and  m.7*k  (25  Hz)  for  case  I.' 
Furthermore,  It  Is  apparent  that  while  neither 
surface  treatment  can  be  regarded  as  resulting  In 
negligible  heating;  at  the  higher  loading  rates 
the  anodized  surface  treatment  produces 
temperature  rises  which  are  significantly  higher 
than  those  where  the  surface  is  painted  with 
ITTRE-S13GL0  paint.  Finally,  It  Is  believed  by 
these  researchers  that  the  nonlinear  nature  of 
the  average  temperature  rise  per  cycle  suggests 
that  the  temperature  rise  asymptotically 
approaches  some  upper  bound,  although  this  belief 
cannot  be  corroborated  at  this  time  due  to  the 
large  computer  times  required  in  the  current 
algorithm. 

Fig.  8  shows  that  the  spatial  temperature 
variation  at  5  Hz  Is  approximately  spatially 
homogeneous.  Apparently,  a  very  thin  boundary 
layer  forms  near  the  end  of  the  component,  and 
this  boundary  layer  has  little  effect  on  the 
temperature  at  x-L/2.  In  fact,  subsequent 
investigations  by  the  authors  have  shown  that,  at 
least  for  the  geometry  and  physical  conditions 
considered  herein,  identical  results  may  be 
obtained  more  efficiently  by  neglecting  spatial 
variations  in  displacement  and  temperature. 


Fig-  3.  Cyclic  Stress-Strain  Curve  at  x-L/2  for 
Case  I  Coating  Loaded  at  5  Hz. 


■M.l  1 - 

•MW 


a. am 

strain 


Fig.  4.  Cyclic  Stress-Strain  Curve  at  x-L/2  for 
Case  I  Coating  Loaded  at  25  Hz. 


Fig.  2.  Cyclic  Stress-Strain  Curve  at  x-L/2  for 
Case  I  Coating  Loaded  at  1  Hz. 
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Conclusion 


The  current  research  has  attempted  to 
predict  the  response  of  a  typical  space 
structural  element  which  is  viscoplastic  and  is 
subjected  to  various  cyclic  loading  conditions  in 
the  presence  of  radiation  boundary  conditions. 
Several  general  conclusions  can  be  made  as  a 
result  of  this  reasearch: 

1 )  significant  temperature  rises  may  occur 
due  to  hysteretlc  loss,  although  the  precise 
amount  depends  on  loading  rate  and  surface 
treatment ; 

2)  the  special  paint  ITTRE-S13GL0  appears  to 
produce  significantly  lower  temperature  rises 
than  anodized  surface  treatment; 

3)  the  temperature  rise  appears  to  be 
approaching  an  upper  bound  which  is  loading  rate 
and  surface  treatment  dependent;  and 

4)  the  thermal  boundary  layer  which  forms 
near  the  end  of  the  member  appears  to  have  little 
effect  on  the  far-field  temperature  rise. 

These  conclusions  indicate  that  future 
research  on  this  subject  should  perhaps 
concentrate  on  spatial  variations  in  the  radial 
direction  rather  than  the  axial  direction.  More 
Importantly,  these  results  indicate  that  an 
inelastic  structural  component  may  undergo 
temperature  rises  during  structural  vibrations 
which  are  so  substantial  that  the  material 
properties  of  the  component  may  be  further 
degraded,  thus  leading  to  failure  of  the 
component  and  perhaps  even  failure  of  the  entire 
structure. 
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ABSTRACT 


Under  the  assumption  that  the  stress  is  a  known  function  of  time,  the 
equations  governing  the  behavior  of  a  homogeneous,  thermoviscoplastic  uniaxial 
bar  can  be  reduced  to  a  single  linear  partial  differential  equation  for  the 
temperature  in  terms  of  stress.  This  is  solved  by  means  of  Laplace  transforms. 
Expressing  the  strain  in  terms  of  stress  leads  to  an  explicit  compliance  relation¬ 
ship  for  Bodner's  model.  The  spatial  dependence  of  the  temperature  can  be  re¬ 
covered  from  the  spatially  homogeneous  case  by  a  convolution  in  time  with  an 
appropriate  heat  kernel. 


Introduction 


A  general  model  for  describing  the  time-dependent  inelastic  behavior  of  a 
thermo-viscoplastic  material  has  been  derived  by  D.H.  Allen  in  [1]. 

The  differential  equations  are  highly  nonlinear,  and  stiff  with  respect 
to  the  time  variable.  Consequently,  a  great  deal  of  computational  effort  is 
needed  to  integrate  the  equations  forward  in  time.  Alternate  methods  for 
solving  the  equations  are  therefore  desirable  from  a  standpoint  of  numerical 
efficiency. 

In  the  special  case  where  the  stress,  a,  is  a  known  function  of  time, 
it  can  be  shown  that  the  equations  reduce  to  a  single  linear  partial 
differential  equation. 

The  modified  tempr  ature 

T*(x,t)  =  exp  (§«UU(T(x,t>-T.) 

P  V 

which  we  introduce  in  this  report,  satisfies  the  heat  equation 

LT*-eZ4  T*  =  F[o(t)].  (1.19) 

at  7 

ax 

In  the  spatially  homogeneous  case,  with  insulated  boundaries,  the  modified 
temperature  satisfies 

|t  T*  =  F[o(t)]. 

The  solution  of  the  partial  differential  equation  (1.19)  is  related  to 

the  solution  of  the  ordinary  differential  equation  above  by 

T  *  (x,t)  =  G(x,t)*  T  *  (f) 

p.d.e.  o.d.e. 

G(x,t)  is  an  appropriate  heat  kernel  incorporating  the  boundary  conditions 
of  the  problem,  and  the  convolution  is  with  respect  to  time. 
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The  spatially  dependent  temperature  solution  can  be  recovered  from  the 
spatially  uniform  case  given  the  kernel  G(x,t),  which  is  derived  in  this  report. 

The  term  e  in  (1.19)  is  typically  on  the  order  of  10  in  nondimension- 
alized  units,  and  causes  a  boundary  layer  to  form  near  the  ends  of  the  bar. 

Away  from  the  ends  the  temperature  is  very  nearly  spatially  uniform.  [2]. 


1.  Stress  Decoupling  in  Bodner's  Model 


As  derived  by  D.H.  Allen  in  [1],  the  coupled  set  of  partial  differential 
equations  describing  a  homogeneous,  thermoviscoplastic,  uniaxial  bar  can 


be  written  in  the  form 

3  {  EA  [3u  -  a.  -  a(T-T0)] 

3x  l  lx  1  R 


(1.1) 


A  [e(3u  -  a,  +  aT„)  3al  +  Ea2T  3T 
1  3x  1  K  3t  3t 


(1.2) 


-  AEaT  3  u  -  Ape  3T  +  A  3_  (k  3T) 
3x3 t  v  3t  3x  3x 


with  the  addition  of  internal  state  variables  a^,  a  ,  . ...  ap  satisfying 

i  =  n.(e,T,a.)  j=l,  n  i=l,  ....  n. 

3t  1  3 

The  strain  e  is  defined  as  e=  3u.  We  also  define  the  stress  a  by 

3x 


o  -  E(e-  al  -  a(T-TR)) , 


For  the  Bodner  model,  we  have  two  internal  state  variables,  a^,  the  axial 
inelastic  strain,  and  c^,  the  drag  stress.  By  assumption  these  satisfy 


*vf-Do  -a-  f-  s±i  f  i 


2n  \  a 


(1.3) 


=  f2(o  2,a)  =  m(zra2)  o^l-AZj  (a2~ZI  ) 

3t 


(1.4) 


where  0Q,  n,  m,  zp  A,  and  r  are  experimentally  determined  parameters 
(constants  ). 
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From  the  definition  of  stress,  we  may  write  the  stress  strain  relation 
as 

e  =  E'1  o  +  oj  +  cl(T-Tr).  (1.5) 

Another  relation  which  easily  follows  is 

E(3u  -  a+aT0)  =  a  +  ,EaT  (1.6) 

3x  R 

This  will  be  used  later  to  simplify  the  Fourier  heat  conduction  equation  (1.2). 
If  we  assume  that  px  =  0,  then  we  can  write  equation  (1.1)  as 

3_  {Act}  =  0  (1.7) 

3x 

which,  for  constant  cross  sectional  area  A,  implies  that 

a  =  o(t).  (1.8) 

Consequently,  no  matter  how  complicated  the  relationship  between  a  and  t, 

(1.8)  will  hold  as  long  as  thereis  no  applied  load  px* 

Instead  of  solving  the  forward  problem,  that  is  determining  o^.o^.u 
and  T  from  (1.1)  -  (1.4),  we  will  regard  a  *  a(t)  as  being  known.  As  we 
will  see,  this  allows  simplification  of  (1.2).  Viewing  a  =a(t)  as  a 
known  function,  we  will  derive  expressions  for  anc*  ^  terms  a» 

which  yields  an  implicit  relationship  between  e  and  a  via  (1.5). 

First  we  derive  a  relation  between  and  d.  If  we  substitute  (1.3) 
into  (1.4),  a  first  order  differential  equation  results. 

£*mU1'a*)8 


^  '  (1.9) 


F2[a2,a] 
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Given  a  =  cr( t),  a2  is  uniquely  determined  by  the  initial  conditions 


a2(x,o)=c^°\x). 

If,  initially  a^\x)=  constant,  then  d2  is  a  function  only  of  the  variable 
t.  Therefore  the  integral  relation 


F2(o(s),a2(s))ds 


(1.10) 


defines  c*2  implicity  as  a  function  of  a.  We  will  occasionally  write 

a2  =  a2(t;a(t)) 

to  emphasize  the  dependence  of  a2  on  a. 

Once  a2  has  been  determined,  by  numerically  integrating  (1.9)  if 
necessary,  we  may  solve  for  a^.  Again,  assuming  that  the  initial  data  for 
dj  is  constant, 

oi^t)  =  +  J  f1(a2(s)  ,a(s))ds.  (1.11) 

0 

Therefore,  knowledge  of  a  =  a(t)  determines  the  internal  state  variables 
oij  and  d2  as  functions  of  time  uniquely.  Mathematically,  we  have  effectively 
decoupled  (1.3),  (1.4)  from  (1.1)  and  (1.2).  Since  a  is  a  function  only 
of  t,  (1.1)  is  automatically  satisfied.  In  this  formulation,  (1.1)  becomes 
a  compatibility  condition  between  the  stress  a,  and  the  displacement,  u. 

Substituting  (1.5)  and  (1.6)  into  (1.2),  we  obtain 

(a+EdT)  9al  +  Ea2T  9T  -  EdT  3  fa  +  d  +  d(T-TDf| 

3t  3T  I^E  1  RJ 

-  pc  3T  +  k  32T  =  0  (1.12) 

9t  3? 
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This  may  be  simplified,  after  cancelling  terms,  to  the  linear  equation 


-pc  £+k32T  +  o(t)  ^1  -oo'  (t)T  =  0 

at  7?  at 


which  we  will  write  as 


aT  -  e2  a2T  +  fctoUt) 

9t  17  L  pcv 


aa^t) 

at 


(1.13) 


p  _c 

where  e  =  k _  is  typically  of  the  order  of  10  in  nondimensional ized  units. 

pCy  A 

The  presence  of  e£«l  causes  a  boundary  layer  to  form,  near  the  ends  of  the 
uniaxial  rod.  Near  the  center  of  the  rod,  for  small  time,  the  temperature 
field  is  nearly  spatially  homogeneous  and  is  in  close  aareement  with  the 
solution  of  (1.13)  obtained  by  setting  e  =  0.  Eventually,  however,  the 
boundary  layer  reaches  the  center  and  influences  the  temperature  field.  The 
investigation  of  this  boundary  layer  phenomenon  wi  1 1  be  the  subject  of  a 
future  technical  report. 

The  boundary  conditions  we  consider  are  those  appropriate  for  a  symmetric 
temperature  distribution,  with  T  initially  equal  to  TR,  and  with  convection 
boundary  conditions  holding  at  the  ends  of  the  rod.  We  normalize  the  domain 
by  setting  the  length  =  2  units. 

T(x,0)  =  Tr 


Tx(l,t)  =  0 


(1.14) 


kTx(0,t)-e(T(0,t)-Tj  =  0 


By  a  simple  transformation,  we  make  the  boundary  conditions  homogeneous. 
Define  the  quantity 

T  =  T(x,t)  -  Tw. 


/’S' ' 
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T  satisfies 


3  T  -  32T  +  fag1  (t)  1  T  =  g(t)  3ai  -LaoVt! 

3t  3X2  [CCv  J  CCv  9t  pcv 


T(x,0)  =  Tr  -  T. 
Tx(l,t)  =  0 
bTx(0,t)  -  8T(0,t) 


(1.15) 


The  boundary  value  problem  (1.15)  may  be  further  simplified  by  introducing 


the  modified  temperature 

T*(xxt)  *  T^t)  exp(cta(t)) 

PCy 

Multiplying  (1.15)  by  the  factor  exp  (  ac(t)  ) 

PC., 


(1.16) 


(essentially  an  integrating  factor),  we  obtain 
3  T*  -  e2  32T*  =  fo(t)  3al  -  T„  aa1 


T  =  g(t)  3al  -  aa1 (t)  exp  aa(t) 
pCy  3t  pcy  pcv 


Assuming  a(0)  =  0,  we  obtain 


T*(x,0)  =  Tr-  T*, 
T*(l,t)=0 

kT*  (0,t)  -  8T*(0,t)  =  0. 


(1.17) 


We  set 


nc  oc  u 

F(t) :  *  [g(t)  -  T«  2SL(tjl  exp  fao(tj] 

I  PCV  3t  PCV  I  PCV  I 

which  is  a  uniquely  determined  function  of  t,  given  a  * 


(1.18) 
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We  have  therefore  derived  the  following  linear,  constant  coefficient, 
heat  equation,  satisfied  by  T* 


ft  T*  -  e2  |!T*  =  F(t)  (1.19) 

T*  (x,0)  =  Tr  -  T. 

Tx*  (l,t)  =  0 

kTx*(0,t)  -  6T*  (0,t)  =  0 


The  temperature,  T.may  be  recovered  by  computing 
T  =  T(x,t)  =  +  T  (x,t) 

-  T.  +  T*  (x,t)  exp  (-ao(t) )  . 


(1.20) 


The  solution  of  (1.19),  by  Laplace  transform  methods,  is  the  subject 
of  the  next  section. 

2.  Solving  for  the  Modified  Temperature  by  means  of  Laplace  Transforms. 


First  we  define  the  quantities 


T(x,s)  =  L[T*(x,t)]  = 


-st 

e  T*(x,t)dt 


and 


F(s)=L[F(t)]. 

Applying  the  Laplace  transform  L  to  both  sides  of  (1.19)  and  to  the 
boundary  conditions,  we  obtain  a  second  order,  linear,  ordinary  differential 
equation  for  T(x,s). 


■'.'TV*?  y*  !■  r  |.||>^  wf 


r 


r 


h 


tx(i,s)  =  o 
kTx(o,s)-0T(O,s)=O. 


A  particular  solution  of  (2.1)  is 
T  =  _I_(  F  (s)+TR-Taj) 


So  the  general  solution  of  (2.1)  may  be  written  in  the  form 

T(x,s)=  l.  ( ^(s)+T^-Tn,]+A(s)exp(-x  YT)+B(s)exp(+x  Vs) 
s  _ 


(2.2) 


for  s>o,  where  A(s)  and  B(s)  are  uniquely  determined  by  the  boundary  conditions. 
If  0^0,  the  expressions  for  A(s)  and  B(s)  are  rather  complicated: 


B(s)  =  j k  V?  (1-exp  2  Vs~  )-0(l+exp  2  Vs) 

e  e  e 

A(s)  =  B(s)  exp  (2  Vp. 

e 


-1 


>]  6  1  [F(s)+TR-TJ 


(2.3) 


Substituting  these  expressions  into  (2.2)  we  obtain 

T(x,s)  =  l  [F(s )+TR-Tojs]  +  B(s)  exp  F(2-x)  VT"]+B(s)  exp(+x  Vs) 
S  e  e 


Therefore 


T(x,s)  =  1  [FfsJ+Tn-.Tc]  +01  [F(s)+Tr-T°o] 


exp(2-x  V?)  +  exp  (x  YT) 

_ c _ e _ 

k  Ys(l-exp(2  Vs)  -  e(l+exp(2_Vs) 
£  e  e 


=  1  [F(s)  +  Tr-T„]  sl-0  [exp  (l^x  VT)  +  exp  (x-1  Vs)  ] 


c 


l 

j 


k  V?[exp(-  Vs)  -  exp(  YsT)]-g[exp(-  Vs")+exp(+  Vsl] 

ec  e  e  e 
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=  I  CF(s)  +  Tr-  TJ  |l-  3  .  cosh  [(1-x)  Vs]  J 


k  Vs~  sinh(_Vs)+  cosh  Vs~ 

e  e  e 


(2.4) 


We  will  find  it  convenient  to  define  the  kernel  function  G(x,t)  by 

G(x,s)=  1-3  .  cosh  [(1-x)  Vs!  [k  VT sinh  VT  +  0  cosh  VT  ]_1 

£  £  £  € 


(2.5) 


The  expression  (2.4)  reduces  to 

R 


T(x,s)  =  1  [F( s )  +  Td-  Tec]  G(x,s) 
s 


which  implies 

T*(x,t)=L_1[T(x,s)]  =  L_1|l  [F(s)+Tr-T„]^  *  G(x,t) 

=  G(x,t)*[  J*  F(x)dt  +  Tr-  T  ] 
o 


(2.6) 


(2.7) 


Therefore,  from  (1.20),  (2.7),  and  (1.18) 

T  «  T,  +  T*(x»t)  exp  (-w(t)) 

pc 

v  t 

=  To  +  exp(-gg(t)  G(x,t)*[  |F(x)dT  +  TR-  Tj 

Pcv  Jo 


Substituting  (1.3)  for  9al  and  cancelling  Tm,  we  finally  obtain  the  result 


■y  s 


f* 


T(x,t)  =  T  +  exp(-  gg(t) )  G(x,t)*  |tr  + 

pc.,  l 


f  ! a( t)  I  2  D  exp(-n+l  /a2^\  )  exp  (a a(T))  dx 

X^Tvr  0 


-  T„  exp  (aa(t) ) 

pc.. 


(2.8) 


=  T  [a2>a  ] 


where  is  given  by  solving  (1.4),  G(x,t)  is  given  by  (2.5),  and  o(t)  is  a 
prescribed  function  of  time. 

/v 

We  note,  that  in  the  case,  8=  0,  the  expression  for  G  simplifies  to. 


G(x,s)=l 


consequently,  for  0=  0. 


T(x,s)=  i  [F(s)+  Tr-TJ 


and  T*(x 


ft 


x)dx  +  T_-  X,,  .  This  implies  that 


T=Tot  +exp(-gg(t))  [Tr-  T« 
pc.. 


J  pcv  3t 


exp  (gg(T))  dx-T^  exp  (ao(t)  +  T„  ] 


exp(-ao(t))  Tr  + 

pc„ 


+  exp  (-aa(t) 

pc.. 


f  0<t)  3al 

J  »cv  “ 


exp  (cxo(t)) 

pc.. 
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r 


=  exp(-ag(t))  It  +  f  |q(x) |  2_  D  exp  (-n+l/a2  \2n).exp(ag(x))dx  exp 

pcv  l  J  Pcv  VI  0  /  pcv 


This  reduces  finally  to 


fc  2  n 

T=  exp(-qq(t)  |t  r+J*  lq(t)l  2Do  exp|~  -^n+l^a2jx  +  ao(-c)j  dr  j  (2.9) 

We  note  that  this  is  a  function  only  of  time,  and  also  that  the  integrand 
is  strictly  positive.  This  implies  that  the  bracketed  quantity  is  monotonically 
increasing  function  of  time.  This  implies  that  even  in  the  periodic  stress 


case,  after  one  cycle,  a  relative  temperature  difference  of 

o  <  f  Mill  exp  [-(n+1)  +  Sil l  3  <It 

Jq  p°v  V3”  2n  \oVcTj  pcv 


(nil)  h{x)Yn  +  2 shl 
2n  \o{tJ)  pcv 


(2.10) 


has  developed;  where  tpdenotes  the  period  of  c(t). 

A  numerical  investigation  of  (2.9)  and  (2.10)  with  a-  cr(t)=  aQ  sin  wt, 
will  be  done  in  a  subsequent  report. 

3.  Nonlinear  Compliance  in  terms  of  Convolutions: 

In  order  to  express  the  strain  as  a  function  of  stress,  we  utilize  (1.5), 
(1.11),  and  (2.8)  to  obtain: 


e  *  E_1a(t)  +  aj  [a(t)]  +  a[T[a2,a]-TR] 
=  E~*  a  +  j fj(a2,a(t)dt  + 


(3.1) 


+  a  )  T^  -  Tr  +  exp  (-  aa(t))  G(x,t)< 
1  PC.. 


Rr  +  f  igjili  i  exp  (-n+1  (a2^)2n  +  aa(x))  dx-L  exp  (aa(x)  )1  1 

L  j  pCv  VT  pcv  pcv  JJ 
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We  will  write  this  in  the  form: 


e  =  E 


ft 

*a(t)+C+  I  H,[a]dT  +a  exp  (-aa(t) ) 

A  pcv 


where  C  =  +  a|X  -  TR], 


Hi [o]  =  fi(a„;a(t)]  -  2-  D_  a  exp  (-n+l/a2^\  ) 

1  12  yj-  0W  TnV“oT“/ 

and 

H2W=Tr  ♦  exp  [-glp^l]2"  *  Mil] 


(3.2) 


-  Tot  exp  .[gg(t)1 

pc„ 


agCt.aCt)]  is  given  by  the  solution  of 


^2  «  f2(a2,a)  =  m  (Zj-ag)  a  ^1  -  I  *l(^l 


(1.4) 


GHSfx.tH'1  |  1-8  C0Sh  (1-x)  V|1 


k  VT~  sinh  (  Vs~ )  +  0  cosh  VT- 


from  (2.5). 


•••  -«y.  ~z+ 


Equation  (3.2)  gives  a  closed  form  expression  for  the  compliance  of 
the  Bodner  model  with  two  internal  state  variables.  We  note  that  the 
contribution  of  is  the  integral 

+  y*  Hj[a]dt  =  +  J'  H^[a(t),  ^(t.oU))]  dt 

0  0 

while  the  effect  of  the  temperature  and  spatial  variation  is  to  introduce  a 

convolution  in  time,  multiplied  by  the  factor  exp(-ac(t)) . 

pcv 

If  cr ' ( t )  =  constant,  then  we  have  a  pure  convolution  with  A  non  constant 
stress  rate  introduces  the  operator 

exp(-aa(t)) .G(x,t)* 

oc 
M  V 

rather  than 

G(x,t)* 

A  detailed  examination  of  the  asymptotic  behavior  of  (2.8)  and  (3.2), 
as  t-*>°,  for  the  cyclic  case 
o(t)-o  sin  tot 

will  appear  in  a  future  report. 

The  function  G  contains  the  diffusion  phenomena  and  thermomechanical 
coupling  within  it.  A  detailed  analysis  of  G=G(x,t)  and  asymptotic  re¬ 
presentations  of  the  solution  to  (1.1)-(1.4)  will  appear  in  a  future  report. 


4.  Conclusions. 


In  the  spatially  uniform  case,  that  is  when  insulated  boundary  con¬ 
ditions  hold  at  the  ends  of  the  rods,  instead  of  (1.19)  we  have 


It T*  * 

T*(o)  -  T^-T< 


(4.1) 


which  has  the  solution 


T*(t) 


-fn 

J  n 


F[a(x)]dT  +  Tr-T*  (4.2) 

which  we  denote  by  T*  (t).  We  therefore  arrive  at  the  important 

o.a.e 

conclusion  that 


T  *  (x,t)  -  G(x,t)*  T*  (t). 
p.d.e.  o.d.e 


(4.3) 


Therefore,  given  the  solution  of  the  problem  in  the  case  of 
insulated  boundary  conditions  and  known  stress  history  a(t),  the 
solution  of  the  problem  for  general  boundary  conditions  is  recovered  by 
a  convolution  in  time  with  an  appropriate  heat  kernel  G(x,t). 

The  effect  of  the  term  in  (1.13)  is  to  force  a  boundary  layer  to 
form  at  the  ends  of  the  rod.  Away  from  the  ends,  and  for  small  time,  the 
solution  is  very  nearly  spatially  independent.  [2], 
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ABSTRACT 

Under  the  assumption  that  the  stress  is  a  known  function  of  time,  the  equations  govern¬ 
ing  the  behavior  of  a  homogeneous,  thermoviscoplastic  uniaxial  bar  (Bodner  Model)  can  be 
reduced  to  a  single  linear  partial  differential  equation  for  the  temperature  in  terms  of  the 
stress:  This  can  be  solved  by  transform  methods. 

In  this  report,  an  explicit  series  representation  of  the  temperature  in  terms  of  the  stress  is 
obtained  by  seperation  of  variables  techniques.  With  the  assumption  of  periodic  stress,  bounds 
on  the  time  averaged  maximum  temperature  increase  are  obtained. 

In  order  to  describe  the  thermal  boundary  layers  near  the  ends  of  the  rod,  an  asymptotic 
model  is  developed  which  uses  a  spatially  homogeneous  solution  as  a  starting  point.  Finally, 
various  mathematical  details  regarding  the  Green’s  function  and  series  solution  are  investi¬ 
gated. 

I 


I 
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Introduction. 

In  [  1  ]  ,  we  showed  that  the  equations  describing  the  behavior  of  a  thermovisco¬ 
plastic  uniaxial  bar  (as  derived  in  [  2  ] ) 


d* 


yc 


(0.1) 


\e  {^-—oli+olT k  )— +£  a2T¥- 

1  a*  y  y 


^--Apcv  *L+A-2-Uc$L)=-Arp, 


-AEaT^--Apcv 

Sx&  y  3*  3* 


y*i  _  2  n  cr 

- —  r  Z)  o  T“ T®^P 

y  V3  ]cr[ 


«  +1 
2ft 


“2 

O’ 


2n 


(0.2) 


(0.3) 


and 


-m  (z  i—a2)<r^^-—A  xz  x 

y  y 


a2-z/ 


(0.4) 


with  prescribed  stress  ,  <r=<r(r  ),  could  be  reduced  to  a  single  linear,  constant  coefficient, 
partial  differential  equation 


±<h-€2*!L =£[o<f)]:=£(r)  (0.5) 

y 

9(x,0)=7,j?— r* 

0x(l,f)=O 

*0x(O,f)-00(O,t)=O 

for  the  modified  temperature 

9=Kx,t  )=(T -T  Jexp( -^2.)  (0.6) 

pcv 


where 


F[o<t  )]= 


o<f)  y*i  _  aT^aXt ) 

■p£  y~  pc3 


exp 


acKf ) 

Pc, 


(0.7) 
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and 

e2=_L-.  (0.8) 

PCv 

In  this  report,  we  solve  (0-5)  by  means  of  separation  of  variables  and  write  down 
an  explicit  power  series  solution  of  (0-5).  This  will  enable  us  to  examine  the  relative 
importance  of  each  term  in  the  series  expansion,  evaluate  various  asymptotic  expressions 
involving  the  solution,  follow  the  formation  of  boundary  layers  in  the  bar,  and  to  exam¬ 
ine  the  effect  of  modifying  the  choice  of  ait ).  Finally,  in  the  appendices,  we  examine 
the  integral  representation  of  (Kx  ,t )  by  means  of  the  Green’s  function  which  was  intro¬ 
duced  in  [  1  ] . 


Section  1.  Series  expansion  of  the  modified  temperature. 
We  write  the  solution  of  the  problem  (0-5)  in  series  form 

<X*,f  )s£M*  (* )' 

n  =1 


Substituting  this  into  (0-5),  we  have  : 


z 


n  *1 


8;&v,uwi,(fv;(x)|-F&) 


£«»(0)*„(*M*-7’oo 

n  *1 


£0n(rVB(l)=O 

n  =1 


£en(f)  k  \fr„  (0)— (0) 

n  =1 


We  solve  (1.1)  by  assuming  that 


(1.1) 

(U) 

(1.3) 

(1-4) 


9  'n  ( t  )+€2u>n20a  (t  )-dn  F  Ct ) 


Therefore, 


tiXx  )+Oi^n  (x  )=0 


£  9„  (t  )tl/n  (x  >-€20n  (f  >//„'(x  )= 


£  0„  (r  (x  )+€2w„20n  (r  tyn  (x ) 


=  £  *Pn  (x )  0„  (r  )+e2w  20„  (f ) 

n  =1 


=  Z  +nU)dnFit) 

n  =1 


=  T.d’i'l>n(x)  F(.t) 

n  =1 


This  will  satisfy  ( 1.1)  if  and  only  if 

£d„</rn(x)= 1  (1.7) 

n  =1 

on  the  interval  0<x  <1.  Equations  (1.3)  ,  (1.4)  ,  and  (1.6)  determine  cu„  uniquely. 


Writing 


we  have 


*pn  (x  )=<*„  sin(l-x  )w„  +bn  cos(l-x  )o>n 


0=^„(l)=-a„  a>„ . 


This  implies  that  a„  =0.  The  second  condition  implies 


0=k  >pn  (Q)—fi\j/n  (0)=-6„  |*  <an  sin<o„  +0cosun 


which  implies  that 


Asymptotically , 


o>„  ->n  v. 

We  note  that  u>n  is  nonzero,  so  that  </Kx)=l  is  not  an  eigenfunction.  Therefore, 
(1.7)  determines  the  coefficients  dn  by 


cos(l— x  )o>n  =1 

n  =1 


(1.9) 


The  eigenfunctions  < iin  (x)=cor(l— x  )oin  can  be  shown  to  be  a  complete  orthogonal  basis, 
consequently 


dn 


f  cos(l— x  )o>n  dx 
i>  __  4sino>„ 

T  ~  "  2o>„  — sin2to, 

J  cosW—x  )oi„  dx 


(1.10) 


Asymptotically  , 


dn 


2sin<on 

<»n 


_  2g 

kn2‘rr2 


It  remains  to  satisfy  (1.2). 


Tr  -rM=  £  9n  (0Mn  (X  )=  £  0n  (0)COS(1-X  )wn 

n  =1  n  =1 

Comparing  this  expression  to  (1.9)  we  see  that  it  suffices  to  choose 

9n(o)=d„(TR-Tj.  (1.11) 

The  remaining  differential  equation 

9„(t)+€2<o2=dnF(t) 

9n(0)=d„  (TX~TJ 


has  the  solution : 


9„  (t  )-dn  (Tr  —T  00>xp(— o>2€2r )+ f  dn  exp(— €2w  2(r  —s  ))F  is  )ds 

o 


-  7  - 


•VW1 


Therefore  the  solution  of  (0.5)  can  be  written  down  concisely  as 


CO  ‘ 

0(x ,  r  )=  £  dn  (Ts—T oJexp (-ta „2€2f  )+dn  f  exp(— €2o>2(t  — j  ))F (r  )dr  co^l-x  >an  (l.l 

71  =1  0 

=(7*  — Too)  £  dn  exp(— to )cos(l—  x  )o)„  +  £  exp(-€2ca2r )  *  F(t )  co^l— x  )o>„ 


("  *  "  signifying  convolution  in  time  ) ,  to„  satisfies 


j(On  =~COt  , 


and  dn  satisfies  (1.10). 


Section  2.  Contribution  of  Series  Terms. 

We  will  examine  two  cases  ,  Fit  )—t  and  Fit  )=1,  in  order  to  examine  the  decay 
of  the  terms  9n  it ). 

Case  1  ,  Fit  )—t .  We  have,  after  some  computations  which  we  omit, 

9n  it  t  +dn  Tr  exp i-ufrh  )  (2.1) 

e2<  e2<  e4w„ 

This  implies  that  each  mode  contributes  a  term  asymptotically  of  order 


and  therefore  each  term  is  of  the  same  order  in  t . 

In  the  second  case.  Fit  )=1,  we  have  (again  after  some  easy  computations) 

Bn  it  )=  1— exp  (-€2o>  2r )  |+d„  iTg  -T  Jexp(-€2ca2r ) 


and  therefore,  as  t  -*co,  each  term  contributes 


As  t  -*oo  we  have  the  following  result 


(2.3) 


<Kx,f)-*£  )tl>B  . 

»*i 

At  the  right  endpoint  then,  the  temperature  is  given  asymptotically  by 

0Jl)=Z4%  (2.4) 

»  =1  €2< 

which  is  bounded. 

Section  3.  Upper  bounds  on  the  temperature  increase. 

Equation  (0.5)  can  be  considered  as  an  ordinary  heat  equation  describing  the  evolu¬ 
tion  of  a  symmetrical  temperature  distribution  in  a  rod,  with  convection  boundary  con¬ 
ditions  at  x  =0  and  x  =2,  starting  from  a  constant  temperature,  with  a  source  strength 
Fit)  which  is  spatially  homogeneous.  From  this  analogy,  or  using  maximum  principles 
of  parabolic  equations,  one  can  see  that  the  maximum  temperature  is  acuieved  at  the 
center  of  the  rod  (for  Fit  )^0)  and  that  if  Fit  )^F  then  the  temperature  distribu¬ 
tion  0max(x  J  )  satisfying  (0.5)  with  Fit )  replaced  by  Fmlx 

—d—€2  = F  ma*  -constant  (3.1) 

Oix,Q)=TR  -T^ 

0i(l,r)=O 

k  9X  (o,  t  y-poio,  t  )=o 

satisfies  the  relation 

(Kx.rXflaaxGc.r ). 

The  steady  state  solution  of  (XI)  is  very  easy  to  compute.  It  is 

0m«(*  )=€-2(-4*2+*  +4^  max 

L  p 

At  the  center  point,  where  the  maximum  temperature  occurs,  we  have 

Ul  ')=e~2ij+j')Fm» 

This  then  bounds  the  solution  of  (0.5)  at  the  point  x  =1,  and  therefore  everywhere. 


(12) 


(3.3) 


(3.4) 


Section  4.  Computing  F  ^  given  ait ). 


If  \<r(t  )|<Ci  and  \crit  )J  < C 2  then  we  can  estimate  F^  as  follows.  In  the  appen¬ 
dix  it  is  shown  that  if  zj  <a-^0)<zl  then  zj  <a2 Ct  )<2t  for  all  t.  This  then  yields  the 
inequality 


< 


2 

73 


D  oexp 


_  n  +1 

2l 

2 n 

Ci 

(4.1) 


and  therefore 


F(t)^Fmtx= 


C  i  aJ  ooC  2 

- L-  3+ - 

PCv  PCV 


exp 


Q(Ci 

pcv 


(4.2) 


This  leads  to  a  tremendous  overestimate  in  general  Since  the  series  solution  (1.12)  is 
analytically  correct  (Le.  not  an  approximation  )  one  can  obtain  the  temperature  at  the 
midpoint  by  substituting  x-l  into  (1.12)  to  obtain 


0(1,  t  )=  £  dn  exp(-o >„2€2r  ) 

n  =1 


T R  ’~Tao+ f  exp(e2o>  2s  )F  is  )ds 


(4.3) 


oo  f 

-*  £  f  exp(— €2o>2(t  — j  ))F  (s  )ds  (4.4) 

n  =1  0 


If  F  (s )  does  not  decay  as  s  -*oo  then  many  terms  will  have  to  be  included  in  (4.4). 
Examples  of  various  choices  of  cr—ait )  will  be  investigated  in  a  subsequent  report. 

Section  5.  Computing  the  time-averaged  mean  (modified)  temperature. 

Since  the  right  hand  side  of  (0.5)  depends  only  on  t,  through  ait ),  we  can  derive  a 
better  asymptotic  bound  on  the  behavior  of  <Xx  jt )  with  respect  to  time  in  the  following 
way.  First,  integrate  (0.5)  with  respect  to  time,  and  then  divide  by  t.  If  we  define 

/)=r_1  /(Kx.tOrfr  (5.1) 

o 


and 


<F(f  F{ai.r))d  r 


we  have  the  following  equation 


t~K(Kx,t )— fKx.O))— €2<fA.x,t  )-<Fit  )> 


<Ml,f)=0 

k<f>x(0,t)-p<tt0,t)=0 

Since  (Kx ,  t )  is  bounded,  t  ~K(Kx ,  t  )—Qix ,  0))-*0  regardless  of  initial  conditions. 
Therefore  the  quantity 

t 

4>(x  )=limr-1  fdix,  r)d  r  (5.2) 

o 

satisfies  the  ordinary  differential  equation 

-€2<I>'  U)=<F>  (5.3) 

<KD=o 

k  ViO)-p<PiO)=0 

where  <F  >-  lim  <F(t)>,  This  can  be  solved  to  yield 

T  —  oo 

4>(x  )=€-2(— -Lx2+x  +  -^-)  (5.4) 

2  p 

'Kx  )  is  the  asymptotic  state  of  the  (time  averaged)  mean  (modified)  temperature.  In  par¬ 
ticular  at  the  center  point ,  x  =1, 

<Kl )=e~Kj+j)<F  >  (5-5) 

If  the  exponential  term  in  (0.6)  is  approximately  unity,  then  <X<x )  is  the  asymp¬ 
totic  mean  temperature  change.  Therefore, 

<r(x)>~7’0O+«Kx) 

In  order  to  compute  <  F  >,  we  note  that  if  F  it )  is  periodic  then 

<F  >  =  \im>*rfF(t)dt=±rfF(t)dt 

f 

where  P  is  the  period  of  F .  If,  in  addition,  J  Fit  \it  =0,  then  < F  >  =0  and 


consequently,  there  is  no  net  mean  contribution  to  the  temperature.  Consider  the  second 
term  of  (0.7).  It  is  a  total  derivative  of  the  quantity 


Therefore  , 


T  ooexp 


aait) 


PCv 


•aTc 


*  PCV  PCy 


=  Taotx  p 


aait ) 


Pcv 


t=p 


=0. 


t  =o 


If  a2  approaches  a  constant  then. 


<F  > 


=rlnnf/ 

T—  oo  i 


1  r<Kr)  3<*i 

{  t>cv  IT 


dt 


(5.8) 


rt  +1 ,  “2 


efr 
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Given  cKr  )  periodic,  one  can  compute  (5.8)  and  then  substitute  into  (5_5)  to  compute  the 
asymptotic  mean  (modified)  temperature  change  at  x  =1.  If  a2 it)  does  not  approach  a 
constant,  then  one  cannot  use  equation  (5.8).  In  many  applications  however,  the  value  of 
a2  approaches  a  constant.  This  is  an  effect  of  the  phenomenon  of  "saturation".  The 
parameter  n  is  seen  to  be  crucial  in  predicting  the  temperature  rise  in  the  purely  one¬ 
dimensional  model. 


i 


In  general,  the  temperature  rise  predicted  by  substituting  (5.8)  into  (3.4)  is  much 
larger  than  physically  realized.  This  is  a  result  of  the  simplification  of  the  problem  into 
one-dimension.  In  reality,  we  cannot  ignore  radial  effects,  and  the  radiation  boundary 
conditions  along  the  longitudinal  surface  of  the  bar.  (  Neglecting  these  is  equivalent  to 
assuming  that  the  rod  is  insulated  along  its  length  ,  which  accounts  for  the  much  larger 
predicted  temperature  increase. ) 


The  one-dimensional  analysis  is  still  important  for  a  study  of  the  various  models 
for  thermoviscoplastic  materials,  and  the  effects  of  various  parameters  on  the  maximum 
temperature  changes.  For  short  times  an  alternate  representation  of  the  solution  can  be 
obtained  by  singular  perturbation  techniques.  This  is  analyzed  in  the  next  section. 


•/■.•vVv'r 


Section  6.  Boundary  Layer  Model. 


In  order  to  examine  the  behavior  of  the  boundary  layer,  we  expand  the  function 
( Kx,t )  in  terms  of  the  small  parameter  e.  We  define  the  stretched  coordinate  £=e_1x . 
This  implies  that  3,  =e_13^  and  that  therefore 

0,-%=F(r)  (6.1) 

<Xe£,  0)=Tr  ~T  m 

e-10£(e_1,r)= 0 

*  €_10^(O,  r ) — /30(o,  t  )=0 

In  the  usual  way,  we  write 

0(£r)=  «„(*.*) 

n  =0 

For  n=0  this  reduces  to 


0o,  *  ~ 0o,  H~F  (t )  ( 6.2) 

0O(. ,  0)=T  *  — F  <*, 

00,^00,0=0 

0O,  ^(0,  r  )=0 

which  implies  that  0o=0o(f ),  and  that 

t 

0o=0o(f  )= JF  (rV  t+F*  (6.3) 

o 

Equation  (6.3)  describes  the  evolution  of  the  modified  temperature  9ix ,  t )  under  insu¬ 
lated  boundary  conditions.  For  n  >  1  ,  we  obtain  a  coupled  set  of  linear  differential  equa¬ 
tions 


-  13- 


We  also  obtain  the  following  important  result,  when  x-0,  that 

0(0/  )— 0o(r  )— €y  —r)d  T  (6.7) 

which  to  leading  order  computes  the  temperature  as  a  function  of  time  of  the  left  end¬ 
point  of  the  bar.  The  first  order  correction  is  in  fact  monotonic  in  x  (if  0o(r)>O  )  ,  as 
expected.  Equation  (6.7)  shows  that  the  first  order  term  overcorrects  and  that  the  next 
term  (of  opposite  sign)  attempts  to  compensate  for  this. 


Appendix  I .  Estimating 


By  examining  (0.3)  and  (0.4)  with  D0  ,  m  ,  2 ,  ,  zj  ,  Ax  ,  and  r  positive  constants, 
and  noting  that  the  product  of  cr  and  is  intrinsically  positive,  we  have  the 
following  conclusions  :  If  initially 

zi  <  <*2(0)  <  Z  j 


dot  2 

then  as  one  approaches  z  t  from  below, - approaches 

92 


—AiZ  ! 


2 1  l 
2  1 


<0 


and  if  one  aproaches  Zj  from  above,  then  - approaches 

92 

m(z Zj  )<r^-> 0 

Consequently,  a£t  )  can  never  cross  these  two  limit  lines. 


Appendix  II .  Integral  representation  of  90c ,  t )  by  means  of  a  Green’s  function. 
In  [  1  1  we  derived  the  following  expression  for  the  Laplace  transform  of  0 


0(x ,  s  )=I  [/  Cr  )+7*  -T„  |  (l-/3G  (x ,  r ) 


where 


G(x,s  )=- 


cosh  |(l— x 


,  Vj  .  ,  Vr a  7  Vj 

k - sinh - — pcosh - 

e  €  6 


We  can  derive  a  series  expression  for  the  inverse  transform  by  means  of  a  contour 
integral 


G  (x ,  t  )= 


e”cosh 

€ 

2vi  i  TTs  .  ,  Vs  ~a  T7T 

c  1  00  It _ ©inh _ — flrncn _ 


ds. 


I 

s 

N 

1 


The  poles  of  the  integrand  are  the  roots  of  the  expression 


k  -^-sinh  ^-~l 3cosh 
€  €  6 


V7 

Setting  ——=ioi,  we  have 


0 =ki  osinh(i  to)— /3cosh(i  co) 


=—2k  usino>— 2/koso> 


Therefore 


k  coso) 

—o>=— — - 

p  sina» 


=-cor  o>„ 


(A2.1) 


which  is  the  same  equation  as  (1.8).  Therefore  s  =— €2o>2  where  oj  satisfies  (A2.1). 

Because  the  integrand  is  even,  there  is  no  contribution  of  the  branch  line  to  the 
contour  integral.  Similarly,  there  is  no  contribution  of  the  branch  point  at  s  =0.  The 
contour  at  infinity  goes  to  zero  as  s  —oo  if  Res  <0.  Consequently,  by  the  residue 
theorem. 


G  (x  /  )=2 '^residues  =  cosh(  1—  x  ) 


y/sn 


n  —0 


=  £  exp (— €2w  2r  )cos(l— X  )d)n 

n  =1 


and 


<Kx/)= 


f  F(.r')d  t+T  r—T  o 


-/3G(x/> 


+7  R  — 7’a0+ f  F  (rXf  r 


=0o(f)-0G(x/)*  0o(r) 


=0o-/3  Z  cos(l-x  )a>n  exp(— €2o>2r )  *  0o(r ). 


n  =1 


This  identity  motivates  the  series  expression  for  the  function  0(x  / )  given  in  sec- 
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ABSTRACT 


Under  the  assumption  that  the  stress  is  a  known  function  of  time,  the  equations 
governing  the  response  of  a  homogeneous,  thermoviscoplastic  uniaxial  bar  can  be  reduced 
to  a  single  linear  partial  differential  equation  for  the  temperature  in  terms  of  the  stress. 
In  this,  the  last  of  a  three  part  study  of  this  problem,  upper  and  lower  bounds  are  deter¬ 
mined  for  the  case  of  a  bar  subjected  to  radiation  boundary  conditions  on  its  longitudinal 
surface  and  in  the  presence  of  a  cyclic  stress  history.  The  temperature  rise,  which  is 
caused  by  inelastic  conversion  of  strain  energy  to  heat,  is  found  to  be  significant  for  some 
example  applications  to  realistic  materials. 


INTRODUCTION 


Large  space  structures  will  require  significant  passive  damping  in  order  to  sustain 
structural  vibrations  in  a  microgravity  field.  One  passive  damping  mechanism  which  has 
been  proposed  is  material  inelasticity.  However,  during  this  process,  a  substantial  propor¬ 
tion  of  the  strain  energy  in  the  structure  is  converted  into  heat  via  hysteretic  loss.  The 
purpose  of  the  current  three  part  research  effort  is  to  determine  if  this  energy  conversion 
process  can  produce  temperature  rises  which  are  large  enough  to  adversely  affect  the 
structural  integrity  of  the  system. 

In  Part  I  it  was  shown  that  assuming  the  stress  to  be  a  known  function  of  time 
decoupled  the  equations  governing  the  internal  state  variables  from  the  heat  conduction 
equation.  This  reduces  the  nonlinear  system  of  partial  differential  equations  governing 
the  motion  of  the  rod  to  a  single  parabolic  equation  for  the  temperature  with  coefficients 
given  in  terms  of  the  stress  ait).  We  then  obtained  an  expression  for  the  strain  e  in 
terms  of  a  convolution  of  the  stress  a  and  a  suitable  kernel  which  was  derived  from  the 
heat  conduction  equation.  This  yielded  the  compliance  relation  for  this  materiaL 

In  Part  n  it  was  noted  that  the  system  could  also  be  solved  by  senes  methods,  and 
that  a  boundary  layer  formed  near  the  ends  of  the  rod.  Series  and  integral  expressions 
for  the  temperature  were  derived  from  a  simple,  spatially  homogeneous  solution.  The 
boundary  conditions  examined  did  not  allow  for  the  radiation  of  thermal  energy  away 
from  the  longitudinal  surface  of  the  rod  ,  so  the  temperature  bounds  derived  in  part  II 
were  not  optimal 

In  Part  HI  we  develop  much  sharper  upper  and  lower  bounds  for  the  asymptotic 
temperature  rise  as  a  function  of  applied  loading,  material  properties,  and  geometry  for  a 
typical  space  structural  member.  Graphical  results  are  given  for  several  example  cases. 

Section  1.  Derivation  of  the  equations  of  Equilibrium.  Constitutive  Laws. 

The  equations  we  use  to  describe  the  quasi-static  response  of  a  thermoviscoplastic 
uniaxial  bar  have  been  formulated  and  discussed  in  [  1  ] .  For  brevity,  we  summarize 
the  main  points  in  this  section.  The  principal  equations  are  conservation  of  momentum 

CTjij=0  (1.1) 

and  the  conservation  of  energy  ,  as  expressed  in  the  modified  Fourier  heat  conduction 


A  jki  (*id  “«u/  +<**;  T  s  )“u  J  + A  jki  (<*,  J  <*u  TT  )— A  jti  (ai  j  T  eu)-pcvT-qjj= 0.  (1.2) 

(Dots  indicating  differentiation  with  respect  to  time.) 

In  addition,  we  have  the  stress-strain  relationship 

O' i  j  “A  jki  (*kl  ~*ti  ~*L )  (1*3) 

The  internal  state  variables  o$  and  e;  obey  the  following  evolution  equations  : 

€i ';=*(€«  (1.4) 

«?,=f?Mki  .T.ai,)  (1-5) 

In  the  case  of  a  homogeneous,  thermoviscoplastic,  uniaxial  rod  we  further  postulate  that 

O’,;  =CT11=(T  (1.6) 

ei)  =€ll~e 

A  jki  =®mi=^ 

=<*1 

=«11=« 

=°Ar  —Tr  ) 

and  finally 

qj=—kT'j  (1.7) 

£  is  a  constant  (  Young’s  modulus  ),  and  T*  is  a  reference  temperature  for  which 
there  is  no  thermal  strain  under  zero  applied  load.  We  also  assume  for  simplicity  that  a, 
k,  part  constant. 

Substituting  the  relations  (1.6)  and  (1.7)  into  (1.2)  we  have  the  following  system  of 
equations 


E  («— a1+a7* R  }ai+E  ae^TT  —E  aT  k—pcv  T  +k  A T  =0 
a9=fIJ(ap  ,T  ,e) 

<T=E[e—a1—a(T  —T  R )] 


(1.9) 


(1.10) 

(1.11) 


and 


€= 


(1.12) 


For  simplicity  we  have  assumed  infinitesimal  deformations;  (1.12).  Since  u  and  tr 
are  now  functions  of  the  axial  coordinate,  x ,  alone,  we  drop  the  tensor  subscript  notation 
and  use  subscripts  to  label  the  different  internal  state  variables  (1.10). 

Given  initial  conditions  for  a q  ,T  ,  and  u,  as  well  as  boundary  conditions  for  T 
and  u ,  this  is  a  well-posed  mathematical  problem.  It  has  been  shown  in  [  2  ]  to  be  ther¬ 
modynamically  consistent  as  well. 

We  make  the  following  observations  at  this  time.  First,  (1.8)  implies  that  tr  is  a 
function  of  time  alone.  Secondly,  in  the  special  case  where  the  f9  do  not  explicitly 
depend  on  T  or  €  but  only  on  {  ,  <r  }  ,  then  given  aitX  (1.10)  is  a  closed  system  of 

equations  ,  and  we  have  ~a9  (t  ;oKt )).  Equation  (1.8)  is  automatically  satisfied  for 
< r-ait )  so  (1.9)  essentially  becomes  a  single,  non-linear  parabolic  differential  equation 
for  the  temperature  T  in  terms  of  cr.  This  analysis  has  been  carried  out  in  [  3 , 4  ]  . 
Generalizing  this  approach,  we  will  obtain  the  asymptotic  behavior  of  T  directly  from 
the  system  (1.8X1.12). 


Section  2.  Radiative  and  convective  boundary  conditions. 

In  order  to  include  the  thermal  radiation  boundary  conditions  along  the  lateral 
surface  of  the  rod,  we  integrate  (1.9)  over  the  (constant)  cross-sectional  area  of  the  rod. 
Let  ( x  ,  r  ,  9  )  denote  cylindrical  coordinates  with  x  measuring  the  distance  along  the 
axis  of  the  rod.  In  this  model,  all  variables  (except  for  T)  are  independent  of  r  and  9. 


AE  (e—a  t+aT R  )dj+ AE  o?TT  —AE  aT  e— A  pcv  T  +kATxx 


(2.1) 


where  qn  is  the  integrated  effect  of  the  normal  component  of  the  thermal  flux  along  the 
lateral  surface  of  the  rod,  and  C  is  the  circumference  of  the  rod-  T  is  now  interpreted 
as  an  integral  average  of  the  temperature  over  the  cross-section  of  the  rod. 

It  is  clear  that  the  term  +Cqn  acts  as  a  "sink."  for  thermal  energy.  In  deep  space, 
we  model  the  thermal  flux  due  to  radiation  as 


qn  -O’,  8CT  4-T£)-Q  =cr,  ST  4-Q  (23.) 

where  T D  =<T  K  is  the  deep  space  ambient  temperature,  and  crs  =5.775x  1CT11  is 
Boltzmann’s  constant.  S  is  an  order  one  constant,  called  the  emissivity,  which  measures 
the  effect  of  the  surface  coating  on  black  body  radiation,  and  Q  measures  the  solar,  earth, 
and  deep  space  thermal  flux  incident  on  the  rod. 

We  approximate  (2.2)  by  the  first  two  terms  of  its  Taylor  series  expansion  in  order 
to  obtain  the  linearized  convective  boundary  condition 

qn  =o-s  S T£+fi(T  -Tr  )-Q  (2.3) 


with 


3-4 <rf  87/ 

In  order  to  estimate  the  actual  temperature  rise  (with  non-linear  boundary  conditions) 
we  choose  3  so  that 

tr,  ST/+3(r  -Tr  )  >  <rs  ST 4  >  <r,  ST^+4cs  STrKT  -Tr  )  (2.4) 

over  the  region  of  interest.  This  implies  that  any  temperature  change  associated  with 
(2.3)  will  bound  from  below  the  actual  temperature  change  associated  with  (2.2).  Choos¬ 
ing  3=4 <r,  ST r  will  bound  the  temperature  increase  from  above,  via  inequality  (2.4). 

Section  3.  Further  reduction  of  the  equations. 
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We  can  further  simplify  (1.9)  by  noting  that  the  left  hand  side  of  (2.1)  can  be 
written  as 

A  Ifa+E ctT  )at+E  o?TT  - olT (&+E  at+E  at  )-pcv  T  +kTxx  )  (3.1) 

“A  j— pcv  T  +kTxx  —or aT  +o-a1  j 

Consequently, 


— pcvr  +*ru  -a&T  +<rat1=+^qn 


S+£[trf8T/+/3(7--7JfH2] 


We  write  this  as 


pc,  T  -*r„  +adT  +££(7  -Tr  |tr,  STj-Q^ 


Letting  (Kx  j.  )=exp(— —Tg  ) ,  we  have 

Pc, 


i!® + -£*  q=f  a ) 
a*  pc,  a*7  pc,  a 


where 


ext*"*0 

<**)*,&  ) 

.  ^ 57,„4+  ^  j 

**  PC, 

P^v 

Ape,  *  Ape,  pc. 

The  modified  temperature  9  has  been  introduced  previously  in  [  4  J.  If  we  integrate 
both  sides  of  (3A)  with  respect  to  t,  and  divide  by  t,  we  obtain 

t-K0(xt)-0(x,O)>— ■^-,ljr0Gc.T)rfTL  +_££_,-»/«*  , r)d r*f ~lfF (r)d r 
PC,  Jo  PC,  A  J0  J0 


Am  t—oo  ,  with  9  bounded,  this  reduces  to 


and 


f 

<9  >(*  )=lim?  -1  /fl(x  ,r)d  r 

I  — °o  •{) 


(15) 


<F>. Ii».-/«xpcf2«)  « -gZLgjL-A-,  L 

*—  *&  pc,  pc,  ApCy  Apcv  pcv 


(16) 


= lim*  _1  j*exp( 

t  a  Acv 


<Kr»,(r)+  Qg__CfLw.  I 

pc,  A  pcv  Apcv 


are  the  asymptotic  mean  values  of  9  and  F  respectively. 

In  the  materials  under  consideration,  <ra1>0  .  Since  the  integrands  are  bounded, 
the  above  limits  exist.  For  simplicity,  consider  the  following  boundary  conditions  for  the 
temperature  T  :  T  (Oj  )=T K  and  TxtLj  )=0.  The  second  boundary  condition  results 
from  symmetry,  for  a  bar  of  length  2 L  units.  Therefore  9  satisfies  the  boundary  condi¬ 
tions  : 


<0>( 0)=0. 


(17) 


and 


<9>'U')*Q 


(18) 


The  solution  of  (14X  subject  to  the  boundary  conditions  (17,18)  is 


<0>(x)=^i</'> 


1- 


coshA  (L  —x  ) 
coshA 


(19) 


where 


A=(C9)i/2 

A/C 


Note  that  if  we  ignore  the  spatial  variation  of  the  temperature,  we  have 


<9>(x)»^l  <F> 


which  is  a  good  approximation  to  (19)  if  A  is  large.  This  accounts  for  the  close  approxi¬ 
mation  observed  between  the  spatially  homogeneous  solution  and  the  spatially  varying 


solution  with  boundary  layer  in  [  5  \ 


Since  p  ,cv  ,$  ,C ,  A ,  L ,  and  k  are  all  known  ,  the  spatial  dependence  of  <0>(x) 
is  determined.  Only  the  magnitude  is  a  function  of  the  loading  history  ,  through  <F  >. 
The  crucial  term  to  estimate  therefore  is  <F  >.  If  (jit)  is  periodic,  and  if  6tt  is  a 
periodic,  then  (16)  reduces  to 


Clio) 


where  P  is  the  period  of  a  and  oti. 

A  case  of  particular  interest  is  the  Bodner  model  [  7  ] ,  for  which 


At  saturation,  the  drag  stress  a2=a2  is  constant,  and  <F  >  becomes  a  function  only  of 
a. 


* 

<F  >=P~*f  exp 
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Since 


<  9  >Gc  )*  limr  SCx  ,r)d  r*  limf  ~1/ exp(  22^1 -XT  —TK)dT 

/-•  <  0  PCy 


we  can  recover  a  mean  temperature  T  by 


f  *<e>(x)|limr-,/expi!!HllldT 

K  rt  PC y 


'0  PCy 


hi 


+TS 


(112) 


For  stress  amplitudes  of  interest,  the  factor  exp 


ocait) 


P^v 


is  very  dose  to  unity. 


Fixing  the  material  parameters  {  a,p, c,  ,D0,n  ,$,kja2  I.  on*  can  compute  the  depen¬ 
dence  of  <9>  and  f  on  ait  X  With  a  periodic,  (112)  simplifies  as  before  to  an  integral 


over  the  period.  We  compare  the  analytical  results  with  experimental  and  numerical 
data  in  the  next  section. 

Section  4.  Numerical  Results. 

In  order  to  verify  the  accuracy  of  the  asymptotic  results  (3.9)  and  0\.li),  we  com¬ 
puted  the  average  temperature  for  the  parameters  given  in  [  5  ]  with  th  e  results  shown 
below  in  fig.l  .  The  values  for  the  emissivity  (with  thermal  equilibrium  at  295  deg. 
Kelvin)  are  0.24 , 0-5 ,  and  0.85  respectively.  The  cyclically  saturated  stress-strain  curve 

used  to  compute  <F  >  for  the  case  <rT..= 336-5  is  given  in  fig.  2.  The  values  for  o~ _ 

and  a2  were  taken  from  numerical  data  supplied  from  reference  [  5  ] . 

The  parameters  in  [  5  ]  describe  a  hollow  cylindrical  bar,  of  uniform  cross-section, 
orbiting  at  an  altitude  of  4,080  km-,  painted  with  a  high-emissivity  coating  (UTRE- 
S13GLO)  white  paint,  with  full  exposure  to  the  sun.  Under  zero  applied  stress;  equili¬ 
brium  of  thermal  flux  occurs  at  72  degrees  fahrenheit  (295  deg.  Kelvin)  for  which  the 
thermodynamic  parameters  of  the  metal  6061-T6  aluminum  were  experimentally  com¬ 
puted  [  6  I 

We  also  can  compute  the  behavior  of  the  thermal  boundary  layers  near  the  ends  of 
the  rod.  From  equation  (3.3)  or  (3*4)  we  compute  that  the  diffusion  constant  is  equal  to 
k  1,2 

( - )  .  This  is  typically  of  the  order  10~3  to  10~5  in  dimensionless  units.  This  leads 

pc* 

to  very  sharp  boundary  layers,  and  exponential  decay  to  the  asymptotic  temperature 
state. 

We  also  mention  that  other  constitutive  laws  than  (111)  have  been  derived  which 
describe  the  behavior  of  metal-matrix  materials.  Since  the  parameters  are  described  by  a 
parameter  fitting  procedure,  essentially  the  same  results  are  obtained  for  the  asymptotic 
temperature  rise. 

Sections.  Conclusions. 

The  results  given  in  reference  [  5  ]  indicate  that  the  cyclic  stress  amplitudes  for 
load  sequences  oflHz.,5Hz.,and25Hz.are  336J  MPa  (484  KSI) .  346.8  MPa  (50.3 
KSIX  and  358.9  MPa  (520)5  KSI)  respectively. 

From  Fig.  1  it  is  found  ,  for  example,  that  a  stress  amplitude  of  3464  MPa  leads  to 
a  minimum  asymptotic  temperature  rise  of  61.9®  K.  This  temperature  represents  an 
extremely  hostile  environment  for  AL  6061-T6  which  could  lead  to  such  serious  losses 
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in  material  properties  as  to  compromise  the  structural  integrity  of  the  system.  It  is 
therefore  concluded  that  the  use  of  material  inelasticiy  to  enforce  passive  structural 
damping  should  be  utilized  with  caution  due  to  the  possibility  of  catastrophic  structural 
heating  in  space  in  the  presence  of  solar  flux. 


MAXIMUM  STRESS 


TEMPERATURE  CHANGE 
(UPPER  AND  LOWER  BOUNDS) 


STRAIN 


STRESS  -  STRAIN 
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